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Abstract 

The work contained herein constitutes a report of the "Beyond the Stan- 
dard Model" working group for the Workshop "Physics at TeV Collid- 
ers", Les Houches, France, 26 May-6 June, 2003. The research pre- 
sented is original, and was performed specifically for the workshop. 
Tools for calculations in the minimal supersymmetric standard model 
are presented, including a comparison of the dark matter relic density 
predicted by public codes. Reconstruction of supersymmetric particle 
masses at the LHC and a future linear collider facility is examined. Less 
orthodox supersymmetric signals such as non-pointing photons and R- 
parity violating signals are studied. Features of extra dimensional mod- 
els are examined next, including measurement strategies for radions and 
Higgs', as well as the virtual effects of Kaluza Klein modes of gluons. 
An LHC search strategy for a heavy top found in many little Higgs 
model is presented and finally, there is an update on LHC Z' studies. 
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Parti 

Introduction 

B. C. Allanach 

The workshop took place at the ecole de physique in the lee of Mont Blanc, and lasted for 
two weeks. Computer systems were installed by the helpful LAPP staff for use by participants 
throughout the workshop. The first two days consisted of plenary talks intended to stimulate 
ideas and the proposition of projects. A couple of further plenary seminars on hot topics occured 
sporadically later in the workshop. The Beyond the Standard Model working group (convened 
by M Battaglia, M Nojiri, T Rizzo, A de Roeck, D Tovey, M Spiropoulu and the author) held a 
meeting on the second evening to set individual projects for small groups of interested parties. 
This report contains a summary of the fruits of participants' labour during and after the work- 
shop on those projects. The projects were phenomenological studies of both supersymmetric 
and non-supersymmetruc models and. The report first discusses the studies of supersymmetric 
models and then those of extra dimensions. We close with an update on Z' studies. 

At the time of the workshop the Tevatron and DESY collider runs were proceeding, and 
the start-up of the Large Hadron Collider (LHC) was eagerly awaited. Many hopes are concen- 
trating upon the production and detection of supersymmetric (SUSY) particles at these colliders. 
Although an immense amount of literature has been accumulated on SUSY phenomenology, 
there is still a large amount of work to do because of its complexity and the abundance of mod- 
els of SUSY breaking. The models predict different cascade chains leading to radically differing 
signals in experiments. In order to facilitate the phenomenological study of SUSY models, we 
need to both calculate the sparticle spectrum and also we need to simulate events. These two 
calculations are typically performed by seperate calculational tools. For each calculation, there 
exist several competing tools performing the same task with different approximations or as- 
sumptions. A common interface between the different tools has clear advantages, and much 
time was spent at (and after) the workshop arguing, debating and negotiating the various con- 
ventions. The write-up of this "SUSY Les Houches Accord" constitutes part|ll|of this report. 
A form-based web tool is presented in part |llll in which one can determine the spectra from 
mSUGRA models via the different public sparticle spectrum generators. The difference be- 
tween the predictions of the different generators gives an idea of the theoretical uncertainties 
involved in the calculation. Many recent works have used the WMAP determination of the 
cold dark matter density VLcdm to vastly restrict the minimal supersymmetric standard model 
(MSSM) parameter space. An initial study on the uncertainty induced from sparticle masses 
upon the prediction of the dark matter density is then entered in part|TVl Supposing supersym- 
metric particles are measured in colliders, a fit to the observables (such as the masses) will 
restrict the SUSY breaking parameter space, and discriminate models of SUSY breaking. A 
tool which enables one to perform this fitting efficiently is presented in part|V| A new code to 
determine the branching ratios and decays of SUSY particles is presented in partlVTl 

The report then turns to issues surrounding the measurement of sparticle mass and mixing 
parameters. The measurement of the lightest chargino mass in a universal minimal supergravity 
(mSUGRA) model is presented in part IVIIi It is pointed out in part IVIIll that without assump- 
tions about SUSY breaking, information coming from a future linear collider facility could be 
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extremely useful when analysing neutralino and chargino signals at the LHC. A bold new re- 
construction technique for SUSY processes at the LHC is proposed in part|K| Normally, it is 
not possible to reconstruct the neutralino momenta involved in R-parity conserving events since 
they remain undetected and the overall energy of the hard collision is not known. In the new 
technique, particularly long SUSY cascade chains are identified and lead to an over-constrained 
system when pairs of events are considered. The idea is pushed further to ensembles of more 
than 20 events in part|Xl If the idea stands up to further scrutiny, the new method would be the 
one yielding the most information about sparticle masses (provided the relevant decay chain(s) 
is(are) present). 

We next consider non-standard supersymmetric signatures, such as (part lXIl) non-pointing 
photons at the CERN LHC, which are often predicted in gauge mediated SUSY breaking mod- 
els. R-parity violation provides an oppurtunity to understand neutrino masses without the need 
for adding gauge singlets, and also to correlate neutrino oscillation observables with SUSY col- 
lider signatures. Part IXIII shows that the scenario predicts a relation between branching ratios 
of lightest SUSY particle (LSP) decay modes given the atmospheric neutrino mixing angle, 
providing a useful test. In part IXIIIl resonant slepton production at the LHC is examined in 
scenarios with ultra-light gravitinos. 

Models and signals incorporating extra dimensions are then considered. Many extra di- 
mension models predict an additional higgs-like scalar (the radion), which stabilizes the branes. 
In part IXIVl the mixing between radion and Higgs' in a two-Higgs doublet model is investi- 
gated. The discovery potential for two decay modes of the radion is determined in part IXVI 
Models with flat dimensions often predict Higgs decays into invisible graviscalars. These de- 
cays would provide a signal for the extra dimensions via an invisible Higgs width, and they are 
investigated in part IXVII Also, certain models have the lightest Higgs as a mixture of brane 
and bulk scalars. This unfortunately would suppress Tevatron Run n or 500 GeV linear collider 
Higgs signals, but would enhance production at the LHC or CLIC. Such issues are examined 
in part IX VlH Part IX Villi examines the sensitivity of the LHC for gluonic Kaluza- Klein states 
by their effects on dijet production. We next present an LHC search study of the heavy supple- 
mentary top quark present in many little Higgs models. An update on Z' studies at the LHC 
is presented in the final part lXXI The analysis focuses on the combination of several measure- 
ments in order to distinguish models. A modification of the leptonic ApB measurement proves 
to be very useful in this respect. 
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Part II 

The SUSY Les Houches Accord Project 

p. Skands, B.C. Allanach, H. Baer, C. Baldzs, G. Belanger, F. Boudjema, A. Djouadi, R. God- 
bole, J. Guasch, S. Heinemeyer, W. Kilian, J-L. Kneur, S. Kraml, F. Moortgat, S. Moretti, 
M. MUhlleitner, W. Porod, A. Pukhov, P. Richardson, S. Schumann, P. Slavich, M. Spira, G. Wei- 
glein 

Abstract 

An accord specifying a unique set of conventions for supersymmetric 
extensions of the Standard Model together with generic file structures 
for (1) supersymmetric model specifications and input parameters, (2) 
electroweak scale supersymmetric mass and coupling spectra, and (3) 
decay tables is defined, to provide a universal interface between spec- 
trum calculation programs, decay packages, and high energy physics 
event generators. 

1. INTRODUCTION 

An increasing number of advanced programs for the calculation of the supersymmetric (SUSY) 
mass and coupling spectrum are appearing [1-5] in step with the more and more refined ap- 
proaches which are taken in the literature. Furthermore, these programs are often interfaced to 
specialized decay packages, [6-8], relic density calculations [9, 10], and (parton-level) event 
generators [11-18], in themselves fields with a proliferation of philosophies and, consequen- 
tially, programs. 

At present, a small number of specialized interfaces exist between various codes. Such 
tailor-made interfaces are not easily generalized and are time-consuming to construct and test 
for each specific implementation. A universal interface would clearly be an advantage here. 
However, since the codes involved are not all written in the same programming language, the 
question naturally arises how to make such an interface work across languages. At this point, we 
deem an inter-language runtime linking solution too fragile to be set loose among the particle 
physics community. Instead, we advocate a less elegant but more robust solution, exchanging 
information between FORTRAN and C(-i-i-) codes via three ASCII files, one for model input, 
one for model input plus spectrum output, and one for model input plus spectrum output plus 
decay information. The detailed structure of these files is described in [19]. Briefly stated, the 
purpose of this Accord is thus the following: 

1 . To present a set of generic definitions for an input/output file structure which provides a 
universal framework for interfacing SUSY spectrum calculation programs. 

2. To present a generic file structure for the transfer of decay information between decay 
calculation packages and event generators. 

Note that different codes may have different implementations of how SUSY Les Houches Ac- 
cord (SLHA) input/output is technically achieved. The details of how to 'switch on' SLHA 
input/output with a particular program should be described in the manual of that program and 
are not covered here. 
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2. CONVENTIONS 

One aspect of supersymmetric calculations that has often given rise to confusion and consequent 
inconsistencies in the past is the multitude of ways in which the parameters can be, and are 
being, defined. Hoping to minimize both the extent and impact of such confusion, we have 
chosen to adopt one specific set of self-consistent conventions for the parameters appearing in 
this Accord. These conventions are described in the following subsections. As yet, we only 
consider R-parity and CP conserving scenarios, with the particle spectrum of the MSSM. 



2.1 STANDARD MODEL PARAMETERS 

In general, the SUSY spectrum calculations impose low-scale boundary conditions on the 
renormalization group equation (ROE) flows to ensure that the theory gives correct predic- 
tions for low-energy observables. Thus, experimental measurements of masses and coupling 
constants at the electroweak scale enter as inputs to the spectrum calculators. 

In this Accord, we choose a specific set of low-scale input parameters, letting the elec- 
troweak sector be fixed by 

1 . The conventional electromagnetic coupling at the Z pole, otcxJ^vriz)'- 

a 

Oiem{mz) = 7^^ , (1) 

1 - Aaicp(mz) - Aal^^^{mz) - Aatop{mz) 

where a is the fine structure constant, A«icp(m^) and Aatop{mz) represent the quan- 
tum corrections coming from leptons and top quarks, respectively (see [20,21]), and 
Aa^J^{mz) is the contribution from the five light quark flavours (see e.g. [22]). 

2. The Fermi constant determined from muon decay, Gp- 

3. The Z boson pole mass, m^. 

All other electroweak parameters, such as mw and sin^ 9w, should be derived from these inputs 
if needed. 

The strong interaction strength is fixed by as{mz)^^ (five-flavour), and the third gener- 
ation Yukawa couplings are obtained from the top and tau pole masses, and from mh{mh)^^, 
see [22]. The reason we take mb{mh)^^ rather than a pole mass definition is that the latter 
suffers from infra-red sensitivity problems, hence the former is the quantity which can be most 
accurately related to experimental measurements. If required, relations between running and 
pole quark masses may be found in [23, 24]. 

It is also important to note that all the parameters mentioned here should be the 'ordinary' 
ones obtained from SM fits, i.e. with no SUSY corrections included. The spectrum calculators 
themselves are then assumed to convert these parameters into ones appropriate to an MSSM 
framework. 

Finally, while we assume MS running quantities with the SM as the underlying theory as 
input, all running parameters in the output of the spectrum calculations are defined in the mod- 
ified dimensional reduction (DR) scheme [25-27], with different spectrum calculators possibly 
using different prescriptions for the underlying effective field content. More on this in section 
1231 
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2.2 SUPERSYMMETRIC PARAMETERS 

The chiral superfields of the MSSM have the following SU{3)c ® SU{2)l ® U{1)y quantum 
numbers 

L: (1,2, -i), E: (1,1,1), g:(3,2,i), f/:(3,l,-|), 

D : (3,1, 1), iJi: (1,2,-i), : (1,2,|). (2) 

Then, the superpotential (omitting RPV terms) is written as 

W = eab [{Ye\jHIL\E, + {YD)^,mQ\Dj + {Yu)^jH',Q'^U, - fiH'^H^] . (3) 

We denote SU{2)l fundamental representation indices by a, 6 = 1,2 and generation 
indices by z,j = 1,2,3. Colour indices are everywhere suppressed, is the antisymmetric 
tensor, with ei2 = e^^ = 1. Lastly, we will use t, 6, r to denote the i = j = 3 entries of mass or 
coupling matrices (top, bottom and tau). 

The Higgs vacuum expectation values (VEVs) are {H^) = Vi/^/2, and tan/9 = V2/V1. 
We also use the notation v = ^/vf + v\. Different choices of renormalization scheme and scale 
are possible for defining tan/3. For the input to the spectrum calculators, we adopt by default 
the commonly encountered definition 

tan/3(m^)^, (4) 

i.e. the tan /? appearing in the input is defined as a DR running parameter given at the scale m^. 
However, an option is included to allow tan /3 to be input at a different scale, tan/3(Minput 7^ 
m^)°^ . Lastly, the spectrum calculator may be instructed to write out one or several values of 
tan /3(Q)°^ at various scales Qj, see [19]. 

Finally, the MSSM DR gauge couplings are: g' (hypercharge gauge coupling in Standard 
Model normalization), g (SU {2)l gauge coupling) and g^ (QCD gauge coupling). 



2.3 SUSY BREAKING PARAMETERS 

We now tabulate the notation of the soft SUSY breaking parameters. The trilinear scalar inter- 
action potential is 



h.c. , (5) 



where fields with a tilde are the scalar components of the superfield with the identical capital 
letter. In the literature the T matrices are often decomposed as 



= Aij ; (no sum over , (6) 

J- in 



where Y are the Yukawa matrices and A the soft supersymmetry breaking trilinear couplings. 
The scalar bilinear SUSY breaking terms are contained in the potential 

Ui^{ml)ijU*^ + diHi^i)ijd*i^ + ei^{mj)iji*^ - {mleabH'^H^ + h.c.) . (7) 
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Instead of ml itself, we use the more convenient parameter niA, defined by: 

< = ■ (8) 

sm p cos p 

which is identical to the pseudoscalar Higgs mass at tree level in our conventions. 

Writing the bino as h, the unbroken SU{2)l gauginos as w"^^^'^'^, and the gluinos as 
~x=i...8^ the gaugino mass terms are contained in the Lagrangian 

£g = i (^Mibb + M2W^w^ + M^g^g^^ + h.c. . (9) 
2.4 MIXING MATRICES 

In the following, we describe in detail our conventions for neutralino, chargino, sfermion, and 
Higgs mixing. Essentially all SUSY spectrum calculators on the market today work with mass 
matrices which include higher-order corrections. Consequentially, a formal depencence on the 
renormalization scheme and scale, and on the external momenta appearing in the corrections, 
enters the definition of the corresponding mixing matrices. Since, at the moment, no consensus 
exists on the most convenient definition to use here, the mixing matrices should be thought 
of as 'best choice' solutions, at the discretion of each spectrum calculator. For example, one 
program may output on-shell parameters at vanishing external momenta in these blocks while 
another may be using DR definitions at certain 'characteristic' scales. For details on specific 
prescriptions, the manual of the particular spectrum calculator should be consulted. 

Nonetheless, for obtaining loop-improved tree-level results, these parameters can nor- 
mally be used as is. They can also be used for consistent cross section and decay width cal- 
culations at higher orders, but then the renormalization prescription employed by the spectrum 
calculator must match or be consistently matched to that of the intended higher order calcula- 
tion. 

Finally, different spectrum calculators may disagree on the overall sign of one or more 
rows in a mixing matrix, owing to different diagonalization algorithms. Such differences do 
not lead to inconsistencies, only the relative sign between entries on the same row is physically 
significant, for processes with interfering amplitudes. 

2.41 NEUTRALINO MIXING 

The Lagrangian contains the (symmetric) neutralino mass matrix as 

£-f- = -l^p'^'M^otp' + h.c. , (10) 

in the basis of 2-component spinors = {—ib, —iw^, hi, h2)^. We define the unitary 4 by 4 
neutralino mixing matrix A^, such that: 

- ^^°'^-^^oV^° = -}-^''^N^N*M^oN\N4>' , (11) 




where the (2-component) neutralinos Xi are defined such that their absolute masses increase 
with increasing i. Generically, the resulting mixing matrix N may yield complex entries in 
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the mass matrix, diag(m^o) j = m^oe'"^\ If so, we absorb the phase into the definition of the 
corresponding eigenvector, Xi* Xi'e*'^'/^, making the mass matrix strictly real: 

diag(m^o) = [N*M^oN^].^ = m^oSij. (12) 

Note, however, that a special case occurs when CP violation is absent and one or more of the 
rriyO turn out to be negative. In this case, we allow for maintaining a strictly real mixing matrix 
A^, instead writing the signed mass eigenvalues in the output. Thus, a negative m^o in the output 
implies that the physical field is obtained by the rotation x° — > x°e"^^/^. 

2.42 CHARGINO MIXING 

We make the identification = {w^ =f iw'^)/\/2 for the charged winos and h^, for the 
charged higgsinos. The Lagrangian contains the chargino mass matrix as 

^mass ^ _ 1 ^-T_^^^^+ + ^.C. , (13) 

in the basis of 2-component spinors ^/^+ = {—iw^, ^2)'^, ^~ = {—iw~, h^)^. We define the 
unitary 2 by 2 chargino mixing matrices, U and V, such that: 

- = , (14) 

X""^ diag(»n-+) X+ 

where the (2-component) charginos xt ^re defined such that their absolute masses increase 
with increasing i and such that the mass matrix, m~+, is strictly real: 

Xi 



diag(m^+) = [UM^^V^]^^ = m^6,, . (15) 
2.43 SFERMION MIXING 

At present, we restrict our attention to left-right mixing in the third generation sfermion sector 
only. The convention we use is, for the interaction eigenstates, that /l and fn refer to the 
SU (2) L doublet and singlet superpartners of the fermion / G {t, b, r}, respectively, and, for the 
mass eigenstates, that fi and /2 refer to the lighter and heavier mass eigenstates, respectively. 
With this choice of basis, the spectrum output should contain the elements of the following 
matrix: 



/i 



Fii F12 
F21 F22 



(16) 



whose determinant should be ±1. We here deliberately avoid notation involving mixing angles, 
to prevent misunderstandings which could arise due to the different conventions for these angles 
used in the literature. The mixing matrix elements themselves are unambiguous, apart from the 
overall signs of rows in the matrices, see above. 



2.44 HIGGS MIXING 

The conventions for /i, vi,V2, v, tan P, and m\ were defined above in sections l2!2l and l231 The 
angle a we define by the rotation matrix: 



cos a 
— sin a 



sm a 
cos a 



HI 



(17) 
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where and H2 are the CP-even neutral Higgs scalar interaction eigenstates, and h and 
the corresponding mass eigenstates (including any higher order corrections present in the 
spectrum calculation), with ruho < rriHO by definition. 

2.5 RUNNING COUPLINGS 

In contrast to the effective definitions adopted above for the mixing matrices, we define the 
gauge couplings, the Yukawa couplings, and the soft breaking Lagrangian terms which appear 
in the output as DR running parameters, computed at a user-specifiable scale Q (or grid of 
scales Qi, see below). 

That the DR scheme is adopted for the output of running parameters is simply due to the 
fact that this scheme substantially simplifies many SUSY calculations (and hence all spectrum 
calculators use it). However, it does have drawbacks which for some applications are serious. 
For example, the DR scheme violates mass factorization as used in QCD calculations [28]. For 
consistent calculation beyond tree-level of processes relying on this factorization, e.g. cross 
sections at hadron colliders, the MS scheme is the only reasonable choice. At the present level of 
calculational precision, this is fortunately not an obstacle, since at one loop, a set of parameters 
calculated in either of the two schemes can be consistently translated into the other [29], see 
also [19] for explicit prescriptions. 

Note, however, that different spectrum calculators use different choices for the underlying 
particle content of the effective theory. The programs Softsusy (v. 1.8), SPheno (v. 2.1), and 
Suspect (v. 2.2) use the full MSSM spectrum at all scales, whereas in Isajet (v. 7.69) a more 
involved prescription is followed, with different particles integrated out of the effective theory 
at different scales. Whatever the case, these couplings should not be used 'as is' in calculations 
performed in another renormalization scheme or where a different effective field content is 
assumed. 

Unfortunately, ensuring consistency of the field content assumed in the effective theory 
must still be done on a per program basis, though information on the prescription used by a par- 
ticular spectrum calculator may conveniently be given as comments, when running parameters 
are provided. 

Technically, we treat running parameters in the output in the following manner: since 
programs outside the spectrum calculation will not normally be able to run parameters with 
the full spectrum included, or at least less precisely than the spectrum calculators themselves, 
an option is included to allow the spectrum calculator to write out values for each running 
parameter at a user-defined number of logarithmically spaced scales, i.e. to give output on 
running parameters at a grid of scales, Qi, where the lowest point in the grid will normally be 
niz and the highest point is user-specifiable. A complementary possibility is to let the spectrum 
calculator give output for the running couplings at one or more scales equal to specific sparticle 
masses in the spectrum. 

3. DEFINITIONS OF THE INTERFACES 

The following general structure for the SLHA files is proposed: 

• All quantities with dimensions of energy (mass) are implicitly understood to be in GeV 
(GeV/c^). 

• Particles are identified by their PDG particle codes. See [19] for lists of these, relevant to 
the MSSM. 
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• The first character of every line is reserved for control and comment statements. Data 
lines should have the first character empty. 

• In general, formatted output should be used for write-out, to avoid "messy-looking" files, 
while a free format should be used on read-in, to avoid misalignment etc. leading to 
program crashes. 

• A "#" mark anywhere means that the rest of the line is intended as a comment to be 
ignored by the reading program. 

• All input and output is divided into sections in the form of named "blocks". A "BLOCK 
xxxx" (with the "B" being the first character on the line) marks the beginning of entries 
belonging to the block named "xxxx". E.g. "BLOCK MASS" marks that all following 
lines until the next "BLOCK" (or "DECAY") statement contain mass values, to be read in a 
specific format, intrinsic to the MASS block. The order of blocks is arbitrary, except that 
input blocks should always come before output blocks. 

• Reading programs should skip over blocks that are not recognized, issuing a warning 
rather than crashing. Thereby, stability is increased and private blocks can be constructed, 
for instance BLOCK MYCODE could contain some parameters that only the program My- 
CODE (or a special hack of it) needs, but which are not recognized universally. 

• A line with a blank first character is a data statement, to be interpreted according to 
what data the current block contains. Comments and/or descriptions added after the data 
values, e.g. "... # comment", should always be added, to increase readability of 
the file for human readers. 

Finally, program authors are advised to check that any parameter relations they assume in their 
code (implicit or explicit) are obeyed by the parameters in the files. For instance, tree-level 
relations should not be used with loop-corrected parameters. 

For the technical specifications of the blocks contained in the SUSY Les Houches Accord 
files the full writeup [19] should be consulted. 

4. OUTLOOK 

The present Accord [19] specifies a unique set of conventions together with ASCII file formats 
for model input and spectrum output for most commonly investigated supersymmetric models, 
as well as a decay table file format for use with decay packages. 

With respect to the model parameter input file, mSUGRA, mGMSB, and mAMSB scenar- 
ios can be handled, with some options for non-universality. However, this should not discourage 
users desiring to investigate alternative models; the definitions for the spectrum output file are 
at present capable of handling any CP and R-parity conserving supersymmetric model, with the 
particle spectrum of the MSSM. Specifically, this includes the so-called SPS points [30]. 

Also, these definitions are not intended to be static solutions. Great efforts have gone into 
ensuring that the Accord may accomodate essentially any new model or new twist on an old one 
with minor modifications required and full backwards compatibility. Planned issues for future 
extensions of the Accord are, for instance, to include options for R-parity violation and CP 
violation, and possibly to include definitions for an NMSSM. Topics which are at present only 
implemented in a few codes, if at all, will be taken up as the need arises. Handling RPV and 
CPV should require very minor modifications to the existing structure, while the NMSSM, for 
which there is at present not even general agreement on a unique definition, will require some 
additional work. 
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Part III 

Web Tool For The Comparison Of Susy 
Spectrum Computations 

B. C. Allanach, S. Kraml 
Abstract 

We present and describe an internet resource which allows the user to 
compare different calculations of MSSM spectra. After providing (cur- 
rently mSUGRA) SUSY breaking input parameters, the spectra pre- 
dicted by the publicly available programs ISASUGRA, SOFTSUSY, 
SPHENO and SUSPECT are output by the resource. The variance and 
range of results is also produced. 

1. INTRODUCTION 

Several publicly-available computer programs exist that calculate the MSSM spectrum consis- 
tent with current data on particle masses and gauge couplings, and a theoretical boundary con- 
dition on SUSY breaking. Given the experimental accuracies that are expected for SUSY anal- 
yses at both the LHC and a future e+e" Linear Collider, theoretical uncertainties in spectrum 
computations are important to consider in the total uncertainty of any fit to a SUSY breaking 
pattern. 

As was pointed out in Ref. [31], important sources of such uncertainties are the treatement 
of thresholds in the renormalization group (RG) running, and SUSY loop corrections to the top 
and bottom Yukawa couplings. There has in fact been much progress recently in improving 
the spectrum calculations in commonly used public codes around 'tricky' comers of the SUSY 
parameter space, such as large tan (3 or large mo. However, depending on the specific parameter 
point chosen, the differences in the results of various state-of-the-art codes may still be of the 
same order as or even larger than the expected experimental accuracies. Differences in earlier 
program versions tend to be significantly larger. 

2. ONLINE SPECTRUM COMPARISON 

A pragmatic approach, which was also used in Ref. [31], is to estimate the to-date uncertainty as 
the spread in the results of the most advanced public codes. As mentioned above, this 'compu- 
tational uncertainty' varies over the SUSY parameter space and should therefore be evaluated 
for each particular benchmark point. There also exist several private RG codes, which their 
authors might like to compare to the available public ones in an easy way. Moreover, it can be 
useful to check the results of older program versions against newer ones. 

For these reasons we have set up a web application which allows to compare the results 
of Isajet [11], Softsusy [1], Spheno [5], and Suspect [3] online. The location is 

http : / / cern . ch/ kraml/ comparison/ 

Here the user can input a mSUGRA parameter point ^ and choose the program versions to 
'At the moment of writing, only the mSUGRA model is supported. Other models may be added at a later stage. 
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Table 1: Sparticle masses as obtained by Isajet 7.69, Softsusy 1.8, Spheno2.1.3 and Suspect 2. 1.0. 2 for SPSla. The 
uncertainty A is calculated as A{x) = 0.5[max(.T) — min(x)]. All values are in GeV. 
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compare. On clicking the submit button he then gets a list of sparticle masses from the 
four codes together with the mean, the range and the variance of the results. Note that for the 
Standard Model input the default values of the various codes are used. 

Figure [l] shows a screenshot of the webpage. The application was set up for the Les 
Houches workshop in June 2003. By 31 Oct 2003, it was used by over 30 different users about 
twice a day on average. 

3. RESULTS FOR SPSIA AND SPS2 

In order to give a concrete example, we list in Table [T] some sparticle masses as obtained by 
today's most recent program versions for the SPSla benchmark point [30] (mo = 100 GeV, 
mi/2 = 250 GeV, Aq = -100 GeV, tan/? = 10, /i > 0, nit = 175 GeV). As can be seen, the 
relative differences amout to about 1-2% at SPSla. 

The agreement is less good for neutralino and chargino masses at SPS2 (mo = 1450 GeV, 
mi/2 = 300 GeV, Aq = GeV, tan/5 = 10, /i > 0, rut = 175 GeV), as shown in Table El 
The differences amout to 3-7% due to the notoriously difficult calulation of the /i parameter 
for large mo. Here note that a variation of the input m^ by 1 GeV has a similar effect on the x° 
and masses. The reason is that large cancellations make /i extremely sensitive to the precise 
value of the top Yukawa coupling. TableElshows, however, an order-of-magnitude improvement 
compared to older program versions, where huge discrepancies have been encountered at large 
mo. 

We note that the effect of going from 2 to 3-loop renormalisation group evolution [32] is 
comparable in size to the differencies we find between the latest 2-loop RGB codes. 
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Table 2: Neutralino masses as obtained by Isajet7.69, Softsusy 1.8, Spheno2.1.3 and Suspect2. 1.0.2 for SPS2 
(to-± ~ m^o, rn^± ~ ™xo)- The uncertainty A is calculated as A(a;) = 0.5[max(a;) — min(a;)]. All values are in 
GeV 
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Comparison of spaiticle mass predictions 



http://ki"aml.home.ceni.ch/kraml/comparison/ 



Comparison of sparticle mass predictions 

On this website you can compare tine mass spectra of four public 
SUSY computational tools: Isajet, Softsusy, Spheno and Suspect. 

Theoretical background: JHEP03{2003)016 [ hep-ph/0302102 ] 
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Figure 1 : Screenshot of the online spectrum comparison webpage. 
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Part IV 

Uncertainties in Relic Density Calculations 
in mSUGRA 

B. Allanach, G. Belanger, F. Boudjema, A. Pukhov, W. Porod 
Abstract 

We compare the relic density of neutralino dark matter within the min- 
imal supergravity model (mSUGRA) using four different public codes 
for supersymetric spectra evaluation. 

1. INTRODUCTION 

One of the most stringent constraints on supersymmetric models with R-parity conservation 
arises from the upper limit on the relic density of dark matter. This is particularly true with the 
recent precise measurements of the cosmological parameters realised by WMAR It is therefore 
crucial to quantify the theoretical uncertainties that enter the calculation of the relic density 
of the lightest supersymmetric particle (LSP) and to see how they reflect on the allowed pa- 
rameter space. We do not attempt to answer this question fully here. We will only consider 
one aspect: the uncertainty introduced by the calculation of the weak scale SUSY parameters 
using renormalization group equations (RGE) within the context of the mSUGRA model. As 
a measure of the theoretical uncertainty on the mSUGRA parameters, we use the four public 
state-of-the-art RGE codes: Isajet? . 69 [33], SOFTSUSYl .8.3 [1], SPHEN02 .20 [5] 
and Suspect 2 . 2 [3], link them to micrOMEGAs 1 . 2 [9] and compare estimates for the relic 
density. At this point no attempt is made to estimate the uncertainties that could arise directly 
in the calculation of the relic density itself. 

2. RGE CODES AND RELIC DENSITY CALCULATION 

A detailed study of theoretical uncertainties on the supersymmetric spectra as obtained by RGE 
codes was presented in [31]. It was shown that differences in masses less than a few percent 
are usually found, although some comers of parameter space are still difficult to tackle and can 
display much larger differences. The discrepancies can be traced back to the level of approxi- 
mation used in the weak-scale boundary conditions. The large tan /? region and the focus point 
region (large Mq) are still subject to large theoretical errors. Both of these regions are pre- 
cisely where one can find cosmologically interesting values for the relic density, ^Ih"^ < .128. 
In the focus point region, the LSP is mainly a Higgsino and annihilates efficiently into gauge 
bosons. At large tan P, even rather heavy neutralinos can annihilate into bb pairs via s-channel 
exchange of a heavy Higgs. The coannihilation region where the Next-to-Lightest supersym- 
metric particle (NLSP) is nearly degenerate in mass with the LSP, is another cosmologically 
relevant region. Although it is a priori not difficult to handle by the RGE codes, the value of 
the relic density depends sensitively on the mass difference between the NLSP and the LSP and 
even shifts of C(l) GeV can cause large shifts in the relic density. The other cosmologically 
viable mSUGRA region, the bulk region, shows a much smaller induced sensitivity upon the 
MSSM mass spectrum. 
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Figure 1: a) flh'^ for Mq = 150 GeV, tan/3 = 10, Ao = 0, a* > for SOFTSUSYl . 8 . 3 (full), 
Isajet7.69 (dashed), Suspect2.2 (dash-dotted), and SPHENO2.20 (dotted). At large A/1/2, 
Isa jet7 . 69 and Suspect2 . 2 give nearly identical results, b) ilh^ vs m^^ — m^o for the same set of pa- 
rameters as a). 



The link between micrOMEGAsl . 2 and the RGE codes is done within the spirit of the 
SUSY Les Houches Accord [19] : common input values are chosen and pole masses, mixing 
matrices, the /i parameter and the trilinear couplings are calculated by the RGE codes. All 
parameters are read by micrOMEGAsl .2 1.2. The annihilation cross-sections are then evalu- 
ated at tree-level. Important radiative corrections to the Higgs widths and in particular the Ami, 
correction are taken into account. 

3. RESULTS 

For the numerical results as default values we have fixed nit = 175 GeV, as{Mz)^^ = .1172 
and mb{mb)^'^^ = 4.16 GeV. This corresponds to mb{Mz)^^ = 2.83 GeV. We concentrate 
on the three regions where the relic density is within the WMAP range and where potentially 
large discrepancies can be observed: the focus point region, the large tan/? region and the 
coannihilation region. 

3.1 Coannihilation 

Mo = 150 GeV, Aq = 0, tan/3 = 10, /i > 

The small M1/2 region corresponds to the so-called bulk region where the bino-LSP annihilates 
into lepton pairs via s-channel Z or Higgs exchange or t-channel slepton exchange. Here one 
finds very good agreement between the values of Qh^us'mg the different RGE codes (see Fig.[lt) 
since the predicted values for slepton and neutralino masses are in good agreement (within a 
few GeV). The exact position of the Z pole (corresponding to the big dip in is slightly 
shifted for SPHEN02 . 2 but the range of values of M1/2 for which Qh'^ < .128 are basically 
identical. Note that the Z pole region is ruled out by the LEP constraints on neutralinos within 
the context of mSUGRA models. 
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Figure 2: a) ilh^ as a function of Mq for AI1/2 ~ 300 GeV, tan/? = 10, = and ^ > and mt = 
175. Same labels as in Fig. 1. b) Dependence of the reUc density on rrit for SOFTSUSYl .8.3 (full) and 
SPHEN02 .20 (dash) . 



As one moves up in M1/2, one reaches the so-called coannihilation region where the f is 
the NLSP and is nearly degenerate with the neutralino, as in Fig.[lj5. Coannihilation with the f, 
and to a lesser extent the selectron and smuon, brings the relic density in the desired range. For 
a given value of M1/2, differences between the codes can reach a factor 2, the largest differences 
are found between SPHEN02 . 2 and SOFTSUSYl .8.3. However very good agreement is 
found between all codes when the relic density is plotted as a function of the mass difference 
between the LSP and the NLSP (here the f). All codes obtain values of fi/i^compatible with 
WMAP for mass differences — m^o ^ 4 GeV (at the extreme left of Fig.Ilt), even though 
the corresponding value of the neutralino mass can differ. The value of M1/2 for which the relic 
density becomes compatible with WMAP varies from 670 GeV (SPHEN02 . 2 0) to 790 GeV 
(SOFTSUSYl . 8 . 3), a 12% difference on M1/2. 

3.2 Focus point 

Ml/2 = 300 GeV, Aq = 0, tan /5 = 10, /i > 

In addition to the small Mq (bulk/coannihilation region) where annihilation into leptons is im- 
portant, the cosmologically relevant region is found at values of Mq well above ITeV As one 
approaches the region where electroweak symmetry breaking is forbidden, the /i parameter ap- 
proaches zero. This means that the LSP is mainly Higgsino. This LSP can then annihilate very 
efficiently into gauge bosons (WW/ZZ) and to a lesser extent into Zh. The parameter fi is how- 
ever very sensitive [34] to the top Yukawa coupling, ht (which is also reflected in a sensitivity to 
the value of the top quark mass) and huge differences between codes were observed [31]. The 
impact on the relic density and on the exclusion region is likewise very significant. 

As can be seen in Fig. |2l all codes agree very well for Mq < ITeV but as one gets to 
large values of Mq, more than one order of magnitude differences in can be found. For 
rrit = 175 GeV, only Isajet finds a large drop in the /i parameter as one moves to Mq ^ 3000 
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Figure 3: a) ^lh^ as a function of A/1/2 for mt ~ 175, tan/? = 52, Aq ~ and /i > 0. Same labels as in Fig. 1. 
b) Dependence of the relic density on mbinib) for SOFTSUSYl .8.3 (full) and SPHEN02 . 2 (dash). 



GeV, this is when fi/i^ drops below the upper limit from WMAP. The other codes do not find 
this drop in fi and do not obtain a cosmologically interesting region for Mq < 4000 GeV. 
These large differences between codes however are just a reflection of the sensitivity to the top 
Yukawa, ht{MsusY) which is proportional to mj. We show in Fig. 2b, the variation of VLh^ with 
mt using SOFTSUSYl . 8 . 3 and SPHEN02 . 20 for Mq = 3000 GeV. The value ^Ih^ = .128 
found in Isa jet7 . 69 formt = 175 GeV can be reproduced in SOFTSUSYl .8.3 (SPHENO) 
by changing the input to = 172.2(172.5) GeV. 

3.3 Large tan /3 

mi/2 = 1500 GeV, Aq = 0,tan/5 = 52 ^ > 

At large tan P the new feature is the annihilation of neutralinos into bb via heavy Higgs ex- 
change. With the current version of the RGB codes, this is observed only for very large values 
of tan p. The crucial parameter here is Af^/2m^o which must be close to unity to provide suf- 
ficient annihilation of neutralinos. Large differences in the value of Ma between the different 
RGB codes occur because of the sensitivity of the RGB to the bottom Yukawa as well as from 
taking into account higher loop effects. 

As Fig. 3a shows, all 4 programs predict a large drop in the relic density when the neu- 
tralino mass gets close to Ma/2 although this drop occurs at much lower values of Afi/2 for 
SPHENO, Ml/2 ~ 1250 GeV than for Isa jet7 . 69 , M1/2 ^ 1750 GeV. However, here again 
the results are very sensitive to the input parameters, in this case the value of the b-quark mass. 
For Ml/2 = 1300 GeV, we find an order of magnitude shift in Qh"^ for mb{mb) = 4 — 4.4 GeV 
with the program SOFTSUSYl. 8. 3. By a slight shift of the b-quark mass we can find perfect 
agreement between SPHENO 2 . 2 and SOFTSUSYl . 8 . 3, as shown in Fig. 3b. 



23 



4. CONCLUSION 

While the predictions for the relic density of neutralinos are rather stable in most of the mSUGRA 
space, it is in the most physically interesting region that large discrepancies can be observed, 
in particular the focus point/large tan P and coannihilation regions. It is however reassuring to 
find that with the newer versions of the codes, the discrepancies in the sparticle spectra tend to 
be reduced. More details on the theoretical uncertainties in the evaluation of the relic density 
arising from the standard model parameters, nib, nit, used as input in a RGE code can be 
found in [35]. 

ACKNOWLEDGEMENTS 

This work was partially supported by CNRS/PiCS-397, Calculs automatiques de diagrammes 
de Feynman. We thank Jean-Loic Kneur for providing an improved version of Suspect 2 . 2 . 



24 



Part V 

SFITTER: A Tool To Determine 
Supersymmetric Parameters 

R. Lafaye, T. Plehn, D. Zerwas 
Abstract 

SFITTER is a new tool to determine supersymmetric model parameters 
from collider measurements. It allows to perform a grid search for the 
minimal and/or a fit of a given model. Currently, the model param- 
eters in the general MSSM or in a gravity mediated SUSY breaking 
model can be tested using a given set of mass, branching ratio and cross 
section measurements. 

1. Introduction 

The most important task for the LHC as well as for any future Linear Collider is to study 
in detail the mechanism which leads to electroweak symmetry breaking. While the Standard 
Model describes all available high energy physics experiments, it still has to be regarded as an 
effective theory, valid at the weak scale. New physics are expected to appear at the TeV energy 
scale. The minimal supersymmetric extension of the Standard Model (MSSM) can provide a 
description of physics up to the unification scale. 

If supersymmetry or any other high-scale extension of the Standard Model is discov- 
ered, it will be crucial to determine its fundamental high-scale parameters from weak-scale 
measurements [36, 37]. The LHC and a future Linear Collider will provide a wealth of mea- 
surements [38], which due to their complexity require a proper treatment to unravel the corre- 
sponding high-scale physics. Even in the general weak-scale MSSM without any unification 
or SUSY breaking assumptions the measurements of masses and couplings are not likely to 
be independend measurements; moreover, linking supersymmetric particle masses to weak- 
scale SUSY parameters involves non-trivial mixing to mass eigenstates in essentially every 
sector of the theory. On top of that, for example in gravity mediated SUSY breaking scenarios 
(mSUGRA/cMSSM) a given weak-scale SUSY parameter will always be sensitive to several 
high-scale parameters which contribute through renormalization group running. Therefore, a 
fit of the model parameters using all experimental information available will lead to the best 
sensitivity and make the most efficient use of the information available. 

If the starting point of the fit is not known and many parameters are involved, the allowed 
parameter space might not be sampled completely in the fit approach. To avoid boundaries 
imposed by non-physical parameter points, which can confine the fit to a 'wrong' parameter 
region, combining the fit with an initial evaluation of a multi-dimensional grid is the optimal 
approach. In the general MSSM the weak-scale parameters can vastly outnumber the collider 
measurements, so that a complete parameter fit is not possible and one has to limit oneselve to a 
subset of parameters. In SFITTER both grid and fit are realised and can be combined, including 
a general correlation matrix and the option to exclude parameters of the model from the fit/grid 
by fixing them to a value. 
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2. SFITTER — Program Structure 

Currently, SFITTER uses the predictions for the supersymmetric masses provided by SUS- 
PECT [3], but the conventions of the SUSY Les Houches accord [19] could be helpful, if pro- 
vided as a common block/C- structure, to ease interfacing other programs. The branching ratios 
and e+e~ production cross sections are provided by MSMlib [39], which has been used exten- 
sively at LEP and cross checked with Ref. [40]. The next-to-leading order hadron collider cross 
sections are computed using PROSPINO [41^3]. The fitting program uses the MINUIT pack- 
age [44]. The determination of includes a general correlation matrix between measurements. 
For unphysical points in supersymmetric parameter space, is set to 10^°. 

2.1 Initialization and Steering 

The program SFITTER is driven by two files: the first one sets up the measurements and the 
corresponding errors. For each measurement one specifies if it is to be used in the grid (G) or 
in the MINUIT fit (M) or in both. 

//set all errors to 0.5% of their central value 
D AT A_ERR = 0.005 

//randomize the measurements around their nominal value 
RANDOMIZE = 1 

//Higgs mass and error to be used in the Fit only 
m_h = 112.6 +/- 0.1 [-/M] 

/ /Neutralinol mass to be used in Grid and Fit 
m_chiO_l = 180.2 +/- 5.1 [G/M] 

//Correlation between two chargino mass measurements 
C0RR(m_chi+_l,m_chi+_2) =0.03 

The second file initializes everything related to the weak-scale or high-scale MSSM model 
parameters. First the model (mSUGRA, pMSSM etc) is specified, then the starting values of 
all MSSM parameters, boundaries, stepsize and the number of points in the grid are specified. 
Moreover, the user defines if a certain MSSM model parameter is included in the grid and in 
the fit: 

MODEL=MSUGRA // use MSUGRA 
// use the GRID (or not) 
GRID=1 

// MO used in grid and fit, grid of 10+1 steps between and 1000. 
M0=500. [M/G] STEP=200. LOW=0 . HIGH=1000. GRID=10 
// AO used only in fit 

A0=0. [M/-] STEP=200. LOW=-1000. HIGH=1000. 

2.2 mSUGRA/cMSSM Parameter Determination 

Assuming that SUSY breaking is mediated by gravitational interactions (mSUGRA/cMSSM) 
we fit four universal high-scale parameters to a toy set of collider measurements: the univer- 
sal scalar and gaugino masses, mo, mi/2, the trilinear coupling Aq and the ratio of the Higgs 
vacuum expectation values, tan p. The sign of the Higgsino mass parameter is a discrete 
parameter and therefore fixed. The assumed data set is the set of all supersymmetric particle 
masses for the SUSY parameter point SPS la [30,45], as computed by SUSPECT. The errors on 
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Table 1: Left: summary of mSUGRA fit in SPSla: true values, starting values, fit values. As in SPSla we fix 
IJ. > 0. All mass values are given in GeV. Right: the (symmetric) correlation matrix of all SUSY parameters in the 
mSUGRA fit. 



the toy mass measurements are uniformly set to 0.5%. The starting points for the mSUGRA pa- 
rameters are fixed to the mean of the lower and upper limit in the fit, i.e. they are not necessarily 
even close to the true SPSla values. The result of the fit is shown in Tab. [T] With SFITTER 
the true parameter values were reconstructed well within the quoted errors, in spite of starting 
values relatively far away from the true ones. The measurement of mo and mi/2 is very precise, 
while the sensitivity of the masses on tan P and Aq is significantly weaker. 

The correlations between the different high-scale SUSY parameters are also given in 
Tab.[T] One can understand the correlation matrix step by step [46]: first, the universal gaugino 
mass mi/2 can be extracted very precisely from the physical gaugino masses. The determina- 
tion of the universal scalar mass mo is dominated by the weak-scale scalar particle spectrum, 
but in particular the squark masses are also strongly dependent on the universal gaugino mass, 
because of mixing effects in the renormalization group running. Hence, a strong correlation 
between the mo and mi/2 occurs. The universal trilinear coupling Aq can be measured through 
the third generation weak-scale mass parameters Ah^t,T- However, the Ah^t,T which appear for 
example in the off-diagonal elements of the scalar mass matrices, also depend on mo and mi/2, 
so that Aq is strongly correlated with mo and mi/2. 

In the SPSla scenario, the pseudoscalar Higgs is heavy and the Higgs masses do not 
show a strong dependence on tan p. Because of the large mass difference between gauginos 
and Higgsinos they essentially decouple, and the neutralino/chargino sector will not yield a 
good determination of tan/3. The stop mixing is governed by At, and not hy fi/ tan/5, while 
the sbottom mixing is small altogether. Only the stau mixing is large and driven by /itan/5 
in the off-diagonal element of the stau mass matrix. The stau mass parameters are dominated 
by mo, in particular the smaller right handed stau mass. Therefore, one expects tan/3 to be 
strongly correlated with mo and less with mi/2. The result from SFITTER as shown in Tab.[T]is 
in agreement with this prediction. Thus, the results obtained with SFITTER can be understood 
from the particular features of the SPSla spectrum. 



2.3 MSSM Parameter Determination 

In total 24 parameters describe the unconstrained weak-scale MSSM. They are listed in Tab.|2j 
tan (3 just like in mSUGRA, plus three soft SUSY breaking gaugino masses Mj, the Higgsino 
mass parameter fi, the pseudoscalar Higgs mass uia, the soft SUSY breaking masses for the 
right sfermions, Mj^ , the corresponding masses for the left doublet sfermions, Mj^ and finally 
the trilinear couplings of the third generation sfermions At fibr- 
in any MSSM spectrum, in first approximation, the parameters Mi, M2, /i and tan/3 
determine the neutralino and chargino masses and couplings. We exploit this feature to illustrate 
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Table 2: Result for the general MSSM parameter determination in SPS 1 a. Shown are the nominal parameter values, 
the result after the grid and the final result. The deviation in the squark sector of 1% is an artefact of differences 
between MSSM and mSUGRA part of the renomalization group code [3]. All masses are given in GeV. 

the option to use a grid before the start the complete MINUIT fit. For testing purposes, the error 
on all mass measurements is again set 0.5%. The starting values of the parameters are set to their 
nominal values, this study is thus less general than the one of mSUGRA. Then we minimize 
on a grid. For this grid minimization the six chargino and neutralino masses are used as 
measurements to determine the four SUSY parameters Mi, M2, /i and tan /5 only. The step size 
of the grid is 10 for tan P and 100 GeV for the three mass parameters. After the minimization, 
these four parameters obtained from the minimum on the grid are fixed and all remaining 
parameters are fitted. Only in a final run all SUSY parameters are released and fitted, to give 
the final results quoted in Tab. |2| 

In Tab.|2lthe intermediate (after the grid evaluation) results, the final results and the nom- 
inal values are shown. The final fit values indeed converges to the correct central values within 
its error. The central values of the fit are in good agreement with generated values, except for 
the trilinear coupling A^. As already mentioned in the discussion of the mSUGRA fit, the mix- 
ing between the two sbottom mass states is very small, so the assumed precision of the 0.5% 
is insufficient to determine the parameter from the mass measurements alone. As At enters in 
the calculation of the lightest Higgs, additional sensitivity for this parameter comes from the 
mass measurement of the lightest Higgs boson. The use of branching ratios and cross section 
measurements should significantly increase the precision in future studies, especially for A^ 
and Ah. 

3. Conclusions 

SFITTER is a new program to determine suspersymmetric parameters from experimental mea- 
surements. The parameters can be extracted either using a fit, a multi-dimensional grid minimi- 
sation, or a combination of the two. Correlations between measurements can be specified and 
are taken into account in the calculation of the x^- SUSPECT, MSMlib and PROSPINO are 
used to calculate the predictions for the masses, branching ratios and production cross sections. 
A more realistic set of the measurements for example assuming the SPSla mass spectrum for 
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the LHC and and a future Linear Collider will be studied as a next step. The impact of corre- 
lations between measurements on the estimated errors of MSSM parameters will be studied in 
detail. In the future public version of the program we will include different generators for the 
calculation of masses and branching ratios. 
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Part VI 

SDECAY: a Code for the Decays of the 
Supersymmetric Particles 

A. Djouadi, Y. Mambrini and M. Muhlleitner 
Abstract 

We present the Fortran code SDECAY, a program which calculates the 
decay widths and branching ratios of all supersymmetric particles in 
the Minimal Supersymmetric Standard Model, including higher order 
effects. The usual two-body decays of sfermions and gauginos as well 
as the three-body decay modes of charginos, neutralinos and gluinos are 
included. Furthermore, the three-body and even the four-body decays 
of top squarks are calculated. The important loop-induced decays, the 
QCD corrections to the two-body widths involving strongly interacting 
particles and the dominant electroweak effects to all processes are eval- 
uated as well. 

1. Introduction 

The search for new particles predicted by supersymmetric (SUSY) theories is a major goal 
of present and future colliders. In the Minimal Supersymmetric Standard Model (MSSM) [47] 
there are still over 20 free parameters even in a phenomenologically viable model. It is therefore 
a very complicated task to deal with all the properties of the SUSY particles once they are 
found. Since their properties will be determined with an accuracy of a few per cent at the LHC 
and a precision at the per cent level or below at future e+e^ linear colliders, the mass spectra, 
the various couplings, the decay branching ratios and the production cross sections have to be 
calculated with a rather high precision, also including higher order effects. The Fortran code 
SDECAY^ [8] which is presented here calculates the decays of SUSY particles in the MSSM, 
including the most important higher order effects. The Renormalization Group Equation (RGE) 
program SuSpect [3] is used for the calculation of the mass spectrum and the soft SUSY- 
breaking parameters. [Of course, SDECAY can be easily linked to any other RGE code.] Due to 
the limited space we refer for details of the notation, the description of the algorithm that is used 
in the code and the various higher order effects that have been included to the user's manual of 
SuSpect. The program SDECAY then evaluates the various couplings of the SUSY particles 
and MSSM Higgs bosons and calculates the decay widths and the branching ratios of all the 
two-body decay modes, including the QCD corrections to the processes involving coloured 
particles and the dominant electroweak effects to all processes. The loop-induced two-body 
decay channels as well as the possibly important higher order decays are included, such as the 
three-body decays of charginos, neutralinos, gluinos and top squarks and the four-body decays 
of the lighter top squark. In addition, the top quark SUSY decay widths and branching ratios 
are implemented. The program will be presented in the following. 

-The code can be obtained at the url: http://people.web.psi.ch/muehlleitner/SDECAY 
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2. The decays of the supersymmetric particles 

2.1 The tree level two-body decays 

The Fortran code SDECAY includes the two-body decays of sfermions into a fermion and a 
gaugino, as well as into a lighter sfermion of the same isodoublet and a gauge boson V = W, Z 
or a Higgs boson ^ = h,H, A, 

h - x^P (1) 

h ^ vjp (2) 

/. ^ ^fP (3) 
For squarks heavier than the gluino the decay into a gluino-quark final state is also possible 

Qi qg (4) 

The heavier neutralino and chargino decays into the lighter chargino and neutralino states and 
gauge or Higgs bosons as well as the decays into fermion- sfermion pairs have been implemented 

Xi XjV (5) 

Xi Xj"^ (6) 

X^ ^ ffP (7) 
For the gluinos the only relevant decay into a squark-quark pair is calculated 

g qqi (8) 

In the case of a GMSB model the decays of the next-to-lightest SUSY particle (NLSP), which 
can be either the lightest neutralino Xi or the lightest sfermion, in general the fi, into a Gravitino 
G and a photon, Z or neutral Higgs boson (for x?) and a r (for fi) are implemented 

X? ^ jG,ZG,^>G (9) 
fi ^ tG (10) 

The masses entering the phase space in the calculation of the widths are the pole masses, but 
when they enter the various couplings they are - for the third-generation fermions - the running 
DR masses at the scale of Electroweak Symmetry Breaking (EWSB). This is also the case for 
all soft SUSY-breaking parameters and the third generation sfermion mixing angles involved 
in the couplings. In addition, we have left the option for the QCD coupling constant and the 
bottom, top Yukawa couplings to be evaluated at the scale of the decaying superparticle or any 
other scale. In this case, only the standard QCD corrections are included in the running [48]. 

2.2 The QCD corrected two-body decays 

The one-loop QCD corrections to the following two-body decays involving (s)quarks and gluinos 
have been implemented using the formulae of Refs. [49,50], [51-54] and [55,56], respectively, 

qi xjq^'^ (11) 

qi ^ (12) 
Qi qg and g qiq (13) 



All the corrections have been included in the DR scheme. The bulk of the electroweak radiative 
corrections due to the running of the gauge and third-generation fermion Yukawa couplings has 
been taken into account by evaluating these parameters at the EWSB scale. 
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2.3 Loop-induced decays 

In case the two-body decays of the next-to-lightest neutralino are kinematically not allowed the 
loop-induced decay into the lightest supersymmetric particle (LSP) Xi ^ photon is calcu- 
lated [57-60] 

X2 ^ xh (14) 

For completeness, the loop-induced decay of a gluino into a gluon and the LSP has also been 
considered [61-63] 

~9 ^ 9Xi (15) 

If the tree-level stop two-body decays are kinematically closed the loop-induced decay into a 
charm and Xi [64] is calculated 

^ cxI (16) 



2.4 Multibody decay modes 

If the two-body decays of the gauginos Eqs. (jSHll) are kinematically forbidden the three-body 
decays into a lighter gaugino and a fermion pair and a gluino and two quarks are calculated 

X^ ^ XjfP (IV) 
X^ ^ M^'^ (18) 

Analogously, the gluino three-body decays into a gaugino and two quarks are considered when 
the gluino two-body modes are closed 

9 xm^"' (19) 

For the calculation of the processes Eqs. (I17H19I) we have used the formulae given in [65-69]. 
Furthermore, the possibly important gluino decay into stop, bottom and a W boson as well as 
the decay into stop, bottom and a charged Higgs boson have been implemented [70,71] 

g iibW- (20) 
g ^ hbH- (21) 

In case the stop two-body decays are not accessible, there are several three-body decay modes 
[72-77] that can dominate over the loop-induced decay Eq. ( fT^ in rather large areas of the 
MSSM: the decays into a bottom, lightest neutralino and a or charged Higgs boson, the 
decay modes into bottom, lepton and slepton, the decays into the lightest sbottom and a fermion 
pair as well as for the heavy stop the possibility of decaying into the lighter stop and a fermion 
pair 



U ^ bW+Xi , bH+Xi (22) 

ti bl~^i>i and/or bl'^i'i (23) 

ii ^ hff (24) 

h ^ hff (25) 
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SDECAY evaluates the three-body decays if the two-body decays are closed, taking into account 
all possible contributions of virtual particles, the radiatively corrected Yukawa couplings of 
third-generation fermions, the mixing pattern for their sfermion partners and the masses of the 
sparticles and gauge/Higgs bosons involved in the processes. Even the masses of the final state 
fermions have been included. The total decay widths of the exchanged particles have not been 
included in the propagators of the virtual particles. 

If the stop three-body decay channels are kinematically forbidden the ii four-body decay 
mode into a bottom, the LSP and two massless fermions can become competitive with the loop 
induced decay into a charm and a neutralino, cf. Eq. (fT6b . so that this channel [78] has also been 
included in the program, 

ii ^ bxlff (26) 

2.5 Top quark decays 

For the top quark the following decays in the MSSM are calculated by SDECAY 

t bW+ (27) 
t hH+ and tix\ (28) 

3. How to use SDECAY 

Apart from the files of the program SuSpect, i.e. suspect 2 . in, suspect 2 . f , subh_hdec . f , 
f eynhiggs . f , hmsusy . f , the program SDECAY consists of three files: 

1) The input file sdecay . in where one can choose the accuracy of the algorithm and the 
various options whether QCD corrections and multibody or loop decays are included or not, 
which scales and how many loops are used for the running couplings and if top and GMSB 
decays are calculated or not. 

2) The main routine sdecay . f where the couplings of the SUSY and Higgs particles are 
evaluated and the decay branching ratios and total widths are calculated. 

3) The output file sdecay . out which gives the results for the branching ratios and total 
widths, as well as the masses of the SUSY and Higgs particles, the mixing matrices and the 
gauge and third-generation Yukawa couplings at the EWSB or a chosen scale. The output is 
given in two possible formats, either in a simple and transparent form or according to the SUSY 
Les Houches Accord [19] which uses the PDG notation for the particles. 

All these files together with a makefile to compile the files can be found on the web page dedi- 
cated to SDECAY at the address: 

http : / / people . web . psi . ch/muehlleitner/ SDECAY 

4. Conclusions 

We have presented the Fortran code SDECAY, which calculates the decay widths and branch- 
ing ratios of all the two-body decays of the SUSY particles in the framework of the MSSM, 
including the QCD corrections to the decays involving strongly interacting particles, the three- 
body decays of the gauginos, gluinos and stops, as well as the four-body decays of the lightest 
top squark. Furthermore, the loop-induced decays of the gluino, the lightest neutralino and the 
lightest top squark, the decays of the next-to-lightest SUSY particle in GMSB models and the 
standard and SUSY decay modes of the top quark have been implemented. The dominant elec- 
troweak corrections due to the running of the gauge and fermion Yukawa couplings have been 
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incorporated. The program which uses the RGE code SuSpect can be easily linked to any 
other spectrum calculator. It is user-friendly, flexible for the choice of options and approxima- 
tions and quite fast. The program is under rapid development and will be updated regularly. 
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Part VII 

Measuring The Mass Of The Lightest 
Chargino At The CERN LHC 

M.M. Nojiri, G. Polesello and D.R. Tovey 
Abstract 

Results are presented of a feasibility study of techniques for measuring 
the mass of the lightest chargino at the CERN LHC. These results sug- 
gest that for one particular mSUGRA model a statistically significant 
chargino signal can be identified and the chargino mass reconstructed 
with a precision ~ 1 1% for ~ 100 fb~^ of data. 

1. INTRODUCTION 

Much work has been carried out recently on measurement of the masses of SUSY particles at 
the LHC [79-84]. These measurements can often be considered to be 'model-independent' in 
the sense that they require only that a particular SUSY decay chain exists with an observable 
branching ratio. A good starting point is often provided by the observation of an opposite-sign 
same-flavour (OS-SF) dilepton invariant mass spectrum end-point whose position measures a 
combination of the masses of the X2, the Xi and possibly also the Observation of end-points 
and thresholds in invariant mass combinations of some or all of these leptons with additional 
jets then provides additional mass constraints sufficient to allow the individual sparticle masses 
to be reconstructed unambiguously. A question remains however regarding how the mass of a 
SUSY particle can be measured if it does not participate in a decay chain producing an OS-SF 
dilepton signature. This problem has been addressed for some sparticles (e.g. for the qji [85]) 
however significant exceptions remain. Notable among these is the case of the lightest chargino 
xf, which does not usually participate in decay chains producing OS-SF dileptons due to its 
similarity in mass to the X2- 

In this paper we attempt to measure the mass of the xf by identifying the usual OS-SF 
dilepton invariant mass end-point arising from the decay via X2 of the other initially produced 
SUSY particle (i.e. not the one which decays to produce the xf)- We then solve the mass 
constraints for that decay chain to reconstruct the momentum of the Xi appearing at the end of 
the chain, and use this to constrain the momentum (via E^^^^) of the Xi appearing at the end of 
the decay chain involving the xf - We finally use mass constraints provided by additional jets 
generated by this chain to solve for the xf mass. The technique requires that both the decay 
chain 

Ql xlo. Ink XiHq 

and the decay chain 

Ql xfQ qW^Xi C[q'c["Xi 

are open with significant branching ratios, and that the masses of the xi, X2, Ir and ql are 
known. No other model-dependent assumptions are required however. 
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2. SUSY MODEL AND EVENT GENERATION 

The SUSY model point chosen was that used recently by ATLAS for full simulation studies 
of SUSY mass reconstruction [86]. This is a minimal Supergravity (mSUGRA) model with 
parameters nio = 100 GeV, mi/2 = 300 GeV, Aq = -300 GeV, tan(/3) = 6 and /i > 0. The 
mass of the lightest chargino is 218 GeV, while those of the ql, the the X2 ^i^d the x? are 
~ 630 GeV, 155 GeV, 218 GeV and 118 GeV respectively. One of the characteristics of this 
model is that the branching ratio of xf — > W^Xi is relatively large (~ 28 %). Chargino mass 
reconstruction involving the decay xf fiUr (BR ~ 68 %) is likely to be very difficult due 
to the additional degress of freedom provided by the missing neutrino. Consequently the 
decay mode must be used. 

The electroweak SUSY parameters were calculated using the ISASUGRA 7.51 RGB code 
[1 1]. SUSY events equivalent to an integrated luminosity of 100 fb^^ were then generated using 
Herwig 6.4 [12,87] interfaced to the ATLAS fast detector simulation ATLFAST 2.21 [88]. With 
the standard SUSY selection cuts described below Standard Model backgrounds are expected 
to be negligible. An event pre-selection requiring at least two ATLFAST-identified isolated 
leptons was applied in order to reduce the total volume of data. 

3. CHARGINO MASS RECONSTRUCTION 

Events were required to satisfy 'standard' SUSY selection criteria requiring a high multiplicity 
of high pt jets, large E^'^^'^ and multiple leptons: 

• at least 4 jets (default ATLFAST definition [88]) with pr > 10 GeV, two of which must 
havepT > 100 GeV, 

• (EtiPTOe*) + ^r^^) >400GeV 

• > max(^100GeV, 0.2 (EtiPT{jet) + ^t "'))' 

• exactly 2 opposite sign same flavour isolated electrons or muons with px > 10 GeV, 

• no b-jets or r-jets. 

Events were further required to contain dileptons with an invariant mass less than the 
expected /^/^ end-point position (100.2 GeV) and at least one dilepton -1- hard jet combination 
(one for each combination of the dilepton pair with each of the two hardest jets) with an invariant 
mass less than the expected /^/^g end-point position (501.0 GeV). The smaller dilepton -1- hard 
jet combination then defined which jet (assumed to be from the decay ql X2Q) would be used 
together with the dileptons to reconstruct the X2 production and decay chain. 

The momentum of the x? at the end of the X2 decay chain was calculated by solving ana- 
lytically the kinematic equations relating the momenta of the decay products (including the x?) 
to the masses of the SUSY particles, which were assumed to be known from conventional end- 
point measurements [79-84]. This process is described in more detail in Ref. [89] and results in 
two solutions for the x? momentum for each of the two possible mappings of the reconstructed 
leptons to the sparticle decay products. In the present analysis just one such mapping was as- 
sumed with no attempt being made to select the correct assignment. Two possible solutions for 
the Xi momentum were therefore obtained for each event. 

The nest step in the reconstruction was to find the jet pair resulting from a hadronic 
decay following production via xf W^^Xi- The potentially large combinatorial background 
was reduced by rejecting jet combinations involving either of the two hardest jets (since these 
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Figure 1: Reconstructed dijet invariant mass distributions for all events (data points) and events not containing 
the decay chain xf — *■ M^^Xi ^ q'<i"Xi selected using Monte Carlo truth. The signal band is labelled T' in the 
figure, while the two sideband are labelled 'IF and TIF respectively. 



were assumed to arise from ql decay) and by requiring that the harder( smaller) of the two jets 
possessed pt greater than 40(20) GeV (i.e. selecting asymmetric jet pairs consistent with a 
significant boost in the lab frame). A further cut was applied on the invariant mass of the com- 
bination of the jet pair with the hard jet giving the larger dilepton + jet mass (assumed therefore 
to be the jet from the ql xfl decay preocess). This invariant mass was conservatively 
required to be less than that of the ql- 

For each event any jet pairs satisfying the above criteria and possessing \mjj — mw\ < 
15 GeV (Fig.[T]), were considered to form W candidates. For each event the candidate with rrijj 
nearest mw was then selected and used together with the momentum of the hard jet identified 
previously and the two assumed x and y components of the Xi momentum (calculated from 
the two solutions for the momentum of the Xi from the X2 decay and E^^'^^) to calculate the 
chargino mass. Each of the two solutions for the momentum gives two possible solutions for 
m-±, the smaller of which is usually physical. Consequently two possible values for m-± were 
obtained from each event (plotted in Fig.|2|). 

Following this procedure significant backgrounds remain from combinatorics in SUSY 
signal events (due to their high average multiplicity), and from SUSY background events (i.e. 
events in which the decay process ql xfq qW'^Xi ~^ '1Q''1"Xi is not present). These 
backgrounds (or at least those not involving a real decay) were removed statistically using 
a sideband subtraction technique similar to that described in Ref . [90] . All jet pairs satisfying all 
the above selection criteria except the \mjj—miv \ requirement were recorded if they satisfied the 
alternative requirement that 15 GeV < \mjj — mw\ < 45 GeV. This requirement then defined 
two side-bands located on either side of the main signal band {\mjj — mw\ < 15 GeV) of 
equal width 30 GeV. The momentum of each jet pair was then rescaled such that the difference 
between its rescaled mass and mw was the same as the difference between its original mass and 
the centre of its sideband (50 or 1 10 GeV respectively). Each jet pair was then given a weight of 
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Figure 2: Reconstructed xf mass distributions showing signal distributions with \mjj — mw\ < 15 GeV (data 
points) and sideband distributions with 15 GeV < \mjj — mw\ < 45 GeV (histograms). The left hand figure was 
obtained by selecting events containing the decay chain W^Xi ~^ q'q"Xi using Monte Carlo truth. The 

central figure was obtained by selecting background events not containing this decay chain. The right hand figure 
was obtained by using all data. 



1 .3 (lower sideband) or 1 .0 (upper sideband) to account for the variation of the background rrijj 
distribution with nijj (Fig.[T]). Values for the chargino mass were then calculated for each jet pair 
and used to create a sideband mass distribution (Fig. El). Finally the sideband mass distribution 
was subtracted from the signal mass distribution with a relative normalisation factor of 0.7 to 
account for the differing efficiencies for selecting sideband events and background events in the 
signal region. 

4. RESULTS 

The sideband subtracted chargino mass distributions obtained from this process are shown in 
Fig. 121 both with and without a selection requirement for xf — > VT^Xi ~^ q'^''^.! obtained from 
Monte Carlo truth. In both cases no events are observed at masses below the kinematic limit 
of 198 GeV (= mw + fn^a) due to the origin of the mass values as solutions to the kinematic 
mass relations. In the case where Monte Carlo truth was used as input a clear peak is seen in 
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Figure 3: Reconstructed xf mass distribution for xf W Xi ~^ q'q"Xi signal events (histogram) and all 
events (points with errors). 



the 200 GeV - 250 GeV bin, corresponding well to the actual mass of 218 GeV. At higher mass 
values the sideband subtraction process has worked well and the distribution is consistent with 
zero. In the case where no Monte Carlo truth signal event selection has been performed (points 
with errors) a clear peak is again seen in the vicinity of the chargino mass, with few events 
at higher values. For 100 fb^^ the statistical significance of the peak is around 3 a indicating 
that more integrated luminosity (or an improved event selection) would be required to claim a 
5 a discovery. Nevertheless it seems reasonable to claim that if this data were generated by 
an LHC experiment such as ATLAS, and that the observed signal were indeed not a statistical 
fluctuation, then the mass of the lightest chargino could be measured to a statistical precision 
~ ± 25 GeV (~ 11 %). More work is needed to determine the likely systematic error in 
this quantity arising from effects such as the statistical and systematic uncertainty in the input 
sparticle masses used when calculating the Xi momentum and xf mass. 

More work is needed to identify the optimum set of selection criteria required to identify 
hadronic decays in this sample, with the efficiency of the tau veto (required to remove xf 
decays via fii'r) in particular needing to be optimised. Possible methods for selecting the correct 
lepton mapping used to calculate the Xi momentum also deserve further study. With these 
improvements and/or more integrated luminosity it should then be possible both to increase the 
accuracy of the chargino mass measurement and to study quantities such as the helicity of the 
xf through measurement of the invariant mass distribution of the and the hard jet produced 
alongside the xf in the decay of the parent qi. 

5. CONCLUSIONS 

A study of the identification and measurement of charginos decaying to H^^Xi produced at 
the LHC has been performed. The results indicate that for one particular mSUGRA model the 
mass of the xf can be measured with a statistical precision ~ 1 1 % for 100 fb^^ of integrated 
luminosity. 
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Part VIII 

Chargino/Neutralino Sector In Combined 
Analyses At LHC/LC 

K. Desch, J. Kalinowski, G. Moortgat-Pick, MM. Nojiri and G. Polesello 
Abstract 

We demonstrate how the interplay of a future e^e~ LC at its first stage 
with ^/s < 500 GeV and of the LHC could lead to a precise determina- 
tion of the fundamental SUSY parameters in the gaugino/higgsino sec- 
tor without assuming a specific supersymmetry breaking scheme. The 
results are shown for the benchmark scenario SPSla, taking into ac- 
count realistic errors for the masses and cross sections measured at the 
LC with polarised beams and mass measurements at the LHC. 

1. INTRODUCTION 

The unconstrained MSSM has 105 new parameters and SUSY analyses at future experiments, 
at the LHC and at a future Linear Collider (LC), will have to focus on the determination of 
these parameters [91,92]. An interesting possibility to explore SUSY is to start with the gaug- 
ino/higgsino particles which are expected to be among the lightest SUSY particles. At tree 
level, this sector depends only on 4 parameters: Mi, M2, /i and tan/3 - the U(l) and SU(2) 
gaugino masses, the higgsino mass parameter and the ratio of the vacuum expectations of the 
two Higgs fields, respectively. 

Some strategies have been worked out for the determination of the parameters M2, Mi, 
/i, tan/5 even if only the light gaugino/higgsino particles, x^, X2 and xf were kinematically 
accessible at the first stage of the LC [93,94]. In this contribution we demonstrate how such an 
LC analysis could be strengthened if in addition some information on the mass of the heaviest 
neutralino from the LHC is available. We consider the cases: (i) stand alone LC data and 
(ii) joint analysis of the LC and the LHC data. The results in the last scenario will clearly 
demonstrate the essentiality of the LHC and LC and the benefit from the joint analysis of their 
data. 

We take the SPSla as a working benchmark [30,45] and assume that only the first phase 
of a LC with a tunable energy up to = 500 GeV would overlap with the LHC running. 
Furthermore, we assume an integrated luminosity of J £ ~ 500 fb^^ and polarised beams with 
P(e~) = ±80%, -P(e"*") = ±60%. In the following ctl will refer to cross sections obtained with 

P{e ) = -80%, P(e+) = ±60%, and an with P{e~) = ±80%, P(e+) = -60%. 

2. THE GAUGINO/HIGGSINO SECTOR 

The mass matrix Aic of the charged gaugino W"^ and higgsino H"^ depends on M2, /i, tan/5. 
The mass eigenstates are the two charginos Xi2- For real Aic the two unitary diagonalisation 
matrices can be parameterised with two mixing angles $i The mass eigenvalues m^± and 

the mixing angles are analytically given by the Susy parameters (see e.g. [95,96]). The cross 
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section a^jij} = a{e^e — xfxf) can be expressed as a function of (cos 2$^ ^j, m^±); the 

' Xi 

coefficients for (T^{11} are explicitly given in [97]. 

The neutralino mixing matrix A4n depends on Mi, M2, ^ and tan/3. Analytic expres- 
sions for the mass eigenvalues m?o and the eigenvectors are e.g. given in [93, 94]. The 

-^1, 4 

characteristic equation of the mass matrix squared, M. ArTWjv, is written explicitly as a quadratic 
equation for the parameter Mi [97]. 

3. STRATEGY FOR THE DETERMINATION OF THE SUSY PARAMETERS 

At the initial phase of future e"'~e~ linear-collider operations with polarised beams, the col- 
lision energy may only be sufficient to reach the production thresholds of the light chargino 
xf and the two lightest neutralino s Xi, X2- Nevertheless the entire tree level structure of the 
gaugino/higgsino sector can be unraveled [93-95]. 

Chargino cross sections measured at ^/s = 400 GeV and 500 GeV with polarised beams 
and the lightest chargino mass are sufficient to determine unambiguously the mixing angles 
cos2$i R. Then the Mi can be obtained from the quadratic equation TWjvA^jy. However, 
using the kinematically accessible cross sections for the neutralino production ^{12} and 
a1 ^{22} leads to a precise determination of the fundamental Susy parameters [97]. 

In the following we perform this strategy for the benchmark scenario SPSla [30,45] 
defined at the electroweak scale: Mi = 99.13 GeV, M2 = 192.7 GeV, = 352.4 GeV, 
tan (3 = 10; the resulting masses are given in tabled 

3.1 SUSY PARAMETERS FROM THE LC DATA 

We use the light chargino and neutralino masses fn^±, m^o ^ and the polarised cross sections for 
the processes e~^e~ xtXi^ X1X2' X2X2 = 400, 500 GeV as experimental input. 

In our scenario the light chargino Xi ^^id also the neutralino X2 decay mainly via fi chains 
producing the final states similar to that of stau pair production, however with different topology. 
Therefore, we assume that the contamination of stau production events can be subtracted from 
the chargino and neutralino production [97]. 

In our analysis we take the production cross sections with statistical errors induced by the 
following uncertainties: 

• The chargino mass measurement has been simulated and the expected error is 0.55 GeV, 
see tabled 

• With J C = 500 fb^^ at the LC, we assume 100 fb~^ per each polarisation configuration 
and we take into account la statistical error. 

• Since the chargino (neutralino) production is sensitive to ni;^^ {rrie^ ^), we include the 
experimental error of their mass determination of 0.7 GeV (0.2 GeV, 0.05 GeV), see 
tabled 

• Concerning the neutralino cross sections we estimate the statistical error based on an ex- 
perimental simulation^ yielding an efficiency of 25% and include an additional systematic 
error {Sa-^g) which takes into account the uncertainty in the background subtraction, for 
details see [97]. 

^M. Ball, diploma thesis. University of Hamburg, January 2003, http://www- 
flc.desy.de/thesis/diplom.2002.ball.ps.gz. 
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Xf X2 


Xl X2 X3 Xi 


eR Cl i>e 


mass [GeV] 
error [GeV] 


176.03 378.50 
0.55 


96.17 176.59 358.81 377.87 
0.05 1.2 


143.0 202.1 186.0 
0.05 0.2 0.7 



Table 1: Chargino, neutralino and slepton masses in SPSla, and the simulated experimental errors at the LC^. 
It is assumed that the heavy chargino and neutralinos are not observed at the first phase of the LC operating at 
y/s < 500 GeV. 





400 GeV 




500 GeV 




(P(e-),P(e+)) 


(-80%, +60%) (+S 


0%, -60%) 


(-80%, +60%) (+S 


0%, -60%) 


a{e+e- xtXi) 


215.84 


6.38 


504.87 


15.07 


^^^total 


7.27 


0.35 


5.28 


0.51 


a{e+e- X1X2) 


148.38 


20.06 


168.42 


20.81 


^^total 


3.0 


1.58 


3.52 


1.57 


a{e+e- X2X2) 


85.84 


2.42 


217.24 


6.10 


^^^total 


3.6 


0.41 


4.3 


0.62 



Table 2: Cross sections a^j^ill} = o-L,/?,(e+e xfx^ ), cr°,ij{12} = (TL.i?(e+e X1X2) and 

CT^_^{22} = crL,R{e+e- X2X2) with polarised beams P{e-) = t80%, P(e+) = ±60% at = 400 
and 500 GeV and assumed errors (in fb) corresponding to 100 fb~^ for each polarisation configuration. 



• The beam polarisation measurement is assumed with an uncertainty of AP(e^) / P(e^) = 
0.5%. 

The resulting errors are listed in table |2| 

From (xfixtXi) at = 500, 400 GeV and a^ixtXi) at = 500 GeV exploiting the 
relation cos2$i(; = /(cos2$l, cr^^jll}) we first predetermine chargino mixing angles as 

cos2$i = [0.62,0.72], cos 2$^, = [0.87,0.91] (1) 

Then using the neutralino cross sections cr^(XiX2)' '^l (X2X2) at \/s = 500, 400 GeV and light 
neutralino masses m^o ^ within their experimental errors, a rather accurate determination of the 

SUSY parameters can be obtained from the Ax^ test defined as Ax^ = J2i K^j ~ Oi)/60i\'^. 
The sum over physical observables Oi includes m^o,m^o and neutralino production cross sec- 
tions 0"° ^{12}, cr° ^{22} measured at both energies of 400 and 500 GeV, Oi stands for the 
physical observables taken at the input values of all parameters, and 50i are the correspond- 
ing errors. The Ax^ is a function of unknown Mi, cos2$l, cos2$/j with cos 2$^, cos 2$/? 
restricted to the ranges given in eqn. dl]) as predetermined from the chargino sector. In fig.Ht 
the contour of A^^ = 1 is shown in the Mi, cos 2$/,, cos parameter space along with its 
three 2dim projections. The projection of the contours onto the axes determines \a errors for 
each parameter. 

Values obtained for Mi, cos 2$^, cos 2$^ together with m-± can be inverted to derive the 
fundamental parameters M2, fi and tan (3. At the same time the masses of the heavy chargino 
and neutralinos are predicted, see table|3l As can be seen in table|3l the parameters Mi and M2 
are determined at the level of a few per-mil. The is reconstructed within a few per-cent, while 
for tan/? the error is of order 15%. 

^H.U. Martyn, LC-note LC-PHSM-2003-071. 
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Figure 1: The A^^ = 1 contour in the {Mi, cos2^l,cos2^]i} parameter space derived a) from the LC data and 
b) from the joint analysis of the LC data and LHC data [97]. 





SUSY Parameters 

Ml Ma /i tan/? 


Mass Predictions 

^xi "^xi "^x2 


LC 


99.1 ±0.2 192.7 ±0.6 352.8 ± 8.9 10.3 ±1.5 


378.8 ± 7.8 359.2 ±8.6 378.2 ±8.1 


LC/LHC 


99.1 ±0.1 192.7 ±0.3 352.4 ±2.1 10.2 ± 0.6 


378.5 ±2.0 358.8 ±2.1 



Table 3 : SUSY parameters with 1 a errors derived from the analysis of the LC data and from the combined analysis 
of the LHC and LC data (with 5m^a — 0.08 GeV and Sm^a — 2.23 GeV derived from the LHC when using the 
LC input of Sm^o ~ 0.05 GeV) collected at the first phase of operation. Shown are also the mass predictions of 
the heavier chargino/neutralinos. 



3.2 COMBINED STRATEGY FOR THE LHC AND LC 

A rather large error of the /i parameter derived above stems from the gaugino-dominated char- 
acter of the light charginos/neutralinos in the SPSla scenario. A significant improvement for 
the /i (and tan /3) is expected with additional information on heavy neutralinos available from 
the LHC. 

The LHC will provide a first measurement of the masses of Xi, X2 Xi^ see^. The 
measurements of X2 Xa will be achieved through the study of the processes Xi ^ ^ 
iiXi, (with i = 2, 4) in which the invariant mass of the two leptons in the final state shows an 
abrupt edge at (m™*^)^ = m|o(l — m|/m^o)(l — m^o/mj). 

With the LHC data, the achievable precision on m^o and m^o is expected be respectively 
4.5 and 5.1 GeV for an integrated luminosity of 300 fb^^. However, since the uncertainty 
on the m^o and m^o depends both on the endpoint determination and also on m^o and m^, 
a much higher precision can be achieved with m^o, rrier, and mg measured at the LC with 
precisions respectively of 0.05, 0.05 and 0.2 GeV, table[T] With this input the precisions on the 
LHC+LC measurements of m^o and m^o become: 5m^o = 0.08 GeV and Sm^o = 2.23 GeV 

A2 A4 A2 A4 

Performing again the A^^ test with this additional input one gets a significant improvement in 
the accuracy, see fig. [1J5 and table |3] for the final results. The accuracy for the parameters /i 
and particularly tan P is now much better, better than from other SUSY sectors [98-100] (and 
references therein). 

^B.K. Gjelsten, J. Hisano, K. Kawagoe, E. Lytken, D. Miller, M. Nojiri, P. Osland, G. Polesello, contribution 
in the LHC/LC working group document, see also http://www.ippp. dur.ac.uk/georg/lhclc/. 
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4. SUMMARY 

We have worked out in a specific example, an mSUGRA scenario with rather high tan /3 = 10, 
how the combination of the results from the two accelerators, LHC and LC, allows a precise de- 
termination of the fundamental SUSY parameters without assuming a specific supersymmetry 
breaking scheme. We have shown that a promising hand-in-hand procedure consists of feeding 
the LSP and slepton masses from the LC to the LHC analyses and injecting back a precise ex- 
perimental determination of the X2 xl masses. It provides a determination of Mi, M2, /i 
at the < 0(1%) level and of (rather high) tan P of the order of < 10%, reaching a stage where 
radiative corrections become relevant in the electroweak sector and which will have to be taken 
into account in future fits. 
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Part IX 

Proposal For A New Reconstruction 
Technique For Susy Processes At The 
LHC 

M.M. Nojiri, G. Polesello and D.R. Tovey 
Abstract 

When several sparticle masses are known, the kinematics of SUSY de- 
cay processes observed at the LHC can be solved if the cascade decays 
contain sufficient steps. We demonstrate four examples of this full re- 
construction technique applied to channels involving leptons, namely 
a) gluino mass determination, b) sbottom mass determination, c) LSP 
momentum reconstruction, and d) heavy higgs mass determination. 

1. INTRODUCTION 

The potential of the LHC for SUSY parameter determination has been studied in great detail for 
the past seven years. One of the most promising methods involves the selection of events from 
a single decay chain near the kinematic endpoint. Information on the masses involved in the 
cascade decay can be extracted from the endpoint measurements. It has been established that 
one can achieve a few percent accuracy for sparticle mass reconstruction using this technique 
with sufficient statistics. 

In this paper we propose a new method for reconstructing SUSY events which does not 
rely only on events near the endpoint. Instead one kinematically solves for the neutralino mo- 
menta and masses of heavier sparticles using measured jet and lepton momenta and a few mass 
inputs. 

To illustrate the idea we take the following cascade decay chain 

g^bb^ xlbb £bb£ Xibbii. (1) 

This decay chain is approximately free from SM background with appropriate cuts. The five 
SUSY particles which are involved in the cascade decay have five mass shell conditions; 

{Pxl + Pei + Pfe + Pfci + Pte ) ^ • (2) 

Of these five masses, m^o,m^o and can be measured at the LHC using first generation 
squark cascade decays with an accuracy of ~ 10% (the mass difference is measured more 
precisely). Moreover, with input from a future high energy Linear Collider these masses might 
be determined with an accuracy ~ 0(1%). We therefore assume for the present work that the 
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masses of the two lighter neutralinos and of the right handed slepton are known, and we ignore 
the corresponding errors. 

For a hMi event, the equations contain six unknowns {nig, and p^o) which satisfy five 
equations. For two hM^ events, we have ten equations while we only have ten unknowns (two 
neutralino four momenta, nig and m^o ). Mathematically, one can obtain the sbottom and gluino 
masses and all neutralino momenta if there are more than two hhU events. 

We call this technique the "mass relation method" as one uses the fact that sparticle masses 
are common for events which go though the same cascade decay chain. Note events need not be 
near the endpoint of the decay distribution to be relevant to the mass determination. In the next 
section we demonstrate the practical application of this method to measurement of the masses 
of the gluino and sbottom. 

As a byproduct of the technique, once the mass of the squark and of all the sparticles 
involved in the decay are known, the momentum of the lighter neutralino can be fully recon- 
structed, and this further constrains the event. 

In SUSY events sparticles are always pair produced and there are two lightest neutralinos 
in the event. If squark decays via g — > X2 ~^ ^ Xi can be identified on one side of the 
event then the neutralino momentum can be reconstructed as described above. The transverse 
momentum of the lightest neutralino in the other cascade decay can then be obtained using the 
following equation 

Pt(x?(2)) = PT(miss) + Pt(x?(1)), (3) 

provided that there are no hard neutrinos involved in the decay. This transverse momentum can 
be used to constrain the cascade decay of the other sparticle. 

For the case where the other squark decays via q ^ Xi(l ^ followed hy W ^ 

q'q" , the chargino mass can be determined by using Eq. © and the following relations, 

P] = ml (4) 

where pj is the momentum of the selected high px jet which comes from the squark decay and 
Pw is the momentum of the two jet system consistent with the W interpretation. The neutralino 
momentum resolution is important for the chargino mass reconstruction and we discuss this in 
section 3. The reconstruction will be discussed in a separate contribution [101]. 

The full reconstruction technique can be extended for higgs mass reconstruction. In sec- 
tion 4, we discuss the heavy higgs mass determination from the process H — * X2X2 followed 
by X2 ^ ^ ^^Xv This process is also useful for discovery of heavy higgs bosons. The four 
lepton momenta and missing momentum can be used to reconstruct the higgs mass assuming 
that the Pt of the higgs boson is very small. 

2. GLUINO CASCADE DECAY 

We first discuss the results of a simulation study of the process where a gluino cascade decays 
into a sbottom at model point SPS la [30]. The relevant sparticle masses for this study are listed 
in Table 1. The events were generated using the HERWIG 6.4 generator [12] [87] and passed 
through ATLFAST [88], a parametrised simulation of the ATLAS detector. 

We study only events which contain the cascade decay shown in Eq.©. We then apply 
the following preselections to reduce backgrounds: 
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m-g 


"^61 (2) 


A,2 




Al 


595.2 


491.9(524.6) 


176.8 


136.2 


96.0 



Table 1: Some sparticle masses in GeV at SPSla. 



• pjp''' > 100 GeV 

• Meff > 600 GeV 

• at least 3 jets with pn > 150 GeV, pt2 > 100 GeV and prs > 50 GeV. 

• exactly two jets with > 50 GeV tagged as 6-jets 

• exactly two OS-SF leptons with pth > 20 GeV, PT12 > 10 GeV, and invariant mass 
40GeV< mil < 78 GeV. 

The solution of Eq. Q can be written in the following form: 

ml = Fo + Fiml±F2D, 
where = + Am? + Dgmf. (5) 

Here Fi and Di depend upon pe^ and ph. and the neutralino and slepton masses. In the event, 
there are two b jets and we assume that the b jet with larger pr originates from the b decay. 
The two leptons must come from X2 and i decay. There are maximally four sets of gluino and 
sbottom mass solutions together with two lepton assignments for each decay, because we cannot 
determine from which decay the lepton originates. To reduce combinatorics we take the event 
pair which satisfies the following conditions: 

• Only one lepton assignment has a solution to the Eq. Q 

• For a pair of events there are only two solutions and there is a difference of more than 
100 GeV between the two gluino mass solutions. 



xVndf 41.79 / 23 




m{gluino) (GeV) 



Figure 1: rrig obtained by using Eq. (|5} for two bbi£ events. 
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In Fig 1, we plot the minimum rrig solution which satisfies the conditions given above. 
The peak position is consistent with the gluino mass, and the error on the peak position obtained 
by a Gaussian fit is around 1 .7 GeV for 100 fb^^. For the events used in the reconstruction, each 
event is used on average five times. Note that the a of the Gaussian fit is large (~ 56.7 GeV) 
and is determined by the resolution on the momentum measurement of the four 6-jets. It is 
worth stressing that the results presented here were produced by using a parametrised simulation 
of the response of the ATLAS detector to jets, based on the results of a detailed simulation. 
Results which crucially depend on the detailed features of the detector response, such as the 
possibility of discriminating the two sbottom squarks (see below) need to be validated by an 
explicit detailed simulation of the detector performed on the physics channel of interest. We 
only attempt here to evaluate the impact of the new technique on sparticle reconstruction. 

Once the gluino mass has been determined one can reconstruct the sbottom mass by fixing 
the gluino mass to the measured value. Here one need only solve Eq.©, which involves only 
two 6-jets in the fit, and therefore errors due to the jet resolution are expected to be less than 
those for the gluino mass reconstruction. 

For each event, there are two sbottom mass solutions mj(soll) and m^(sol2), each sensi- 
tive to the gluino mass input. The difference between the gluino and sbottom mass solutions is 
however stable against variation in the assumed gluino mass. The mass itself may have a large 
error in the absolute scale, but the mass differences are obtained rather precisely, as is the case 
in the endpoint method. 

In Fig. 2 (left), we plot the solutions for all possible lepton combinations in the nig — 
(soil) rrig — (sol2) plane. Here we use the fe-parton momentum obtained from generator 
information. One of the solutions tends to be consistent with the input sbottom mass. Moreover 
the two decay modes g bib and are clearly separated. 

We can compare the results from the previous analysis with those from the endpoint anal- 
ysis [85], where one uses approximate the formula 

This formula is correct only at the endpoint of the three body decay ^2 ~^ Xi^^^ but is never- 
theless approximately correct near the edge of ^2 ^ ^-^ -^^Xi for SPSla. The sbottom mass 
obtained by using Eq.® is shown in Fig. 2(right). For this case, the 62 peak at 70.6 GeV is not 
separated from the 61 peak at 103 GeV. 

The bi mass, or the weighted average of the sbottom masses, is easily obtained. The b 
jet resolution is not sufficient however to clearly separate the bi and 62- This can be seen in 
Fig. 3 where the plots show the distributions corresponding to Fig.2(left) and (right) but now 
with the b parton momenta replaced by b jet momenta. For the endpoint analysis (Fig. 3 right), 
a correct evaluation of the sbottom masses would require a fit taking into account the shape of 
the response of ATLAS to b-jets. In order to approximately evaluate the achievable statistical 
precision, a naive double gaussian fit was performed on the distribution shown in Fig. 3 right, 
which corresponds to / dtL = 300 fb^^. The resulting statistical uncertainties are ±1 GeV 
(±2.5 GeV) for the nig — mi_^(mg — nii^ ) peak positions respectively. Additional systematic 
uncertainties, not yet evaluated, as well a 1% error due to the uncertainty on the jet energy scale 
should also be considered. These numbers are obtained assuming the presence of two gaussian 
peaks in the data. 

For the mass relation method the number of events available for the study is larger by a 
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m(gluino)-m(sbottom,sol1) mgl-msb(GeV) endpoint approximation(parton level) 

Figure 2: The distribution of rrig — mi calculated using the parton level h momentum by solving Eq.(|3 (left) and 
using the approximate relation Eq. |6jright). 



factor of 2 because events away from the endpoints can be used. We also use the exact formula 
for the mass relation method. Although the analysis is more complicated due to the multiple 
solutions, we believe it to be a worthwhile technique for use when attempting to reconstruct the 
hi and 62 masses. 



3. NEUTRALINO MOMENTUM RECONSTRUCTION 

In this section, we discuss the reconstruction of the momentum of the lightest neutralino. As we 
have discussed already, the mass shell condition can be solved for long decay cascades, such as 
9 ~^ X2Q ~^ ~^ XiQ^^- For this process we have two neutralino momentum solutions for 
each lepton assignment. One may wonder if the solutions for the neutralino momentum might 
be smeared significantly, because of the worse jet energy resolution as compared to leptons, and 
the jet pt is generally much larger than the neutralino momentum for the cascade decay. In Fig. 
4(left) we show the distribution of pT(reco)/pT(tmth) for the point studied in [101]. Here we 
choose the correct lepton combination using generator information, and take the solution which 
minimizes |]?T(reco)/pT(truth)— 1|. Except for the case where we took the wrong jet as input 
the reconstructed pr is within 20% of the true neutralino momentum. The result for the gluino 
cascade decay into sbottom Eq.(l) is similar. 

In Fig. 4(right) we show a similar reconstruction for the gluino cascade decay, but unlike 
Fig.4(left), we use both lepton combinations. We fix the gluino mass to the input value"* and 
take events where one of the four sbottom mass solutions is consistent with the input sbottom 
mass such that |m^^ — m|;(best) | < 10 GeV. We then take the solution where the sbottom mass is 
closest to the input m^^ . There are still two p^o solutions, and we choose the one which minimize 
|min(pT(reco)/pT(truth),pT(truth)/pT(reco)) — 1|. The neutralino momentum resolution is 
worse than that obtained using the correct lepton assignments only. Nevertheless a significant 
fraction of events are reconstructed with 0.8 < |p(reco)/p(truth)| < 1.2. 

'^Here we adopt an event selection which makes use of the true (input) gluino and sbottom mass values, although 
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m(gluino)-m(sbottom,sol1) m(xbb)-m(xb) (GeV) 



Figure 3; As for Fig. 2 but with the 6 jet momentum used instead of the b parton momentum. 



4. HIGGS MASS RECONSTRUCTION 

A promising decay for the observation of heavy and pseudo-scalar higgs bosons in the difficult 
region with intermediate tan /5 is the decay into two neutralinos. When both neutralinos decay 
through the chain 

the resulting signature consists of events with four isolated leptons (paired in opposite-sign 
same-flavour pairs) and no jet activity. The main SM backgrounds to this signature are ti pro- 
duction, where both the 6-jets and the Ws decay into leptons and Zbb production. The key 
element for the rejection of these backgrounds is the fact that the leptons from b decays are not 
isolated. A detailed study of the performance of lepton isolation in the detector is needed to 
assess the visibility of the signal. Additionally there is an important SUSY background, includ- 
ing irreducible backgrounds from direct slepton and gaugino decay. Full background studies 
as a function of the SUSY parameters were performed by the ATLAS and CMS Collabora- 
tions [79, 102]. We propose here, along the lines of the previous sections, a technique for the 
complete reconstruction of the higgs peak, based on the knowledge of the masses of 
and Xi- In this case one has 8 unknown quantities: the 4-momenta of the two LSP's, and 8 
constraints: six on- shell mass constraints (3 for each leg), and the two EJp'^'^^ components. 

To demonstrate the power of the method, we apply it to Point SPSla, for which the mass 
of the A and of the i/ is ~ 394 GeV. The BR into X2X2 is 6% (1%) for the A(H). We perform 
the study on 1000 events for 

corresponding approximately to the expected statistics for 300 fb^^ We simply require 2 iso- 
lated leptons with pt > 20 GeV and 2 further isolated leptons with pt > 10 GeV, all within 
\r]\ < 2.5. The efficiency of these cuts is ~ 60%. 

We have not performed any background simulations because at this stage we only wish to 
explore the viability of the full reconstruction technique. The main problem for the reconstruc- 



in practice fitted values would be used. 
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p-,-(reco)/Pj(truth) |p|(reco)/|p|(truth) for the best solution 

Figure 4: Left; The calculated Xi transverse momentum divided by the true transverse momentum. The de- 
cay ql iqi XiQ^^ is studied for the model point used for the chargino study: toq = 100 GeV, 
"^1/2 = 300 GeV, Aq = —300 GeV, tan/3 = 6, and yu > 0. Only the correct lepton choice is used. Right: 
|p(reco)/p(truth)| for the decay chain Eq.Q for SPSla. 



tion is the correct assignment of the leptons to the appropriate decay chain. The first selection 
is based on requiring a unique identification of the lepton pairs coming from the decays of the 
two X2^- We therefore require that either of the following two criteria is satisfied: 

• the flavour configuration of the leptons is e+e^/i+zi" 

• the lepton configuration is either e+e'e+e^ or but for one of the two possible 
pairings the invariant mass of one of the pairs is larger than 78 GeV, i.e. above the lepton- 
lepton edge for the X2 decay. 

The total efficiency after these cuts is ~ 30%. At this point, on each of the two legs there is 
still an ambiguity due to the fact that each lepton can be either the product of the first or of 
the second step in the decay chain. This gives 4 possible combinations. Furthermore, the full 
reconstruction results in a quartic equation which can have zero, two or four solutions. We show 
in Fig. ^ the distribution of the calculated A mass for all of the retained combinations as a full 
line. The dashed line shows the combinations with the wrong lepton assignment. A clear and 
narrow peak emerges over the combinatorial background. The width is approximately 6 GeV, 
determined by the resolution of the measurement of the momentum of the leptons. 

5. CONCLUSIONS 

In this contribution, we have described a novel technique for reconstructing the mass and mo- 
menta of SUSY particles. This technique does not rely on any approximate formulae nor on 
endpoint measurements. All events contribute to the sparticle mass determination and decay 
kinematics reconstruction, even if they are away from the endpoint of the distribution. The 
method may be particularly useful when the SUSY mass scale is large. In that case the statistics 
can be so low that the endpoint cannot be seen clearly while the SUSY sample itself is very 
clean. 
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Figure 5: Reconstructed mass of the A. 



The method applies most effectively when we know some of the sparticles' masses ex- 
actly, because the number of unknown parameters in e.g. Eq.© is reduced. In the particular 
case where some of the sparticle masses are measured at a LC the sparticle cascades may be 
solved completely and study of the decay distributions and higher mass determination becomes 
possible at LHC. 

When all the sparticle masses are known the neutralino momentum can be reconstructed if 
four sparticles are involved in the cascade decay. The sparticles would be pair produced, and if 
we can identify both of the cascade decay chains in the events then we only need six sparticles in 
the cascade decay to solve both of the neutralino momenta on account of the missing momentum 
constraint. The reconstruction of sparticle momenta provides us with an interesting possibility 
for studying the decay distribution at the LHC. 

On the other hand, our method is not valid when some of the particles in the cascade 
produce hard neutrinos. This is unfortunately the case when the chargino decays into (s)leptons, 
when a f is involved in the decay, or when a 11/ is produced and decays leptonically. If such 
SUSY decay processes dominate then this method may not be useful. 
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PartX 

Building on a Proposal for a New 
Reconstruction Technique for SUSY 
Processes at the LHC 

C.G. Lester 

Abstract 

There has recently been interest in "a new reconstruction technique for 
SUSY processes at the LHC". The primary intention of this note is to 
describe a modification to the way the technique is used. This modifi- 
cation suggests the method is much more powerful than originally pro- 
posed. We demonstrate that, in principle at least, the method does not 
need to rely on input from other experiments. We show that the method 
is capable of standing on its own, and is able to measure the masses 
of all the sparticles participating in the relevant decay chains. Results 
from other experiments such as a future linear collider may easily be 
incorporated if desired. 

1. INTRODUCTION 

The authors of [103] propose "a new reconstruction technique for SUSY processes at the LHC". 
Their method is described in detail in their article elsewhere in these proceedings, and so only 
an outline of their method will be provided here. The reader is strongly encouraged to read their 
article before reading this one. 

The authors of [103] refer to their technique as the "mass relation method". I wish to 
narrow the meaning of this phrase, as I want to draw a distinction between (1) the idea that 
makes the whole method work, and (2) any particular implementation of that idea. I will use the 
phrase "mass relation method" to describe any method which, for its success, is forced to rely 
on the extraction of information from two or more independent events which are related only 
by their sharing similar or identical particle content. 

2. OUTLINE OF ORIGINAL IMPLEMENTATION 

In [103] it is suggested that "when several sparticle masses are known, the kinematics of SUSY 
decay processes observed at the LHC can be solved if the cascade decay contains sufficient 
steps". Therein, an original implementation is described which accomplishes this solution. The 
efficacy of this implementation is demonstrated using the decay chain 

g^bb^ xlbb Ibbl XibbU. (1) 

In order to complete the demonstration, the authors make the assumption that "... the masses 
of the two lighter neutralinos and the right handed slepton are known, and [they] ignore the 
corresponding errors". By ignoring the corresponding errors, they in effect demonstrate their 
method in a scenario in which the slepton and neutralino masses are already known to something 
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like the one-percent level - something not usually assumed to be possible with LHC data alone. 
Having made these assumptions, the authors set out to measure the masses of the two remaining 
sparticles. Firstly, on an event-pair by event-pair basis, they obtain an estimate of their gluino 
mass with an accuracy of about ±60 GeV (approximately 10%). By histogramming the results 
from all these event-pairs an overall measurement of the gluino mass is obtained with an error 
of a few GeV. This measurement is then fed back into the events which are analysed in a second 
pass in order to make a measurement of the sbottom mass.^ 

2.1 A small part of the original implementation in more detail 

As already mentioned in Section [Lj the authors begin by assuming prior knowledge of m^o, 
nil "^X2 ' ^^^^ masses of the two lightest neutralinos and the mass of the slep- 

ton participating in the decay chain shown in Equation ([T]). This leaves in any one event the 
six unknown real quantities comprising: uig, rrii and the four components of p^o. These six 
unknowns, are however constrained to satisfy the five mass constraints: 
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in which pi-^, pi^, p^^ and pi,.^ are the four-momenta of the emitted standard model particles.^ 
Since the number of unknowns (six) exceeds the number of constraints (five) it is not possible 
to conclude much from one event. 

Taking a second event together, however, the number of unknowns rises by only four, 
namely the four components of the x^'^norn^iitum in the new event. As ever, the new event, 
like the old, will satisfy another five mass constraints of the form shown in Equation Q. With 
two events, then, the number of unknowns and number of constraints have each risen to ten. 
So in principle, with only two non-degenerate events, it is now possible to determine all the 
unknowns. This amounts to full reconstruction of both events and determination of rrig and 
m^. The interested reader is directed to [103] to see how the authors handle choice-ambiguities 
that arise from (a) the solution of simultaneous quartic and quadratic equations, and (b) lack of 
knowledge of which of the two observed leptons is li and which of the two observed 6-tagged 
jets came from bi. 

The implementation of [103] proposes that one should do exactly as described above: 
namely consider events in pairs. 

3. MOTIVATIONS FOR BUILDING ON THE ABOVE 

A natural reaction on seeing the original implementation is to ask: 

^In practice, as in most of these analyses, the mass difference between the sbottom and the gluino is measured 
more accurately than the absolute value of either of their masses. 

^Strictly speaking these mass constraints apply only to the true rather than the measured momenta of the emitted 
standard model particles. However, for the purposes of the "original implementation" this distinction did not need 
to be drawn, and the measured momenta were used "as if" true. The resulting smearing of the answer was accepted 
as a source of reconstruction error. There are differences between this method and that of my proposal. 
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Since two events are better than one, why not consider even more? 

The majority of the rest of this note tries to address the above question. The motivation for 
building on the above is that every new event adds five more constraints but only four more un- 
knowns. Put another way, for every additional event that is acquired, one can either answer one 
new question, or else better constrain any answers that one already has. This note concentrates 
on the last of these two possibilities. 

4. PROPOSAL FOR A "NEW IMPLEMENTATION" OF A MASS RELATION METHOD 

4.1 General comments 

The new proposal is to do nothing more than consider all the relevant events simultaneously. 

In this note, we will not address the important question of whether, in a real LHC ex- 
periment, it would be possible to satisfy the preconditions for the success of any mass relation 
method, namely the requirements that it be possible to construct sufficiently pure samples of 
the appropriate standard model samples from chains of sufficiently similar or identical particle 
content. This needs to be addressed in further papers, and has already been considered in part 
by [103]. The intention of this note is only to look at what may be achieved ?/such selection 
were possible. 

4.2 Detailed description of proposal 

Ideally we would like to know the masses of the sparticles in our events. Realistically, we can 
only expect to find the masses within some finite precision or error. Bearing correlations in 
mind, the best we can expect to determine is the relative probability of any particular combina- 
tion of the five masses m = (m^o, m^-, m-j^o, m^, nig) given the data. In short, we would like to 
plot p(m I data). 

By Bayes' theorem, j9(m|data) oc p(data|m)po(m)- The first factor, the likelihood, will 
be determined purely from the events considered and the mass relation method itself and is 
thus the objective "result" of this experiment. The last factor, the prior, incorporates all exist- 
ing knowledge gained from other experiments (if you should wish to include them) and any 
subjective preferences you might have. Because we choose here to use a non-informative prior 
(uniform in the hierarchical sparticle masses^) the reader may view the results at his or her dis- 
cretion as either (a) simple plots of the objective likelihood distribution ]?(data|m), or else, (b) 
indicative of the results which a single experiment would provide p(m|data) in the absence of 
data from other experiments. 

The only thing remaining to be defined is p(data|m). As the data consists of many inde- 
pendent events, we have j» (data |m) = p(eventi|m).p(event2|m). ■ ■ ■ .p(event„|m), and we are 
left needing to evaluate p(eventi|m) for a given event i. 

Evaluating the event-likelihood p(eventj|m) properly and efficiently is the hard part. In 
principle there is only one right answer, which you would obtain by taking the square of the 
matrix element for the observed final state and integrating it over all the unknowns in the prob- 
lem, namely the measurement errors and the momentum distribution of the unobserved chain 
progenitor. The answer you obtain will thus depend on which model assumptions you wish to 
make, e.g. whether you would choose to model the differing spins of the sparticles. 

^Additionally, the cosmetic constraints 30 GcV < m-^o < 300 GeV and rrig < 1000 TcV are incorporated 
into our prior in order to frame all plots nicely. 
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In order to meet the time constraints imposed by the submission of this note, however, it 
was necessary to implement the event-likelihood using an ad-hoc approximation.^ It is hoped 
that the approximation to the event-likelihood described later is sufficiently similar to the full 
form of p(eventi|m) that the basic features of the proposed technique can be demonstrated. 

We construct the approximation to the event-likelihood used in this note as follows. For 
a given chain momentum hypotheses H consistent with the given mass hypothesis m, we can 
define a "distance" S{H) between the observed and the true momenta of the visible particles 
produced in that chain by 

6iH) = ef,/af + ef,/af + el^/af + eln/al (3) 

where (for example) efn/af is square of the number of standard deviations by which the mea- 

sured momentum of li differs from a hypothesised true value pf^. 

We can then perform a least-squares minimisation of 6{H) over all possible chain momen- 
tum hypotheses H consistent with the given mass hypothesis m. We now make the assumption 
that this least-squared minimisation will have provided us with the momenta which are most 
consistent with the observation and the mass hypothesis m. Finally, then, we can approximate 
the event-likelihood by the simple process of evaluating the probability for the observed mo- 
menta to deviate as far as their observed values, assuming that the "true" momenta are given by 
the result of the fit. This approximation thus depends crucially on a good understanding of the 
measurement errors associated with the observed standard model particles. 

The key features of this approximation to the event-likelihood are that it will be large when 
p(eventj|m) is large and small when p(eventj|m) is small. ^ This gives us confidence that the 
maxima and minima of the event-likelihood will be well approximated. It will most probably 
not be the case, however, that the widths of the resulting distributions or the finer shape details 
(for a particular number of events) can be completely relied upon. This, though undesirable, is 
not too great a problem as the widths naturally scale with the number of events analysed. For 
this reason, the approximation may be thought of as resulting in an uncertainty in how many 
events are necessary to achieve a given reconstruction error, rather than an uncertainty in the 
quality of the reconstruction itself. 

5. RESULTS OF NEW PROPOSAL 

In order to demonstrate the potential of the above technique, it is necessary to generate some 
events for analysis. A toy montecarlo was used to generate these events only. For the reasons of 
Section |4~T] it simulated only relativistic kinematics and decays were according to phase space 
only. In effect all particles were treated as scalars. Furthermore, only the chain described in 
Equation ([T]) was simulated. No extraneous particles were produced, and nor is there an "other 
side of the event". 

'^A more in-depth paper currently in preparation deals with the evaluation of the full form of p(eventi |m). 

''Unlikely events (for a given m) are clearly those in which, no matter how hard you try, you find a huge 
disagreement between the momenta of the particles you see in the event and the momenta you would expect to 
have seen considering all hypothesised chain momenta that would have been consistent with the masses m. Put 
another way, you will know you have a very unlikely event when you cannot hypothesise a set of chain momenta 
in which the visible particles have momenta "close" to those observed in the event. If then we discover that 5{H) 
is large, no matter what hypothesis H we choose consistent with m, then we know that that particular event is 
unlikely given that particular m. Conversely, when m is close to the right answer, the true chain momenta H^^ 
will lead to a small value for 6{H*'^). 
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Events were simulated events using the following arbitrary values for the sparticle masses: 
m^o = 150 GeV, m^- = 200 GeV, m^o = 300 GeV, = 500 GeV and m-g = 650 GeV. 

Measurement errors were simulated by randomised rescaling of the four momenta of the 
observed particles in a manner similar to that described in [104], but without rj and (p depen- 
dence in the resolutions. Lepton momenta were smeared by 1% and jet momenta by 5%. The 
reconstruction part of the analysis was handed the two lepton momenta in a random order and 
the two 6-jet momenta in a random order so that it could not know which of them was which. 

After generating 100 events of the form described above, the results shown in Figures [T] 
were obtained. Here we choose to plot p(m|data) by sampling from it using a Metropolis 
Markov-chain sampler [105], although this particular choice is unimportant.'*' 

Figure[l]shows how the reconstructed mass distributions of the heavier sparticles correlate 
with the reconstructed mass distribution of the lightest neutralino. The correlation show that at 
this point the mass differences tend to be measured more accurately than the absolute masses. 
The mass of the lightest neutralino is well reconstructed (input value 150 GeV) with an error 
of order approximately 17%. The reconstructed mass distributions of the remaining sparticles 
are similar to that of the lightest neutralino as the mass differences are constrained better than 
the masses themselves. 

6. CONCLUSIONS 

This note aims to convey the message that should it be possible at the LHC to isolate clean 
samples of the decay products of a hundred or so sparticle decay chains containing enough 
sufficiently similar or identical particles, then by using a mass relation method of the type 
proposed, one should be able to reconstruct the masses of all the sparticles in these decay chains 
to precisions of at least some tens of percent, depending on the number of these events. 

Further work is needed to establish that the necessary sample purity is achievable, and 
to ascertain the effect that a better model for the event-likelihood would have on the widths of 
the reconstructed distributions. If the identified collections of outgoing particles are sufficiently 
pure and large, it does not seem unreasonable to believe that one might expect a precision on 
the reconstructed masses which is competitive with any other independent method found so far. 

It seems likely that the usefulness of this method will be limited by the ability to produce 
the necessary pure samples of decay-chain products.'^ 
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'^Finding good start points for the Markov-chain sampler is not a difficuh task, just a very time consuming one 
due to the high dimensionality (five) of the space to be probed. I have avoided the computational overhead of doing 
such pre-scans, and have opted to simply start the Markov-chain sampler off in the vicinity of the "correct" answer 
It is possible that as a result, the Markov-chain sampler may have failed to find other islands of high probability 
separated from the correct area by a deep ocean of improbability. Experimentation with less well constrained 
situations (e.g. only 20 events instead of 100) suggests that the other islands are not commonly a problem, but may 
be able to admit extra solutions where two particles in the chain are close to being mass degenerate. 

Remember that there is an uncertainty associated with the interpretation of this and the other reconstructed 
widths for the reasons discussed in Section l4!2l 

'^although there is no reason why a Bayesian analysis could not include extra hypotheses for additional simul- 
taneous chains 
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Part XI 

Study of non-pointing photons at the 
CERN LHC 

K. Kawagoe, M.M. Nojiri, G. Polesello and D. Prieur 
Abstract 

Measurement of non-pointing photons is a key issue to study the gauge 
mediation models at the CERN LHC. In this article we study the 
resolution of non-pointing photons with the ATLAS electromagnetic 
calorimeter, and discuss the impacts to the study of the gauge media- 
tion models. 

1. GAUGE MEDIATION MODELS AND NON-POINTING PHOTONS 

Origin of the SUSY breaking in the hidden sector and its mediation to the MSSM sector are 
key features of SUSY models. When hidden sector SUSY breaking is expressed by the order 
parameter F and the scale of the mediation to the MSSM sector by A/, the mass scale of MSSM 
sparticles Msusy is of the order of XF/M, where A is the coupling of the hidden sector to the 
MSSM sector. The SUSY breaking mediation may be due to renormalizable interactions, such 
as the gauge interaction. This is called "gauge mediation" (GM) models. In the GM models M 
and F are arbitrary and we expect M ^ Mp\. 

When M < Mpi, the lightest SUSY particle (LSP) is the gravitino (G) in the GM models. 
The next lightest SUSY particle (NLSP) is a particle in the MSSM sector which decays into a 
gravitino. If the lightest neutralino (xi) is the NLSP, the dominant decay mode is Xi ~^ iG. 
The neutralino lifetime cr is a function of Fq and rn^a, and the neutralino may be long-lived. 
Therefore it is an important subject to study non-pointing photons at the CERN LHC. 

In a paper [106], a procedure is proposed to solve the gravitino momentum and x^ decay 
position for the cascade decay i x^i — > Gi'y using the ATLAS detector at the LHC. It 
is shown that one can determine the mass and lifetime of Xi- To this purpose, one need to 
measure the photon momentum and arrival time very precisely. A toy simulation is made under 
the following assumptions to the photon momentum resolution; 

• A good angular resolution of ae = 60 mrad/A/E is assumed, where 9 is the polar angle of 
the photon momentum with respect to the beam axis and E is the photon energy measured 
in GeV. This resolution is based on a simulation for pointing photons. 

• The azimuthal angle of the photon momentum is only poorly measured by the electro- 
magnetic (EM) calorimeter. The resolution is good only for the photons converted in 
the transition radiation tracker (TRT) located in the EM calorimeter. 

The analysis takes two steps. First, events with converted non-pointing photons are se- 
lected. As the photon momentum is precisely measured, the events can be used to determine 
the mass of Xi- In solving the kinematics, the direction of the gravitino momentum is sensi- 
tive to the photon momentum. The mass resolution at GM point Gl with cr^o = 100cm is 
Am^o = 3.5 GeV for 10^ generated SUSY events, which corresponds to an integrated luminos- 
ity of 13.9 fb^^ 
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Once the mass is determined precisely, the good (p resolution of the photon momentum 
is not required to solve the decay kinematics. The decay kinematics can be solved for all £7 
events with non-pointing photons, from the 6 component of the photon momentum, the arrival 
time and position at the ECAL, and the lepton momentum. Therefore, all £7 events with or 
without photon conversion can be used to determine the lifetime. The lifetime resolution for the 
GM point Gl is /\ct/ct = 0.045 for the 10^ generated SUSY events. 

Although the result looks nice, some of the assumptions in the paper may be too opti- 
mistic, and should be studied more realistically. Among them, the 6 resolution obtained for 
pointing photons are used to estimate that for non-pointing photons. In this article we study the 
9 resolution for non-pointing photons by a full simulation and discuss the impacts to the GM 
study. 

2. PARAMETRISATION OF ANGULAR RESOLUTION OF ATLAS BARREL ELEC- 
TROMAGNETIC CALORIMETER 

In this part we try to obtain a more refined parametrization of the angular resolution for non- 
pointing photons with the ATLAS EM calorimeter. So far the resolution has been studied only 
for nearly pointing photons, i.e. photons coming from the ATLAS interaction point, and it is 
of the order of 60 mrad/v^ [79, 107]. In the following we will try to see which resolution is 
achievable for non-pointing photons using a full simulation of the ATLAS detector. 

The EM calorimeter [79] is a projective calorimeter with a good granularity to perform 
precision measurements of the shower position. It is longitudinally divided into three compart- 
ments: strip, middle and back compartment. The strip section is segmented along rj into very 
thin cells of At] = 0.003125, leading to a resolution on rj position with pointing photons of 
0.30 X 10"^. The middle compartment has a wider 77 granularity of Ar] = 0.025 and it is de- 
signed to contain most of the shower energy. It has a resolution on 77 position of 0.83 x 10~'^. 
By combining the measurement of the r] position in the first two compartments, it is possible to 
determine the shower direction in 77. 

For this study, different samples of single photons have been generated. Each of these 
samples consist of 20000 photons of pt = 60 GeV, randomly triggered with 77 from —1.4 to 1.4. 
Pointing photons were generated at ATLAS origin with a spread on the position of the genera- 
tion vertex of 5.6 cm along Z axis and 15 fim along radial axis, as it should be in ATLAS final 
setup. For non-pointing photon samples, the generation vertex has been shifted along ATLAS 
Z axis with values from 10 cm to 150 cm. No spread on the generation vertex position has been 
applied for these dataset. Finally all of these events were fully simulated using Dice/Atlsim 
(v3.2.1), the GeantS ATLAS detector description [108]. 

Here we focus only on the barrel part of the EM calorimeter. The reconstruction of all events 
has been done using ATLAS standard reconstruction software (Athena v6.5.0). No electronic 
noise or pile-up have been added in the reconstruction. The noise will certainly contribute to 
degrade the resolution, however this has not been studied yet. 

With the ?7 position, in each layer of the calorimeter (771,772) and using a parametrization of the 
shower depth for each layers (Ri{'r]i),R2iV2)) [107], we are able to reconstruct the shower axis 
Vpointing uslug the foUowlug relation: 
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The standard reconstruction uses a sliding window algorithm to find regions of interest in 
the calorimeter, then 3 x 3 (77 x 0) clusters are made in order to compute the t] barycenter of 
the shower in each layer. This position is the average of cells rj position weighted by energy in 
each cell: 

The angular resolution achieved using the standard reconstruction is shown in Fig. [l]as a func- 
tion of the position of the generation along ATLAS Z axis. The resolution is significantly 
degraded from 60 mrad for Z^^rtex = cm to about 800 mrad for Z^ertex = 100 cm. This is due 
to several effects. First, the S-shape corrections, that were tuned for pointing photons, are no 
longer valid and tend to degrade the resolution for large deviation from pointing. The S-shape 
effect is a distortion of the reconstructed 77 position due to the finite cluster size. Then the 3x3 
clusters are no longer sufficient to contain all the shower and some energy leakage outside the 
cluster is possible. Another point is that Eq. (O is no more true for non-pointing photons and 
gives rise to a systematic shift in computing 77 position for each layer. Finally, the shower depth 
parametrization, tuned from pointing photons, is also no more valid for large deviation from 
pointing. 

In order to improve the resolution, some changes have been made to the standard recon- 
struction algorithm. First, the cluster size has been extended to 5 x 3 and we do not apply 
the S-shape corrections. For each layer the systematic shift observed in the 77 position recon- 
struction has been parameterized as a function of the generation vertex position. Using this 
last parametrization we have made an iterative algorithm which corrects the rj position in each 
layer. The convergence is obtained in about 3 iterations. Results of this correction on angular 
resolution is show on Fig. [B For small vertex shift (Z < 30 cm), the standard reconstruction 
gives the best resolution. This is mainly due to the absence of S-shape corrections. For larger 
vertex shift {Z > 30 cm), the reconstruction algorithm with iterative correction gives better 
results than the standard one. 

This study has shown that standard reconstruction algorithm is not well suited for non- 
pointing photons and that a specific treatment is necessary. The next step to improve the resolu- 
tion would be to try other clustering algorithms such as NearestNeighbor. Another possible way 
would be to study the dissymmetric shower profile of non-pointing photons and try to extract 
an alternative method for computing 77 barycenter from this information. 



3. EFFECT ON THE GAUGE MEDIATION STUDY 

First of all, the 9 resolution affects the selection efficiency of non-pointing photons and back- 
ground contamination from pointing photons. Second, the resolution of the decay kinematics 
would be reduced if the photon angle resolution is reduced. We study this in terms of the angle 
tjj defined as the opening angle of the gravitino direction and the photon direction. This can be 
calculated by using the following formula 



cosip 



e 



ct^-L cos a 

where 4 = . , (3) 

L sm a 

where is the photon arrival time at the barrel EM calorimeter, L is the distance between the 
interaction point (O) and the point where the photon arrives at the EM calorimeter (A) , and 
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Figure 1 : Angular resolution of ATLAS barrel EM calorimeter as a function of the position of the production 
vertex. 
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cr6i(mrad) x y^E{GeV) 


60 100 200 300 400 500 


cr^(mrad) 


36.5 38.2 40.6 46.9 52.6 56.6 60.2 


RMS^, (mrad) 


46.2 48.0 50.9 63.0 78.8 93.0 105.5 



Table 1 : Resolution of the angle ip for various asumptions of 9 resolution. 

a is the angle between the photon momentum and the position vector oA. Changing the 9 
resolution and keeping all the other resolutions same as those assumed in Ref. [106], we show 
in Table[l]the resolution of the angle 'tp for point Gl with cr = 100 cm and / Cdt = 13.9 fb"^ 
Here, cr^ is obtained by Gaussian fit using center part of A^(= ip — iptYue) distribution, while 
RMS^, is obtained using the whole distribution. The error increases as ag increases. This 
affects the error of the Xi mass determination. 
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Part XII 

Measuring Neutrino Mixing Angles At 
LHC 

W. Porod and P. Skands 
Abstract 

We study an MSSM model with bilinear R-parity violation which is 
capable of explaining neutrino data while leading to testable predictions 
for ratios of LSP decay rates. Further, we estimate the precision with 
which such measurements could be carried out at the LHC. 

1. INTRODUCTION 

Recent neutrino experiments [109-112] clearly show that neutrinos are massive particles and 
that they mix. In supersymmetric models these findings can be explained by the usual seesaw 
mechanism [113-115]. However, supersymmetry allows for an alternative which is intrinsically 
supersymmetric, namely the breaking of R-parity. The simplest way to realize this idea is to 
add bilinear terms to the superpotential W: 

W = Wmssm + eiLiHu (1) 

For consistency one has also to add the corresponding bilinear terms to soft SUSY breaking 
which induce small vacuum expectation values (vevs) for the sneutrinos. These vevs in turn 
induce a mixing between neutrinos and neutralinos, giving mass to one neutrino at tree level. 
The second neutrino mass is induced by loop effects (see [116-118] and references therein). 
The same parameters that induce neutrino masses and mixings are also responsible for the decay 
of the lightest supersymmetric particle (LSP). This implies that there are correlations between 
neutrino physics and LSP decays [1 19-122]. 

In this note we investigate how well LHC can measure ratios of LSP branching ratios 
that are correlated to neutrino mixing angles in a scenario where the lightest neutralino x? is the 
LSP. In particular we focus on the semi-leptonic final states kq'q (k = e, jj, r). There are several 
more examples which are discussed in [120]. In the model specified by Eq. O the atmospheric 
mixing angle at tree level is given by 

tan6'atm = 7^ (2) 
Ai = eiVd + fiVi (3) 

where Vi are the sneutrino vevs and Vd is the vev of H^. It turns out that the dominant part of 
the Xi-W-/i coupling Of is given by 

Of = Ai/(Mi, M2, /i, tan/5, Vd, v^) (4) 

where the exact form of / can be found in Eq. (20) of ref. [120]. The important point is that 
/ only depends on MSSM parameters but not on the R-parity violating parameters. Putting 
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everything together one finds: 



atm 



A. 



A, 



BR{x\ r±WT) BR{x\ r±qq' 



(5) 



where the last equality is only approximate due to possible (small) contributions from three 
body decays of intermediate sleptons and squarks. The restriction to the hadronic final states of 
the W is necessary for the identification of the lepton flavour. Note that Eq. (|5ll is a prediction 
of the bilinear model independent of the R-parity conserving parameters. 



2. NUMERICAL RESULTS 

We take the SPSla mSUGRA benchmark point [30] as a specific example, characterized by 
mo = lOOGeV, mi = 250GeV, Aq = -lOOGeV, tan/3 = 10, and sign(/i) = V\ The low- 
energy parameters were derived using SPheno 2.2 [5] and passed to Pythia 6.3 [16] using 
the recently defined SUSY Les Houches Accord [19]. The R-parity violating parameters (in 
MeV) at the low scale are given by: ei = 43, 62 = 100, 63 = 10, Vi = —2.9, V2 = —6.7 
and V3 = —0.5. For the neutrino sector we find Aml^^^ = 3.8 ■ 10~^ eV^, tan^ 6'atm = 0.91, 
Am^Qj = 2.9 ■ 10^^ eV^, tan^ ^soi = 0.31. Moreover, we find that the following neutralino 
branching ratios are larger than 1%: 



BR(W±/i^) = 2.2%, BR(W±rT) = 3.2%, BR(qq>T) = 1.5%, 
BR(qqVT) = 2.1%, BR(qqi^i) = 4.7%, BR(bbz/i) = 15.6%, 

BR(e±rTz/i) = 5.9%, BR(/i±r=Fz/i) = 30.3%, BR(r+r-i^i) = 37.3%, 

where we have summed over the neutrino final states as well as over the first two generations of 
quarks. Moreover, there are 0.2% of neutralinos decaying invisibly into three neutrinos. In the 
case that such events can be identified they can be used to distinguish this model from a model 
with trilinear R-parity violating couplings because in the latter case they are absent. 

We now turn to the question to what extent the ratio, Eq. could be measurable at an 
LHC experiment. The intention here is merely to illustrate the phenomenology and to give a 
rough idea of the possibilities. For simplicity, we employ a number of shortcuts; e.g. detector 
energy resolution effects are ignored and events are only generated at the parton level. Thus, 
we label a final-state quark or gluon which has p± > 15GeV and which lies within the fiducial 
volume of the calorimeter, |?7| < 4.9, simply as 'a jet'. Charged leptons are required to lie 
within the inner detector coverage, |?7| < 2.5, and to have p± > 5GeV (electrons), p± > 6GeV 
(muons), or p± > 20GeV (taus). The assumed efficiencies for such leptons are [79] 75% for 
electrons, 95% for muons, and 85% for taus decaying in the 3-prong modes (we do not use the 
1-prong decays), independent of p±. 

For SPSla, the total SUSY cross section is ususy ~ 41 pb. This consists mainly of gluino 
and squark pair production followed by subsequent cascades down to the LSP, the With an 
integrated luminosity of 100 fb~^, approximately 8 million x? decays should thus have occurred 
in the detector. 

An important feature of the scenario considered here is that the Xi width is sufficiently 
small to result in a potentially observable displaced vertex. By comparing the decay length, 

'^Strictly speaking, the SPS points should be defined by their low-energy parameters as calculated with ISAJET 
7.58. 
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mode 






^gen 


f rcc 


A^,ec(100 fb"') 




I2W - 
tW - 


'Hj— prong 11 


235000 
51600 


0.10 
0.054 


12500 
1400 



Table 1: Statistical sample, estimated reconstruction efficiencies, and expected event numbers. 



CT = 0.5 mm, with an estimated vertex resolution of about 20 microns in the transverse plane 
and 0.5 mm along the beam axis, it is apparent that the two neutralino decay vertices should 
exhibit observable displacements in a fair fraction of events. Specifically, we require that both 
neutralino decays should occur outside an ellipsoid defined by 5 times the resolution. For at 
least one of the vertices (the 'signal' vertex), all three decay products (yuqq' or rqq') must be 
reconstructed, while we only require one reconstructed decay product (jet in the inner detector 
or lepton in the inner detector whose track does not intersect the 5a vertex resolution ellipsoid) 
for the second vertex (the 'tag' vertex). 

Naturally, since the decay occurs within the detector, the standard SUSY missing E± 
triggers are ineffective. Avoiding a discussion of detailed trigger menus (cf. [123]), we have 
approached the issue by requiring that each event contains either four jets, each with p± > 
lOOGeV, or two jets withp± > lOOGeV together with a lepton (here meaning muon or electron) 

withpx > 20GeV, or one jet with > 1 00 GeV together with two lep tons with > 20GeV. 
Further, since the Standard Model background will presumably be dominated by tt events, we 
impose an additional parton-level b jet veto. 

To estimate the efficiency with which decays into each channel can be reconstructed, a 
sample of 7.9 million SUSY events were generated with Pythia, and the above trigger and 
reconstruction cuts were imposed. To be conservative, we only include the resonant decay 
channels, where the quark pair at the signal vertex has the invariant mass of the W. The number 
of generated decays into each channel, the fractions remaining after cuts, and the expected total 
number of reconstructed events scaled to an integrated luminosity of 100 fb~^ are given in ta- 
ble [T] The comparatively small efficiencies owe mainly to the requirement that both neutralino 
decays should pass the 5a vertex resolution cut. Nonetheless, using these numbers as a first esti- 
mate, the expected statistical accuracy of the ratio, i? = BR{xi f^'^W^) / B R^xl t'^W'^), 
appearmg in Eq. © becomes ^ ~ 0.028. 

3. CONCLUSIONS 

We have studied neutralino decays in a model where bilinear R-parity violating terms are added 
to the usual MSSM Lagrangian. This model can successfully explain neutrino data and leads 
at the same time to predictions for ratios of the LSP decay branching ratios. In particular 
we have considered a scenario where the lightest neutralino is the LSP. In this case the ratio 
BR{xi f^'^W^) /BR(xi t'^W^) is directly related to the atmospheric neutrino mixing 
angle. Provided R-parity violating SUSY is discovered, the measurement of this ratio at col- 
liders would thus constitute an important test of the hypothesis of a super symmetric origin of 
neutrino masses. 

We have investigated the possibility of performing this measurement at a 'generic' LHC 
experiment, using PYTHIA to generate LHC SUSY events at the parton level and imposing 
semi-realistic acceptance and reconstruction cuts. Within this simplified framework, we find 
that the LHC should be sensitive to a possible connection between R-parity violating LSP de- 
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cays and the atmospheric mixing angle, at least for scenarios with a fairly light sparticle spec- 
trum and where the neutralino decay length is sufficiently large to give observable displaced 
vertices. Obviously, the numbers presented here represent crude estimates and should not be 
taken too literally. A more refined experimental analysis would be necessary for more definitive 
conclusions to be drawn. 
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Part XIII 

Resonant slepton production at the LHC 
in models with an ultraUght gravitino 

B. C. Allanach, M. Guchait, K. Sridhar 
Abstract 

We examine resonant slepton production at the LHC with gravitinos in 
the final state. We investigate two cases: (i) where the slepton undergoes 
gauge decay into neutralino and a lepton, followed by the neutralino 
decay into a photon and a gravitino, and (ii) direct decays of a slepton 
into a lepton and a gravitino. We show how to accurately reconstruct 
both the slepton and neutralino masses in the first case, and the slepton 
mass in the second case for 300 fb^ of integrated luminosity at the 
LHC. 

1. INTRODUCTION 

This letter is devoted to the study of the signals at the Large Hadron Collider (LHC) due to a 
supersymmetric generalisation of the Standard Model (SM) which (a) violates i?-parity, and (b) 
has an ultra-light gravitino in its spectrum. The anomalous events in the CDF experiment in the 
production rate of lepton-photon-^^j. in pp collisions were explained [124-126] in the frame- 
work of a i?-parity violating supersymmetric model with dominant L-violating coupling, 
and an ultra-light gravitino of mass ~ 10"^ eV. 

The resonant production of a smuon via the i?- violating coupling, its decay into neutralino 
and a muon and, finally, the decay of the neutralino into a gravitino and a photon leads to the 
/i7^y final state studied in the CDF experiment. The range of smuon and neutralino masses 
reljevant to the explanation of these anomalous observations of the the CDF experiment is 
such that most of this range will be explored at the Run II of the Tevatron. In the event that 
this signal is not seen at Run II it will rule out the model at the lower end of the neutralino 
and smuon masses. For heavier smuon and neutralino masses (above 250 GeV, roughly), the 
aforementioned Run I signal would be a statistical fluke and will probably disappear in Run 
II data. In that case, experiments at the LHC can be expected to discover and measure the 
sparticles. Here, we perform a study of the ability of the LHC to perform these two tasks, 
identifying the sensitive observables. 

2. THE MODEL 

We assume a single dominant R violating coupling, for example. If the i?-violating cou- 
pling is small, the existence of an ultralight gravitino in the mass range of lO^'^ eV drastically 
alters the decay mode of the slepton. The slepton overwhelmingly decays into a lepton and a 
(bino-dominated) neutralino, with the latter decaying into a photon and a gravitino resulting in 
a I'-f^T final-state. The Feynman diagram for the process is shown in Fig.[T] We should also ex- 
pect signals from sneutrino production. The background to the 7-^^ final-state that this would 
give rise to depend crucially on cosmic ray events which are difficult to estimate. Therefore 
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we have not studied the signal from sneutrino production. All sparticles except the neutralino, 
gravitino and slepton are set to be arbitrarily heavy in our analysis. 

3. SIMULATION RESULTS 

For our study of the process shown in Fig.[l]at the LHC {pp collisions at ^ = 14 TeV), we have 
chosen to work with the following default set of model parameters (unless indicated otherwise): 

• Gravitino mass, = 10"^ eV, 

• i?- violating coupling A' = X'211 = 0.01, 

• tan/? = 10, 

• sparticle masses (m^o , m;-)=(120 GeV,200 GeV) or (200 GeV,500 GeV) GeV ("low mass" 
and "high mass" scenarios) respectively. 

The choice of using X'211 rather than some other flavour combination is arbitrary and can be 
easily generalised to other i?,- violating couplings. We have checked that the chosen value for 
A211 is quite consistent with the existing bound [127, 128]. By selecting rather low values for 
-R-parity coupling and gravitino mass, we avoid significant rates for the possible _R-violating de- 
cays of Xi l^jj or ujj. Xi iG is the dominant channel. We use the ISASUSY [129] 
to generate the SUSY spectrum, branching ratios and decays of the sparticles selecting a repre- 
sentative point tan /3=10, Af,r,b = along with large values of ji and other flavour diagonal soft 
supersymmetry breaking parameters. 

The signals have been simulated using HERWIG6 . 4 [130] and the W'-) SM background 
has been simulated using pythia [16]. In our simulations, both the signal and background, we 
have used only selection cuts of Et, J^t > 25 GeV on the transverse energies of the muon, the 
photon and missing energy. We have used the following cuts on the rapidity of the photon and 
the muon: \ri^^^i\ < 3. There is an isolation cut between the photon and other hard objects o in 
the event of v^(^/7^-"^/^)M-707^-"^oF > 0.7. Since the signal is hadronically quiet, we veto 
events with jets reconstructed with Et > 30 GeV and rjj < 4. Initial and final state radiation 
effects, as well as fragmentation effects are included in the background simulation. 

The transverse mass distributions of final state particles along with show a clear dis- 
tinction between signal and backgrounds. The MTifi-y^T) Mt{'J^t) distributions are 
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Figure 2: Mt distributions of (a) nj^j-, (b) for slepton. 300 fb^^ integrated luminosity at the LHC is 
assumed. The purple (lighter) histograms display W'^ SM background, the red (darker) histograms show signal 
plus background for {m^o , TOj)=(120 GeV,200 GeV), whereas the blue (dotted) histograms display the signal plus 
background distributions for {m^o , m;-)=(200 GeV,500 GeV). In (a), the dashed black line displays a log-linear fit 
to the background distribution for Mt = 150 — 400 GeV. In (b), the insert shows a linear scale magnification of 
an area of the plot. 



displayed in Figs.|2h,b for the simulated high and low mass points and W'j simulated SM back- 
ground. In Fig.|2k, sharp peaks which are expected in the MrifJ^'j^T) distributions are clearly 
visible at values of the smuon mass and will be detected above the SM background, leading 
to the accurate measurement of the mas. Fig.|2jD shows that the signal peaks in Mt(7-^t) (prc- 
dicted to be at the neutralino mass) should be able to provide a measurement of the neutralino 
mass. 

In order to calculate the search reach, we use the signal S in the 4 highest peak bins (cov- 
ering 20 GeV) of the signal A<fT(/i7^y) peak. The background distribution in these four bins 
is estimated by fitting a simple function B = 4 exp [aMxifJ^'y^'r) + ^] ^ Mxi'jfi^'r) between 
150-400 GeV in Fig.EK- Using purely ^/B statistical errors, we obtain to a = —0.018 ± 0.001, 
b = 9.25 ± 0.22. B is displayed in Fig. |2^ as the dashed black line. We show the region of 
parameter space corresponding to S/ ^/B > 5 and 5 > 10 for 300 fb^^ luminosity option, as 
a function of smuon mass and R-parity conserving coupling in Fig.|3ji. 

We now turn to the decay I Gl. We ignore sneutrino production in this case because 
it would lead to an invisible final state. We have calculated the production matrix element 
and the branching ratio and implemented in a parton-level Monte Carlo. We have used cuts 
in our analysis on the muon and missing transverse energy identical to the 'jfi^r analysis, 
i.e. ^T^El^ > 25 GeV and \r]^\ < 3. Fig.EK displays the Mrifi^r) distribution for the W 
background plus signal in the cases = 1.2 and 1.5 TeV respectively and two different values 
of the A'. For A' = 0.01 there are not enough signal events to be seen, but for higher values (eg 
0.1), a clear mass peak should be seen in the tail of the Mt distribution of the W. We define the 
search reach by the criteria that for S signal and B background events, S/ VB > 5 and S* > 10 

''^The statistical uncertainties on fitted a and b parameters make a negligible difference to the final numerical 
results. 
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Figure 3: Search reach for the nj^x signal (as defined in the text) for 300 fb ^ integrated luminosity at the LHC. 




Mj{\l ET(miss)) (GeV) HI- (eV) 

Figure 4: (a) Mp distribution of the /iG final state for two different values of the smuon mass and 300 fb~^ 
integrated luminosity at the LHC. The left(right)-most three peaks are for smuon masses of 1.2 TeV and 1.5 TeV 
respectively, (b) Search reach (as defined in the text) for the fiG final state and 300 fb~^ integrated luminosity at 
the LHC and A211 = 0.1. The search reach is contained to the left of the curve. 
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in the 5 GeV signal peak bin of MxifJ^^T) ■ It is displayed for 300 fb ^ of integrated luminosity 
in Fig.|4j). 

CONCLUSIONS 

Resonant slepton production and its decays into I'jG or IG can be discovered at the LHC for 
slepton masses into the multi-TeV region, depending upon the Rp violating coupling and pro- 
vided that the gravitino is ultra-light (with a mass less than 0.1 eV). Various Mt distributions 
will allow the accurate measurement of sparticle masses involved. For full details, the new 
matrix element and additional results (for example including the case of a smuon NLSP), see 
ref. [131]. 
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Part XIV 

Radion Mixing Effects In The 
Two-Higgs-Doublet Model 

J.L. Hewett and T.G. Rizzo 
Abstract 

We begin an examination of the effects of mixing between the radion 
of the Randall- Sundrum (RS) model and the Higgs fields of the Two- 
Higgs-Doublet model as would be motivated by, e.g., supersymmetry. 
Preliminary results for the shifts in various particle masses and cou- 
plings are obtained. 

1. INTRODUCTION 

The RS model [132] provides an interesting solution to the hierarchy problem which can be 
tested experimentally [133-135] at future colliders. One prediction of this model is the exis- 
tence of a relatively light scalar radion which can mix with other scalars such as the Higgs boson 
of the Standard Model (SM). Such mixing can lead to substantial modifications in the expected 
properties of both the Higgs and the radion and has been extensively studied [136-139] in the 
literature. Here we extend this study to the case of two Higgs doublets as would be expected in 
a number of scenarios, e.g., supersymmetry (SUSY). Although not necessary for solving the hi- 
erarchy problem within the RS scenario, SUSY may have other model building [140, 141] uses 
such as coupling constant unification or radius stabilisation [142]. The expectation from previ- 
ous analyses of the single doublet model is that the properties of the mass eigenstate CP even 
neutral fields would substantially differ from those predicted in either the SM or the Minimal 
Super symmetric Standard Model (MSSM). The preliminary discussion of our findings given 
here supports these expectations though further study is required to understand the breath of the 
possible modifications. 

2. ANALYSIS 

With two Higgs doublets, mixing arises from the TeV brane action 

Smi. = [ d^x^/d^g R{g) {^MHi + 6i^li^2 + ^12HIH2 + h.c.) , (1) 

JTeV 

where g symbolises the induced metric on the TeV brane, R, the induced curvature arising from 
g. Hi are the two Higgs doublets and are dimensionless, order one parameters which we take 
to be real, thus assuming CP conservation for simplicity. The possible complexity of ,^12 may 
lead to interesting phenomenology. (In what follows ,^1 = ^2 will be assumed since it is unlikely 
that gravitational interactions distinguish between these two Higgs doublets. This assumption, 
however, may be incorrect.) Thus, in the unitary gauge, the CP-odd field. A, as well as the 
charged scalars are not directly affected by modifications to the neutral CP even sector. 
(Of course, their couplings to the CP-even fields will be modified.) Before mixing with the 
radion, tq, we denote the usual two CP-even Higgs by /iq, Hq which have been obtained from 
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the weak interaction eigenstate fields via rotations by the angles a,/3 as usual. Assuming tq 
obtains a mass from the stabilisation procedure, the above action can be expanded to quadratic 
order in the CP-even neutral fields from which one obtains the following effective Lagrangian 
generalising the notation of Csaki etal. [137]: 

1 1 

^ = -^^od^ho - -ml^hl + (/iq ^ i^o, '^o) - 3a7Vo9Vo + 6'j{Thhod'^ro + /iq ^ Hq) , (2) 

where 7 = v/^/GAt, ~ 0.05 — 0.10, with v the SM vev, and A^^, of order a few TeV or so, 
being the 'TeV scale' of the RS model. The parameters a, t^^h are functions of the and the 
usual mixing angles a, (3. The kinetic mixing in the above Lagrangian can be removed by a set 
of field redefinitions, i.e., ho, Hq h', H' + Q^yTh^ur' /Z and tq r' /Z. To obtain the mass 
eigenstate basis after the above redefinitions are employed a further orthogonal transformation, 
{h\ H'r') = 0{h, H, r), must be performed. The elements of the orthogonal matrix, O, as well 
as the corresponding mass eigenvalues can then be determined analytically. 




-10 10 -4-2024 



fl fl 

Figure 1 : Constraints from the absence of ghosts/tachyons(left) and perturbative unitarity(right) for different values 
of the ratio R in the range between -2 and +2. The allowed region lies between the curves. 

The number of parameters in this model is unfortunately rather large making it difficult to 
analyse in all generality. In order to show a specific example with somewhat fewer parameters 
we assume that the spectra, couplings, etc., of the unmixed two-doublet-model Higgs sector to 
be that given by SUSY, including the effects of radiative corrections [143], with tan/3 = 10, 
Ma = 500 GeV, At = = -yu=l TeV and Af| = (M^ + Af2)/2 = 1 TeV^ where the Mt^ 
are the stop masses. (The remaining parameters are ^1, m^o, 7 and the ratio R = ^12/^1-) These 
parameter choices yield nih ~ 125 GeV. The resulting parameter space can be further restricted 
by noting that both ,^1 and R are of order unity, 7 is anticipated to be near the range described 
above and is expected to be of order the weak scale. Further numerical restrictions on 
the parameter ranges can be obtained by demanding that, e.g., there are no ghosts or tachyons 
in the spectrum arising from the diagonalisation process and that W^W^ scattering satisfies 
perturbative unitarity [144] up to the scale A^; samples of such constraints can be seen in the 
figure above where we have assumed that —2<R<2. Note that for small values of 7 the 
constraints from unitarity are more restrictive than the requirement that no ghosts or tachyons 
be present; for larger 7 both constraints are found to be of comparable strength. 



76 



3. PRELIMINARY RESULTS 

We now turn to a brief sample of our results which survey only a Small region of the allowed 
parameter space. The first thing to examine is the effect of mixing on the masses of the physical 
states h, H, r as is shown in the next set of figures. (The curves are cut off at large values 
of iCil by the no ghost/tachyon requirement.) Here we see the typical result that the mass 
'levels' of the three states repel each other due to mixing. In all of our examples, uiho is the 
largest mass parameter and thus the mass of the H is itself raised by mixing. Note that these 
upward shifts can be enormous near the parameter space boundaries. Though made heavier, 
H's couplings will be seen to grow as well. When nih^ < {>)mra we see that the h(r) mass is 
pushed downward while that of r{h) is pushed upward. These shifts are not as large as those 
experienced by the H. Thus while we might expect the light Higgs to have a mass < 130 GeV 
in the MSSM, mixing with the radion allows it to be larger provided the radion itself is less 
massive. 

Next we consider the shifts due to mixing in the squares of the couplings of the various 
fields to either bb or ti(cc); it is important to note that the corresponding shifts for the couplings 
to WW/ZZ are found to almost identical to those for bb in almost all cases so we do not 
present those results separately here. Note that we have scaled these couplings shown in the 
figures either by their value in the MSSM using the input parameters above, as in the case of 
h, H, or by their unmixed values, as in the case of r. The range of coupling shifts, especially 
for H and r, are truly impressive being orders of magnitude in some cases. Not only are large 
enhancements seen for some parameter ranges, it is also important to note that there are regions 
of the parameter space, not necessarily near the boundaries, where couplings can completely 
vanish. For the light Higgs we see that while the mixing effects are not as large, for the cases 
at hand they always lead to a reduction in the coupling strengths in comparison with MSSM 
expectations. 

For the light Higgs, h, it is particularly important to examine the variations in the cou- 
plings to the gg and 77 final state since these control the dominant signal rate at the LHC; these 
are shown in the next set of figures. While the couplings of h to quarks and massive vector 
bosons was generally reduced via mixing, we see here that the loop-induced processes can be 
either enhanced or suppressed depending upon the region of parameter space we happen to be 
sitting in. The shift in the couplings of H, r to gg and 77 will be presented elsewhere. 

4. CONCLUSIONS 

We have begun a preliminary examination of the effects of mixing between the radion and 
the two CP even Higgs fields present in the two-doublet-model. As a result of this mixing 
the masses and couplings of all of these fields are found to be substantially modified from the 
expectations of the MSSM. Further analysis into the details of these mixing effects and the 
corresponding shifts in the various particle widths is on-going. 
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Figure 2: Typical shifts in light Higgs(top), heavy Higgs(middle) and radion(bottom) masses due to mixing for 
different values of the model parameters. Figures in a given column are case correlated. 
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Figure 4: Shifts in gg(top) and 77(bottom) couplings of the h due to mixing. Figures in a given column are 
correlated. 
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Part XV 

Search For The Radion Decay ^ hh 
With 77+bb, rr+bb And bb+bb Final 
States In CMS 

D. Dominici, G. Dewhirst, S. Gennai, L. Fanb, A. Nikitenko 

1. INTRODUCTION 

The Randall Sundrum model (RS) [132, 145] has recently received much attention because 
it could provide a solution to the hierarchy problem, by means of an exponential factor in a 
five dimensional non-factorizable metric. In the simplest version the RS model is based on 
a five dimensional universe with two four dimensional hypersurfaces (branes), located at the 
boundary of the fifth coordinate y. By placing all the Standard Model fields on the visible 
brane at y = 1/2 all the mass terms, which are of the order of the Planck mass, are rescaled 
by the exponential factor, to a scale of the order of a TeV. The fluctuations in the metric in 
the fifth dimension are described in terms of a scalar field, the radion which in general mixes 
with the Higgs. This scalar sector of the RS model is parametrized in terms of a dimensionless 
parameter ^, of the Higgs and radion masses mh, m^f, and the vacuum expectation value of the 
radion field A,^. The phenomenology of the Higgs and radion at LHC has been the object of 
several studies [136, 138, 139, 146-148] concentrating mainly on Higgs and radion production. 
In general the Higgs and radion detection is not guaranteed in all the parameter space region. 
The presence in the Higgs radion sector of trilinear terms opens up the important possibility 
of ^ hh decay and h (jxp. For example for mh = 120 GeV/c^, = 5 TeV and 
~ 250 - 350 GeV/c^ the BR{<p hh) ~ 0.2 - 0.3. In this paper we estimate the CMS 
discovery potential for the radion (0) in two Higgs decay mode (0 — hh) with 77+bb, rr+bb 
and bb+bb final states. 

2. ANALYSIS 

One point of m0=3OO GeV/c^ and mh=125 GeV/c^ was taken and the observability in the 
A^) plane was evaluated. Signal events were simulated with the PYTHIA [149] MSSM gg 
H — > hh process when values of mn and mh were set to 300 and 125 GeV/c^. Fh was set to ~ 
1 GeV/c^, thus the variation of the radion width in the (^, A,^) plane was neglected. However, 
the effect of changing of the radion width will be not visible due to the fact that F,^ is small in 
comparison with the detector resolution when the radion mass is reconstructed. 

2.1 77bb final state 

Signal events were processed with the full detector simulation and reconstruction. Events were 
required to pass the Level- 1 and High Level Trigger (HLT) selections for the di -photon stream 
[150, 151] with HLT thresholds on photons of 40 and 25 GeV/c. Photons were required to 
be isolated with the tracker and the electromagnetic calorimeter. The two highest Et jets of 
Et > 30 GeV and |?7| <2.4 were reconstructed with the calorimeter and were taken as b-jet 
candidates from h ^ bb decay. At least one of these two jets has to be tagged as a b jet. The 
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efficiency of single b tagging is 0.61 per event. Further selections require the di-jet mass, Mbj, 
to be in the window mh± 30 GeV/c^ (efficiency 65 %) and the di-photon mass, M^^, to be in the 
window mh± 2 GeV/c^ (efficiency 78 %). Finaly, the M^^bj niass should be in the window of 
m<^± 50 GeV/c^ (efficiency 95 %). Figure [T] shows Mbj and M^^bj distributions with arbitrary 
normalization after selections. The signal efficiency of the whole selection chain is 3.7 %. For 
= 1 TeV and ^ = the expected number of signal events with 30 fb"^ is 41 (with a(gg 0) 
= 40.8 pb, Br(0 ^ hh) = 0.33, Br(h ^ bb) = 0.61, Br(h ^ 77) = 0.00225). 




Figure 1: Reconstructed b j (left) and 77bj mass for the signal. 



Irreducible di-photon backgrounds were generated with CompHEP [14] for the 77]] (j=u,d,s,g) 
process and with MadGraph [152, 153] for the 77CC and 77bb processes with the factorization 
and the renormalization scales set to Mz and CTEQ5L PDF. The generator level preselections 
are PT,i„ax (min) >35 (20) GcV/c, p^^ >20GeV/c, \ri\ <2.5, AR^^ >0.3, AR^j >0.3, ARjj >0.3. 
Cross sections are shown in Table [T] PYTHIA was used for the hadronization. Initial and final 
state radiation in PYTHIA (ISR, FSR) were switched on. A fast detector simulation with the 
realistic resolution of the photon and jet energies and track momentum was used. The track 
reconstruction efficiency of 0.9 was taken into account in the tracker isolation criteria. The 
efficiency of b tagging, 0.5, for b jets and the mistagging probability of 0.01 (0.1) for u,d,s,g 
(c) jets were used. These numbers correspond to what was obtained with the full detector 
simulation [154] using the impact parameter tagging method. The efficiency of the selections 
and the expected number of the background events with 30 fb~^ after all selections including b 
tagging are shown in Table[T] Statistical errors on the expected number of events are also shown. 
The number of background events was then multiplied by 0.92 and by 0.90 to take into account 
of the Level- 1 e/7 trigger and the calorimeter isolation efficiencies which were not taken into 
account in the fast simulation. These efficiencies were obtained from a full simulation of the 
signal events. 

The CMS discovery reach was obtained in the A^) plane. Figure |21 (left plot) shows the 5 
a discovery contour in the the A^) plane when the irreducible background only (6.9 events 
with 30 fb"^) was taken into account. Theoretically excluded regions are also shown in the 
plot. Dashed line contours present the discovery reaches when the irreducible background cross 
sections were calculated for the renormalization and factorization scales set to 0.5 x /ig and to 
2x/io, where /iq = Mz. The background cross section uncertainty due to the scale variation 
found to be of ~ 40 % for 77bb. It was guessed that the cross section variation for 77CC and 
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Table 1: Background cross sections, efficiency, number of events with 30 fb ^ after all selections including b 
tagging^ 





77jj 


77CC 


77bb 


cross-section, fb 


13310 


778 


76 


selections 


efficiency 


E|''' > 40, 25 GeV, |r/| 


< 2.5 


0.446 


0.466 


0.487 


tracker isolation in cone 0.3 


0.328 


0.345 


0.379 


two jets Er > 30 GeV, 


\r]\ <2.4 


0.127 


0.125 


0.133 


M^^ window 4 GeV/c^ 


0.00278 


0.00263 


0.00410 


Mjj window 60 GeV/c'-^ 


0.00086 


0.00096 


0.00144 


M^^^jj window 100 GeV/c^ 


0.00045 


0.00061 


0.00123 


N events after all selections including b tagging 


4.2 ±0.8 


2.0 ± 0.6 


2.0 ± 0.6 



77jj production may be only slightly different, thus 40 % variation was applied to the cross 
sections of 77CC and 77]] processes as well. 
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Figure 2: 5 a discovery contours for ^ hh ^ 77-i-bb channel (m0=3OO GeV/c^, mh=125 GeV/c^); Left plot : 
with the irreducible background only; Right plot : with the total background assuming the ratio of the reducible to 
the total background of 0.4. Dashed line contours present the discovery reaches when the irreducible background 
cross sections were calculated for the renormalization and factorization scales set to 0.5x/io and to 2x/iQ where 
Mo = Mz. 



Reducible backgrounds from 7 + three jets and four-jet processes still have to be evaluated. It is 
expected from the inclusive h ^ 77 studies that the reducible background will be about of 40 % 
of the total background, thus the total background is expected to be of 1 1 .5 events with 30 fb^^ 
Figure |2l (right plot) shows the 5 a discovery contours in the the A^) plane with the total 
background taken into account. The experimental systematics uncertainty of the background 
estimated as ~ 5 % hardly affects the discovery reach due to the signal to background ratio in 
the 5 a region of the A^) plane is bigger than two. 
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2.2 rrbb final state 

The signature when one r lepton decays leptonically and another r lepton decays hadronically 
(producing a r jet) was considered. The highest signal cross section times the branching ratios 
of 0.96 pb was obtained for ^ = —0.35 and A,^ = 1 TeV. The background processes consid- 
ered in the analysis are shown in Table|21with the NLO cross sections taken from [155-157]. 
Background Z+jets (W+jets) were generated with px > 20 (80) GeV/c. Signal events were pro- 



Table 2: Trigger and total efficiency for the signal and the backgrounds; expected number of events with 30 fb 



samples 


a X BR (pb) 


efficiency (%) 


number of 
events 


trigger 


trigger -i- off-line 


(j) — * rrbb 


0.96 


6 ±0.2 


0.35 ±0.06 


102 


tt ^ 1 + z/ + jets + bb 


180 


0.57 ± 0.02 


(1.6 ±0.2)xlO-3 


111 


tt ^ 1 + z/ + r jet + bb 


15 


3.1 ±0.2 


(7.7 ±0.3)x 10-3 


66 


Zbb ^ rr + bb 


5.4 


1.4 ±0.2 


0.009±0.003 


21 


Z -|- jets ^ rr -|- jets 


306 


0.35 ± 0.02 


(3.3 ±0.5)x 10-4 


36 


W+jets 1 + 1/ + jets 


175 


0.039 ± 0.002 









cessed with the full detector simulation and reconstruction, while the background was processed 
with the fast detector simulation package CMSJET [158] . The combined electron(muon)-plus-r 
jet trigger [151] was used in this analysis. The Level- 1 trigger threshold is 21 GeV for electrons 
and 45 GeV for the r jet. The inclusive muon threshold is low enough (14 GeV/c) to allow 
a good efficiency however, to increase the background rejection a r jet with Et >35 GeV is 
required at the Level- 1 trigger. The trigger efficiency is shown in Table El for the signal and 
background samples. Missing transverse momentum and b tagging were used to reconstruct r 
leptons and to identify b jets coming from the two Higgs bosons. In order to increase the signal 
statistics it was necessary to tag at least one jet. The off-line selections are the following: 

- A(p between lepton and r jet direction >0. 1 ; 

- Et of the b-tagged jets >30 GeV, and Et of the most energetic jet > 55 GeV; 

- transverse mass of the lepton and missing momentum < 35 GeV/c^; 

- 75 < M,, <165 GeV/c2, 100 < Mbj <150 GeV/c^, 265 < M,,bj < 350 GeV/c^. 

Table |2l shows the signal and the background efficiencies for the off-line selections and the 
number of the expected events with 30 fb^^. Figure 13 shows the reconstructed rrbb mass after 
all selections (left plot). The total number of the background events after all selections is 234 
with 30 fb-^ The biggest background is tt (177 events) while the W ± jets background is 
negligible. Estimating also the contribution of tt when both W bosons decay into r the total 
number of background events increases up to 254. For the maximal signal cross section of 
0.96 pb, (^= —0.35 and = 1 TeV) about 102 signal events are expected. Signal significance 
(S/\/B) at this point is 6.4. Figure |3] (right plot) shows a 5 cr discovery contour in the A^) 
plane. The two contours correspond to the uncertainties of the background cross section values 
at NLO due to the scale variation and different PDFs [155-157]. The experimental systematics 
uncertainty of 3 % for the total background was taken into account. 
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Figure 3: Reconstructed rrbj mass for the signal and background (left) and 5 a discovery contour (right). The 
two contours correspond to the maximum and minimum k-factor value derives from the NLO order calculation of 
background cross section. 



2.3 bbbb final state 

The signal cross section times branching ratios for the gg ^ ^ hh ^ bbbb process is 
10.3 pb for = 1 TeV and ^ = — 0.35 The main QCD multi-jet background was generated by 
PYTHIA in different px bins. Other backgrounds considered are tt, ttjj and Zbb. In Table |3l 
cross sections and expected numbers of events with 30 fb^^ are summarized. A rejection factor 



Table 3: Signal and background events with 30fb . 





cross section 


events in 30 fb ^ 


Signal 


10.3 pb 


3.1 X 10^ 


QCDp^(30-l70) 


0.2257 mb 


6.79 X 10^2 


tt 


615 pb 


1.8 X 10^ 


ttjj 


507 pb 


1.5 X 10^ 


Zbb 


349 pb 


1.0 X 10^ 



on background higher than 10° is needed to reach a 5cr statistical evidence of the signal. A fast 
detector simulation with the CMS JET package [158] was used for both signal and background 
samples. Dedicated trigger selections were developed to keep the QCD multi-jet background 
rate at the acceptable level whilst maintaining a high efficiency for the signal. At Level- 1 
multi-jet triggers were used with the thresholds taken from Table 15-13 presented in [151] and 
restricted in pseudorapidity, |?7| <0.8. At the High Level Trigger at least 4 jets were required 
within the restricted pseudorapidity range, |?7| <0.8. Two jets must be b-tagged with the impact 
parameter tagging method (2 associated tracks with significance on the transverse impact pa- 
rameter > 2). The output QCD rate after these selections is ~ 5 Hz. In off-line selections all 
possible di-jet invariant masses were calculated from the 4 highest Et jets. The two jet pairs 
were chosen minimizing the value of mij-mk i, the same jets were then used to reconstruct the 
radion mass. The mean values (and cxfit) of the di-jet and four-jet effective masses reconstructed 
in this way are: 120 (39) GeV/c^ and 313 (76) GeV/c^. A 1.5 cr window in mass around mh and 
m<^ was used to select signal and background events. Efficiencies for the signal, background 
and the expected number of events with 30 fb"^ are summarized in Table IH The di-jet and 
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Table 4: Trigger and total selection efficiency and expected number of events with 30 fb 





^trigger 


^ total 


events 


signal 


0.038 


0.031 


9.57-10^ 


QCD pt(80 - 120) 


110-5 


7-10-6 


7.5-105 


QCD pt(120- 170) 


1-10-4 


6.6-10-5 


1.1-10^ 


tt 


0.015 


0.010 


1.84-105 


ttjj 


0.056 


0.026 


1.8-105 


Zbb ^ 4b 


0.002 


8-10-4 


1.2-10^ 



four-jet invariant mass for the background and the signal at A,^ = 1 TeV and C, = —0.35 point are 
shown respectively in the left and the right plots of Figured For this point it may be possible 
to achieve a signal significance (S/VB) of 5.5 if the background shape of the four-jet mass 
distribution is well understood (with ~ 0.1% uncertainty). 
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Figure 4: Reconstructed di-jet (left) and four-jet mass (right) for signal and background. 



3. CONCLUSION 

We estimated the CMS discovery potential for the radion into two Higgs decay mode (0 — > hh) 
with 77-i-bb, rr-i-bb and bb-i-bb final states. One point of m^=300 GeV/c^ and mh=125 GeV/c^ 
was taken and the observability in the A^) plane was evaluated. It was found that the 77-i-bb 
topology provides the best discovery potential. The bb-i-bb final state requires the dedicated 
High Level Trigger with the double b tagging and an excellent understanding of the background 
shape. 
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Part XVI 

The Invisible Higgs Decay Width In The 
ADD Model At The LHC 



Abstract 

Assuming flat universal extra dimensions, we demonstrate that for a 
light Higgs boson the process W*W*+X Higgs, graviscalars+ 
X — > invisible + X will be observable at the 5 a level at the LHC over 
the portion of the Higgs-graviscalar mixing and effective Planck 
mass (Md) parameter space where channels relying on visible Higgs 
decays fail to achieve a 5 u signal. Further, we show that for some val- 
ues of 6 and ^ ~ 0(1) the invisible decay signal can probe values of 
Md up to and possibly above those probed by the -independent) jets/7 
-I- missing energy signal from graviton radiation. We also discuss vari- 
ous effects, such as Higgs decay to two graviscalars, that could become 
important when nih/MD is of order 1. 



1. INTRODUCTION 

In several extensions of the Standard Model (SM) there exist mechanisms that modify the Higgs 
production/decay rates in channels that are observable at the LHC. One example is the Randall 
Sundrum model where the Higgs-radion mixing not only gives detectable reductions (or en- 
hancements) in Higgs yields, but also allows the possibility of direct observation of radion 
production and decay [139, 147]. It is also possible for the Higgs rate in visible channels to be 
reduced as a result of a substantial invisible width. For example, this occurs in supersymmetric 
models when the Higgs has a large branching ratio into the lightest gravitinos or neutralinos. 
Invisible decay of the Higgs is also predicted in models with large extra dimensions felt by 
gravity (ADD) [159, 160]. In ADD models the presence of an interaction between the Higgs 
H and the Ricci scalar curvature of the induced 4-dimensional metric gi-nd, generates, after the 
usual shift H = (^^, 0), the following mixing term [136] 



Above, Mp = (SttGjv)"^/^ is the Planck mass, 6 is the number of extra dimensions, ^ is a 
dimensionless parameter and is a graviscalar KK excitation with mass rni = 47r^n^/L^, L 
being the size of each of the extra dimensions. (Note that with respect to [136] our normalization 
is such that we have taken only the real part of the 0^. fields, writing </)g. = :^(sn + ^^n) and 

using (pQ = [(j^cPY to restrict sums to n > 0, by which we mean that the first non-zero entry of 
n is positive.) After diagonalization of the full mass-squared matrix one finds that the physical 
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eigenstate, h' , acquires admixtures of the graviscalar states and vice versa. Dropping C(e^) 
terms and higher, 



h' 



m>0 



irrihTh - m: 



2 '-'m 
fh 



mt 



imhTh - m: 



-h 



2 

m J 



(3) 



In computing a process such as WW ^ h' + J2m>o -^'m ~^ normalization and admixture 
corrections of order that are present must be taken into account and the full coherent sum 
over physical states must be performed. The result at the amplitude level is 



A{WW 



gWWhQhF 



mt 



(4) 



where -F(p^) = — e^Re X] 



^m>0 p2_^2_ 



and Gip^ 



-e^Im 



E 



m>0 



. Taking the am- 



plitude squared and integrating over dp in the narrow width approximation gives the result 



a{WW h'- 



" 7 > '^m 

m>0 



F) = asM{WW ^h^F) 



[m 



h ren/ 



where ml^^^ - ml 



+ F{m 



■T/i ~l~ T/i-^graviscalar 

(5) 



2 1 
hrcn) 



and we have defined m/iF/i^graviscaiar = G{m 



2 ^ 
hvcru 



We 



will argue that for a light Higgs boson both the wave function renormalization and the mass 
renormalization effects will be small. In this case, the coherently summed amplitude gives the 
Standard Model cross section suppressed by the ratio of the SM Higgs width to the sum of the 
SM Higgs width and the Higgs width arising from mixing with the graviscalars. 



2. INVISIBLE WIDTH 

As described, there is a decay of the Higgs arising from the mixing (or oscillation) of the 
Higgs itself into the closest KK graviscalar levels. These graviscalars are invisible since they 
are weakly interacting and mainly reside in the extra dimensions whereas the Higgs resides 
on the brane. The mixing width F/j^^graviscaiar ~ G{ml)/mh thus corresponds to an invisible 
decay width. The equation for G{m\) below eq. © shows that it is calculated by extracting the 
imaginary part of the mixing contribution to the Higgs self energy. The result is [136, 161] 

^ ^ 3(5 — 1) m^^^ 

T{h graviscalar) = T{h-^22 ■^") = 27r^^t;^ jjtfs^^^^ 

^ (:»«...^(^)-(^)%, 

where Ss-i = 27r'^/^/r(5/2) denotes the surface of a unit radius sphere in 6 dimensions while 
Md is related to the D dimensional reduced Planck constant by Md = (27r)'^/(^+'^)Mz). 
Our eqs. ^ are a factor of 2 larger than those presented in refs. [136, 161]. 



2.1 The wave function renormalization factor and mass renormalization 

A simple estimate of the quantity -^'("^Lcn)' appearing in the wave function renormalization 
factor found in eq. Q, suggests that it is of order where A is an unknown ultraviolet cutoff 
energy presumably of order A ~ Mo [162]. Assuming this to be the case, F' will provide a 
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correction to coherently computed LHC production cross sections that is very probably quite 
small for the rrih ^ Mb cases that we are about to explore. However, one must keep in 
mind that a precise calculation of F' is not possible. Similarly, the mass renormalization from 
-^("^Lcn) should be of order ^^m^/M^ and, therefore, small for mh -C M^. There are other 
incomputable sources of v^/M^ corrections lurking in the theory beyond these sources, and the 
results presented here are computed using the first, and perhaps only, calculable terms in the 
perturbation series. 



2.2 Contribution to the invisible width from direct two graviscalar decay 

In addition to decay by mixing, one expects also a contribution to the invisible width of the 
Higgs from its decays into two graviscalars. This can be evaluated by using the transformation 
of eq. ^ between the physical eigenstate h' and the unmixed h to derive the relevant trilin- 
ear h'skSi vertices. These are used to compute the corresponding matrix element. The final 
expression for T{h' — > graviscalar pairs) can be written as 



3+2(5„,2 / <- ^ N 2 r A/0 n 2 



rf/i' ^ graviscalar pairs) = — — ^ . , „ . f ^ 



5 





r Ti'/' 1 


(5+2) 


[n6/2)\ 



I, (7) 



where / is an integral coming from the sum over all the possible kinematically allowed h' — > 
SkSi decays. The integral / decreases rapidly as 6 increases. As a result, T(h' — > graviscalar pairs) 
is only significant compared with T(h graviscalar) if 5 < 4. The ratio of the two widths is 
given by: 

r(/i' graviscalar pairs) 3(6 — 1) 2/^ tt''/^ 



r(/i ^ graviscalar) 2TT'^{6 + 2y \Md J r(5/2)' 

From this result, we immediately see that even for small 5 the pair invisible width will be smaller 
than the mixing invisible width unless nih is comparable to M/j. 

To lowest order in ^'^{rrih/ M£,Y'^^ , decays of other states nearly degenerate with the h' 
can be neglected in the computation of a cross section obtained by coherently summing over 
the h' and the nearly degenerate s'^^ states. Thus, to this same order of approximation, V{h' 
graviscalar pairs) should simply be added to V{h — ^ graviscalar) in the expression for the 
narrow-width cross section of eq. Q. 

In Figure [H we show an extreme case corresponding to 5 = 2 and mh = 1000 GeV. 
Depending on the values of the parameters ^ and Mo, the pair invisible width can be a signif- 
icant correction to the invisible width from direct mixing. More generally, for rrih > Md the 
graviscalar-pair invisible width can provide a 3% to 20% correction to the direct-graviscalar- 
mixing invisible width. However, if rrih is substantially smaller than Md, then the graviscalar 
pair width is not important. For example, for 5 = 2, rrih = 120 GeV and = 500 GeV, 
T(h' — ^ graviscalar pairs)/r(/i — > graviscalar) < 0.0015 for ^ < 2. Therefore, in the follow- 
ing analysis, where we will assume a light Higgs, we can safely neglect the contribution to the 
invisible width from the decay into two graviscalars and use the expression given by eq. ©. 



3. MEASUREMENTS AT THE LHC 

For a Higgs boson with rrih below the WW threshold, the invisible width causes a significant 
suppression of the LHC Higgs rate in the standard visible channels. For example, for 5 = 2, 
Md = 500 GeV and rrih = 120 GeV, T{h graviscalar) is of order 50 GeV already by 
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Figure 1: In the left-hand plot, we display the total invisible width of a 1 TeV Higgs boson into one and two 
graviscalars as a function of Mjj for various values of ^ = 1 solid, ^ = 2 dashed, ^ = 3 dotted). For these 
plots we have fixed 5 = 2. The plot on the right shows the ratio of the two-graviscalars decay width to the 
one-graviscalar decay width for the same choices of parameters. 



,^ ~ 1, i.e. far larger than the SM prediction of 3.6 MeV. Even when nih is greater than the 
WW threshold, Figure[T]shows that the partial width into invisible states can be substantial even 
for M/p values of several TeV; therefore, for any given value of the Higgs boson mass, there is a 
considerable parameter space where the invisible decay width of the Higgs boson could be the 
first measured phenomenological effect from extra dimensions. 

Detailed studies of the Higgs boson signal significance, with inclusive production, have 
been carried out by the ATLAS [79] and CMS [163] experiments. If 115 GeV < nih < 
130 GeV, the /i — ^ 77 channel appears to be instrumental for obtaining a > Sa signal at low 
luminosity. The tih, h ^ bb and h — > ZZ* — > 4 £ channels also contribute, with lower statis- 
tics but a more favorable signal-to-background ratio. Preliminary results indicate that Higgs 
boson production in association with forward jets may also be considered as a discovery mode. 
However, here the background reduction strongly relies on the detailed detector response. 

In the ADD model, these results are modified by the appearance of an invisible de- 
cay width suppressing the Higgs signal in the standard visible channels. Here, we fix irih = 
120 GeV and perform a full scan of the ADD parameter space by varying M d and ^ for differ- 
ent values of the number of extra dimensions 6 and demonstrate that there are regions at high 
^ where the significance of the Higgs boson signal in the canonical channels drops below the 
5 cr threshold. However, the LHC experiments will also be sensitive to an invisibly decaying 
Higgs boson through l^H^-fusion production, with tagged forward jets. A detailed CMS study 
has shown that, with only 10 fb~^, an invisible channel rate of rinv/r=0. 12-0.20 times the SM 
WW ^Higgs production rate gives a signal exceeding the 95% CL significance for 120 GeV 
< nih < 400 GeV [163, 164]. Given that the effective (from the sum over the h state and 
nearby degenerate states) WWh coupling is of SM strength, this defines the region in the ADD 
parameter space where the Higgs boson signal can be recovered through its invisible decay . 

Figure |21 summarizes the results for specific choices of parameters. In the green (light 
grey) region, the Higgs signal in standard channels drops below the 5 a threshold with 30 fb"^ 
of LHC data. But in the area above the bold blue line the LHC search for invisible decays in the 
fusion channel yields a signal with an estimated significance exceeding 5 a. It is important to 
observe that, whenever the Higgs boson sensitivity is lost due to the suppression of the canonical 
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decay modes, the invisible rate is large enough to still ensure detection through a dedicated 
analysis. 



8 = 2 M„ = 120GeV 8 = 3 M„ = 120GeV 




Figure 2: Invisible decay width effects in the - Md plane for Mh = 120 GeV. The green (grey) regions indicate 
where the Higgs signal at the LHC drops below the 5 a threshold for 30 fb^^ of data. The regions above the blue 
(bold) line are the parts of the parameter space where the invisible Higgs signal in the 14^M^-fusion channel exceeds 
5 a significance. The vertical lines show the upper limit on Mo which can be probed by the analysis of jets/7 with 
missing energy at the LHC. The plots are for different values of 5: 2 (upper left), 3 (upper right), 4 (lower left), 5 
(lower right). 

The analysis of jet/7+ missing energy is also sensitive to the ADD model over a range of 
the Md and 5 parameters [165]. The invisible Higgs decay width appears to probe a parameter 
space up to, and possibly beyond, that accessible to these signatures (see Figure |2ll. Further, 
the sensitivity of these channels decreases significantly faster with 5 than that of the invisible 
Higgs width if ,^ ~ 1. Finally, it is interesting that, in the region where both signatures can be 
probed at the LHC, a combined analysis will provide a constraint on the fundamental theory 
parameters. A TeV-class e+e^ linear collider will be able to further improve the determination 
of the Higgs invisible width. Extracting the branching fraction into invisible final states from the 
Higgsstrahlung cross section and the sum of visible decay modes affords an accuracy of order 
0.2-0.03% for values of the invisible branching fraction in the range 0.1-0.5. But the ultimate 
accuracy can be obtained with a dedicated analysis looking for an invisible system recoiling 
against a Z boson in the e"'"e~ hZ process. A dedicated analysis has shown that an accuracy 
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0.025 < 5BR/BR < 0.04 can be obtained for 0.1 < BR < 0.5 [166]. This accuracy would 
establish an independent constraint on the Md, i and 5 parameters. 
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Part XVII 

Determining the extra-dimensional 
location of the Higgs boson 

A. Aranda, C. Baldzs, J. L. Diaz-Cruz S. Gascon- Shotkin and O. Ravat 
Abstract 

In the context of a TeV^^ size extra dimensional model, we consider 
the lightest Higgs boson as an admixture of brane and bulk scalar fields. 
We find that at the Tevatron Run 2 or at the LC the Higgs signal is 
suppressed. Meanwhile, at the LHC or at CLIC one might find highly 
enhanced production rates. This will enable the latter experiments to 
distinguish between the extra dimensional and the SM for up to 
about 6 TeV and perhaps determine the extra-dimensional location of 
the lightest Higgs boson. 

1. INTRODUCTION 

Extra dimensional models have been used recently to address a wide class of problems in par- 
ticle physics, such as the hierarchy, unification and flavor problems [159, 160, 167, 168]. In this 
work, we examine a TeV^^ size extra dimensional model, in which the lightest Higgs boson 
emerges as an admixture of brane and bulk scalar fields. This model predicts a suppression 
of the Higgs production cross section at LEP and the Tevatron, while it promises a significant 
enhancement of the signal at the CERN Large Hadron Collider (LHC) and possibly at a multi 
TeV linear collider (CLIC). We present results for the cross section of the associated production 
of Higgs with gauge bosons at the LC and LHC. 

2. THE MODEL 

In this section we present the general features of the extra dimensional model (for a detailed 
description see [169]). We work with a five dimensional (5D) extension of the SM that contains 
two Higgs doublets. The SM fermions and one Higgs doublet ($„) live on a 4D boundary, the 
brane, while the gauge bosons and the second Higgs doublet (^d), are all allowed to propagate in 
the bulk. The constraints from electroweak precision data [168] show that the compactification 
scale can be of C(TeV) (3-4 TeV at 95 % C.L.). The relevant terms of the 5D SU{2) x U{1) 
gauge and Higgs Lagrangian are given by 

= -\ {Flmf - \ {Bmn? + \DM^d? + \D,^^\'5{x'') , (1) 

where the Lorentz indices M and N run from to 4, and n runs from to 3. 

After spontaneous breaking of electroweak symmetry one obtains the following 4D La- 
grangian: 



^4 ^ ^^(^hsm{/3-a) + Hcos{/3-a))Z^Z^' 
2 c VI/ 
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Figure 1: SM, THDM and XD cross sections for e+e Zh. Each plot corresponds to a different set of values 
for a and (3 all with m^j = 120 GeV and with a compactification scale M^, = 4 TeV. 

+ V2^^(/isin/3cosa + i7sin/5sina) 

+ gMw {h sin(/3 - a) + H cos(/3 - a)) W+VT"^ 

oo 

+ V2gMw (/i sin/? cos a + sin /5 sin a) ^ (W+iy-^")^ + . (2) 

n=l 

where /i and if are the CP-even Higgses (m/i < mu), d is the mixing angle that appears in the 
diagonalization of the CP-even mass matrix, and tan (3 is the ratio of vevs. 



3. HIGGS PRODUCTION AT FUTURE COLLIDERS 

We present results for the associated h + Z production cross section obtained from Eq.© at 
linear colliders and at the LHC. Fig. [l] shows the results for the e+e" hZ cross section. asM 
stands for the standard model cross section, THDM labels the (4D) two Higgs doublet model, 
and the results from the extra dimensional model are denoted by XD. The three plots correspond 
to three different choices of the parameters a and (3. It can be observed that the SM cross section 
dominates in all cases up to ^ ~ 2 TeV. This is understood from the fact that the heavier KK 
modes, through their propagators, interfere destructively with the SM amplitude thus reducing 
the XD cross section. Moreover, as Fig. [l] shows, once the center of mass energy approaches 
the threshold for the production of the first KK state, the cross section starts growing. For 
instance, with Mc = 4 TeV, asM — <^xd for -y/s — 2 TeV. However, one would need higher 
energies in order to have a cross section larger than that of the SM, which may only be possible 
at CLIC [170]. Based on this, we also conclude that at the Tevatron the luminosity required to 
find a light Higgs boson is higher than in the SM case. 

The Higgs discovery potential in this model is more promising at the LHC. We illustrate 
this in Fig. El showing the pp hZ differential cross section as the function of the hZ invariant 
mass Mhz- The typical resonance structure displayed by Fig[T]is preserved by the hadronic 
cross section. The resonance peak is well pronounced when M^z ~ Mc. This leads to a 
large enhancement over the SM (or THDM) cross section. The singularity at Mc = M^z is 
regulated by the width of the KK mode, which is included in our calculation (cf. Ref. [169]). 
Depending on the particular values of a and /? the enhancement is more or less pronounced. 
For an optimistic set a = and P = 7r/2, the XD production cross section is considerably 
enhanced compared to the SM at Mhz = Mc. This enhancement may be detectable up to about 
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Figure 2: Higgs production cross section in association with a Z boson at the LHC as a function of the compacti- 
fication scale for selected values of the mixing parameters. 

Mc = 6 TeV. We estimate that with 100 fb~^ for Mc = 6 TeV there are about 20 hZ events in 
the bins around M^z ~ M^,. As Fig. |2l shows, in the SM less than one event is expected in the 
same M^z range. It is needless to say that similar results hold for hW^, which further 

enhances the discovery prospects. 

Based on these results, we conclude that in the Bjorken process alone the reach of the 
LHC may extend to about = 6 TeV, depending on the values of a and p. The Zh (or W^h) 
production cross section determines only a specific combination of the mixing angles a and 
p. In order to determine the individual angles, one has to also measure the production cross 
sections of the heavier Higgs boson in association with a Z (or W^). Fortunately, these cross 
sections are also enganced by the same amount as the ones for the lightest Higgs boson. 

4. CONCLUSIONS 

We presented Higgs production cross section calculations in the framework of an extra dimen- 
sional model that releaves the tension between the low mass predictions for the SM Higgs and 
the missing Higgs signal at LEP. We found a suppresson at the LC and at the Tevatron, but an 
enhanced signal at the LHC or CLIC. The fact that the lightest Higgs boson is an admixture 
of brane and bulk fields (that is it has a non-trivial bulk location) is the key ingredient in the 
suppression-enhancement mechanism for the signal. This may enable the LHC and CLIC to 
determine this location. 
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Part XVIII 

The sensitivity of the LHC for TeV scale 
dimensions in dijet production 

C. Baldzs, M. Escalier, S. Ferrag, B. Laforge and G. Polesello 
Abstract 

In this work, we present results for dijet distributions at the LHC with 
the assumption of a TeV size extra dimension. In our calculation, we 
included the virtual effects of gluonic Kaluza-Klein state exchanges, 
as well as the modified running of the strong coupling constant (but 
restricted our numerical study to the case of standard as evolution). 
Computing the transverse momentum distribution of dijets, we found 
that the LHC is able to discover a single extra dimension up to Mc ~ 15 
TeV. 

1. Introduction 

String theory is the most promising candidate for a unified framework of matter and interactions. 
Among the predictions of string theory are extra, compact space dimensions (XDs) which, 
depending on their sizes R, play a role in determining physics close to the weak scale. The 
string arguments of Refs. [171, 172] do not prevent the standard gauge bosons from penetrating 
the bulk. If the compactification scale = 1/R is higher than (9(TeV), phenomenology 
and present experiments do not conflict with this scenario either [167, 173-177]. These new 
dimensions, on the other hand, may be probed at near future particle accelerators, in particular 
at the CERN Large Hadron Collider (LHC). 

Since the LHC produces strongly interacting particles abundantly, if the gluons propagate 
in XDs then dijet production at the LHC is a sensitive discovery channel for TeV scale XDs. 
The main effect of these dimensions on the dijet production cross section is twofold. On the 
one hand, the gluonic Kaluza-Klein (KK) excitations enhance the dijet distributions in the high 
invariant mass (s) and transverse momentum (pr) region [178]. On the other, the modified 
evolution of the strong coupling (as) further distorts these distributions [179]. These competing 
effects are entangled and has to be taken into account simultaneously in order to predict the 
discovery potential of the LHC. 

In this work, we computed dijet distributions for the LHC including both of these ef- 
fects. We assumed that the standard model (SM) gauge bosons, especially gluons, propagate 
in a single TeV^^ size compact dimension. We implemented gluonic KK excitations and the 
modified running of as in the Monte Carlo event generator PYTHIA [15], for dijets. Then, we 
used PYTHIA to calculate dijet px distributions, determined the enhancement at large pr, and 
estimated the significance of a potential discovery. 

2. Dijet production at the LHC 

The formalism that we use is described in detail in Refs. [178] and [179]. When calculat- 
ing the dijet production cross section, besides the SM ones, we include the tree-level diagrams 
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shown in Fig. [T]'^ In these diagrams (7* signals that a KK tower of virtual gluons is exchanged. 
This means that in the SM diagrams we replace the gluon propagators by 



N 



(1) 



n=0 



Here 



Dnip) 



12 



(2) 



is the propagator of the n gluon KK resonance with p'. 
2asmn. The SM gluon is identified with the zero mode and cq 
When calculating the cross section, ^ ^ , 

it is necessary to evaluate amplitude squares 
which will contain products of propaga- 
tors of the form 
1 
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(We note that in Eq.©, we corrected a 
typo which is present in v.3 of hep-ph/0012259.)) 
As it was noted in Ref. [178] the sum in 
Eq. dSl) converges rapidly. We checked that 
for the LHC (with ^f(S) = 14 TeV) and 
for Mc > 1 TeV, choosing = 50 (or 
equivalently a cutoff scale of Ms > 50 
TeV) leads to a satisfactory numerical pre- 
cision. We implemented the effective prop- 
agators given by Eq.© in PYTHIA, mod- 
ifying the parton level processes represented 
by Fig.©. Finally, we checked the im- 
plementation against the numerical results Figure 1: Feynman diagrams for dijet production involving 
given in Ref. [178] and found a good agree- Kaluza-Klein excitations of the gluons. The indices i and j 
ment. represent distinct (i ^ j) quark flavors. 

In our calculation, we also include the modified running of as, as described in Refs. 
[167, 177, 179]. Above the compactification scale, we implemented the modified running of 
as given by 

5 




hi - hi 
2ti 



in 



rfj, 



dt 
T 



it 



(4) 



where i = 1,2,3 labels the gauge groups of the SM. The coefficients of the usual one loop beta 
functions 



{hi,h2,h-s) = (41/10,-19/6,-7) 



(5) 



'^We note that five-momentum conservation forbids internal gluonic KK excitations in any tree-level dijet dia- 
grams involving external gluons, that is the KK excitations do not affect the process qq gg, for example. 
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jets at LHC, 50 KK excitations 
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Figure 2: Left panel: number of dijet events vs. the dijet pr calculated in the SM and for various values of Mc- 
Right panel: KK contribution to the cross section as the function of pTmin- 



axe supplemented by new contributions from the KK towers 

(61, 62, h) = (3/5, -3, -6) + 7] (4, 4, 4) (6) 

(where, for simplicity, we set 77 = 0). In the last term of Eq. Q, denotes the elliptic Jacobi 
function and 

. = .TO-^" with = (7) 

We note that in Refs. [167, 177] an approximate expression is used to calculate the running of 
the couplings, but we implemented the exact formula ^ in PYTHIA. 

3. Numerical results 

Our goal is to quantify the sensitivity of the LHC to a TeV size XD. To this end, following 
Ref. [178], we compute the dijet transverse momentum (pr) distribution. As a first step, we 
include 50 KK excitations of the gluons but keep the standard evolution of a^. We use PYTHIA 
version 6.210 with the modifications described in the previous section. On the final state, we 
apply the following kinematic cuts: 

Pt > PTmin, \y\ < 2.5, pTjet > 100 GeV (oH each jet), (8) 

where y is the rapidity of the two jet system. In PYTHIA, we also turn the initial and final state 
radiation (ISR and FSR) on. 

The left panel of Fig. |2lshows the results as the number of dijet events against the dijet pt 
assuming 100 fb^ integrated luminosity. In this and in the subsequent computations we used 
the CTEQ6L1 parton distribution function (PDF) with the dijet invariant mass as the factoriza- 
tion scale. It is shown that there is a significant enhancement in the high for = 8 TeV, 
and there is still a detectable excess for = 15 TeV. The right panel of Fig. |2l shows the KK 
contribution to the cross section 



0"KK — CTtotal — CTSM, 



(9) 
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Figure 3: Statistical significance of the XD dijet signal in units of a. 

as the function of PTmin- Considering the difference in the PDF, slightly different cuts and the 
PYTHIA effects (ISR, FSR, etc.), our resuks still reasonably agree with Fig.3 of Ref. [178]. 

Finally, Fig. |3] shows the statistical significance 

S =\NsM-NxD\/^/N^, (10) 

of the XD dijet signal in units of a plotted against the dijet px- Here Nx is the number of events 
predicted by model X. This plots shows that by measuring the dijet px distribution in the 1-3 
TeV region the new dimension can be easily discovered even if it is as small as 15 TeV"^. 
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Part XIX 

Little Higgs Model: LHC Potential 

K.Mazumdar 

Abstract 

A recent idea of solving the 'little hierarchy' between physical mass 
of the Higgs boson and present limit of relatively low cut off scale of 
about 10 TeV has a vector- like new heavy quark Top as one of the addi- 
tional particles. The potential of CMS experiment at LHC to discover it 
has been studied using the available phenomenology for collider exper- 
iments. 

1. INTRODUCTION 

The Little Higgs model is an alternative solution for the hierarchy problem in the Standard 
Model (SM) of strong and electroweak interactions which has passed stringent experimental 
verifications upto the electroweak scale. The precision measurements indicate that the physical 
mass of the Higgs boson mn < 219 GeV/c^ to 95% CL. Considering the SM as an effective 
theory valid upto an energy scale ~ few TeV, any new physics appearing below a cutoff of 
about 10 TeV has to be weakly coupled to be consistent with precision measurements. The gap 
between electroweak scale and this cutoff value is called the 'little hierarchy'. 

This hierarchy problem is solved in Little Higgs model by requiring mn to be safe from 
only one-loop divergences. The Higgs fields are considered as Nambu-Goldstone-Bosons of a 
global symmetry which is spontaneously broken at a higher scale, but the physical Higgs boson 
states continue to remain light due to approximate global symmetry. In the minimal scenario 
both gauge and Yukawa interactions are necessary to break all the global symmetry which 
protect the Higgs boson mass and hence pushes the cutoff scale to ~ 10 TeV. The additional 
particles postulated in the model are three heavy gauge bosons, a heavy top-like quark (7), a 
heavy Higgs triplet and a Higgs doublet. The lightest neutral scalar is identified with the SM 
Higgs boson. Interestingly, the most important divergences are cancelled between loops of 
particles with the same statistics. 

The masses of the new particles need to be light enough to avoid fine tuning, whereas 
the precision electroweak constraints push them higher. Hence parameters of the model, eg. 
couplings, need to be tuned. Present expectations of 1-2 TeV mass range make these particles 
accessible at LHC and this motivates us to study their prospect of discovery. 

2. COLLIDER PHENOMENOLOGY 

The collider phenomenology has been worked out in [180]. The heavy gauge bosons can be 
produced in LHC via Drell-Yan type process. The heavy quark T can be produced via QCD 
(gluon-gluon fusion) in a model independent way but with increasing mass the cross-section 
falls off rapidly. It can also be produced singly, through t-channel fusion process W'^b T, 
and the cross-section, though larger, depends on the parameters in the model, relating the Top 
mass Mt to the scale, or equivalently to the Higgs condensate / (~ TeV). 
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The tree-level Feynman diagram for the production of single T and the cross-section at 
LHC as a function of its mass is given in Fig.[ll The single T production cross-section is taken 
as 0.25 pb for Mt = 1 TeV/c^. 

3. HEAVY TOP SEARCH IN CMS 

We have investigated the potential of CMS experiment at LHC to discover some of these new 
particles predicted in the model. We report in the following on detection capability and possible 
discovery of singly-produced T decaying through channels T ^ th and T tZ, each with a 
branching fraction of 25%, where t, h, Z are the top quark, the Higgs scalar and the neutral 
vector boson in SM. The other mode T hW , though have 50% branching ratio, is expected 
to be overwhelmed at LHC by backgrounds and hence we have not studied it. 

3.1 Simulation method 

We have used Pythia subprocess no. 83 for 4*^ generation heavy quark production in t-channel 
process q.-Jj Qkfi (Qk = T). Pythia decay table is also modified according to model 
predictions. Background events are various SM processes and they are generated using Pythia 
or Alpgen (multi-parton final states) generator. In all cases, CTEQ5L structure function is used. 

The detector response is emulated according to CMS experimental characteristics and 
using a fast monte carlo package CMSJET where detector effects are parametrised using full 
detector simulation. The jet reconstruction and identification is assumed to be possible upto 
pseudorapidity |?7| < 4.5 and above transverse energy Et > 40 GeV. The b-tagging of jets 
plays a crucial role and it is effective upto |?7| < 2.5 with an efficiency of about 60%. While 
considering possible bakground events, we note that the signal may be mimicked due to the 
mis-identification of a light quark jet as a b-jet, though the probability is about 1% only. 

We have first done a preliminary study of the kinematics of the signal event at particle 
level after event generation through Pythia as exemplified in Fig. |2| The source of the initial 
state b-quark is gluon splitting and being a t-channel process, both the initial state quarks are 
reasonably forward going. In the final state the particles are highly energetic as expected due 
to the heavy mass involved; the rapidity distributions of the decay products from T are central 
and the outgoing light quark, though not central, has |?7| < 4.5 for about 87% of the cases. 
Hence this light jet can be used for tagging the event selection as will be discussed later. Note, 
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Figure 1 : The tree-level Feynman diagram for the production of single T and the cross-section at LHC as a function 
of its mass. 
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Figure 2: Some of the kinematic distributions of the initial and final state particles in qb q'T. 

since these are preliminary studies our selections do not include any threshold on the missing 
transverse energy in the event, though it is reasonably high for signal events. 

4. HEAVY TOP SEARCH IN T ^ /it, hi MODE 

In this mode we have investigated the situation where the W from top-quark decays leptonically, 
i.e., T ^ ht, h ^ bb, t ^ bW, W ^ where £ = e, yU. For a higgs boson mass of 200 
GeV/c^, the decay branching to bb pair is about 65%, so the rate of the final state, given by 
the signal production cross-section times the relevant branching ratios, is 8.9 fb. The final 
state contains three b-jets, one isolated lepton, missing transverse energy and one light-jet. For 
background we have considered the following SM processes: 

• pp ^ ti production with subsequent decays of top-quarks to bW mode. We consider here 
the leptonic decay of one W and the hadronic decay of the other. We have used Pythia to 
generate such events and the leading order cross-section is 435 pb. 

• pp ^ Wbbjj with a total cross-section of 3.95 pb, including BR(W^ iu). 

• pp ^ W + Ajets with a total cross-section of 146 pb, including BR(11/ iv). 
For processes (2) and (3) we used Alpgen generator and are limited by statistics. 

4.1 Possibility of Trigger 

To judge the trigger viability of the signal events, we resort to CMS -designed set of criteria 
for triggering interesting physics events where the threshold values are decided by the DAQ 
capability. Triggering with an isolated lepton (inclusive) is a suitable condition in this case due 
to the leptonic decay mode of W. The reconstruction efficiencies for electron and muon events 
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are different and are taken into account to determine the signal selection efficiency according to 
the following set of kinematic criteria. 

• Electron: pr > 29 GeV/c, \r]\ < 2.5 with exclusion of 1.442 < \r]\ < 1.566 

• Muon: pr > 14 GeV/c, \r]\ < 2.1 

• Total selection efficiency for electron and muon events are 0.607 and 0.852 resp.ly. These 
numbers are also independent of mn in the interesting region. 

4.2 Event Selection 

We now discuss the main features of signal vis-a-vis the main background of ti events after 
they have passed the detector simulation and the trigger selection criteria. The momentum 
distribution of the lepton from the signal T ^ ht ^ 3 b-jets + 1 lepton + (light-jet) event is 
harder compared to that from ti bb, Wi lu, W2 jj event. The b-tagged jets are both 
central and also more energetic. The 4-th jet is typically due to a light quark and is mostly in 
the higher pseudo-rapidity region as expected for signal events. We have applied the following 
set of criteria to select signal events; the succesive efficiency factors are also quoted. 

• 1 isolated lepton at trigger level, efficiency = 0.73 . 

• Total no of jets with transverse energy Et > 60 GeV in the event = 4, among which three 
most energetic jets are b-tagged, efficiency = 0.007 . 

• First 3 jets should have minimum transverse energy > 200, 125, 80 Gev, |?7| < 2.3 
and 4th jet should be in the forward/backward region |?7| > 3.0, efficiency = 0.29 . 

• Transverse momentum of the isolated lepton, p^ > 40 GeV/c, efficiency = 0.82 . 

Combining all we get a grand efficiency of 0.12%, whereas the signal cross-section is 8.9 fb. 
Hence we conclude that the signal rate is too small to be observed even with large accumulated 
luminosity. 

The above selection criteria could be relaxed and optimised by demanding, for example, 
less number of b-tagged jets. Evidently with a demand of three b-tagged jets very few events 
from the background processes could survive due to jet mistagging, but none survived all the 
selections. This situation may change if we require only one or two b-tagged jets. We should 
also generate events of type ti +jets using Alepgen. 

To make our exercise complete, the mass of the higgs boson from the two b-jets with min- 
imal separation, AR = ^J\Aiff + A0^), is reconstructed and subsequently the transverse mass 
of T also can be reconstructed. This can finally lead to the estimation of the model parameter. 
We note here that the distance between the two b-quarks from the higgs boson decreases with 
increasing mn due to higher boost. 

5. HEAVY TOP SEARCH IN T ^ Zt, Z ^ l+t (i = e, //) MODE 

The complete final state considered is q b ^ q' T, T ^ Zt, Z e+£-,t ^ bW, W ^ iu 
and so there are 3 isolated charged leptons, one b-jet, one light jet and missing transverse energy. 
The cross-section is 0.92 fb. The main backgrounds are SM events of type WZ-hjets, with a 
cross-section of 91 fb, where both Z and W decay leptonically and the light-jet is mis-tagged as 
b-jet. As expected the leptons from the signal event have harder momentum distribution. 

5.1 Trigger 

The signal can be triggered in di-lepton mode with the following criteria. 
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• Electron: pr > 17GeV/c, \'q\ < 2.5 with exclusion of 1.442 < \ri\ < 1.566 

• Muon: pr > S GeV/c, \r]\ < 2.1 

• For two leptons of same flavour and opposite signs, to satisfy the criteria, the efficiency 
for electron type events is 0.77 and for muon type it is 0.92 . 

5.2 EVENT SELECTION 

The event selection criteria are discussed below with corresponding signal efficiencies for sub- 
sequent requirements. 

• 3 isolated leptons with any two passing di-lepton trigger criteria = 0.83 

• Maximum 2 jets in the event with Et > 40 GeV, among which one is b-tagged, efficiency 
= 0.049 

• Highest Et jet has transverse energy E^ > 100 GeV, efficiency =0.86 

• Invariant mass of 2 leptons, of same flavour and opposite sign Mu > 80 GeV/c^, effi- 
ciency = 0.98 

The WZ background, though large, is tamed by the requirement of accompanying hard 
b-jet and so no event finally survives the full set of requirements. The grand efficiency is 0.035 
and corresponding to a cross-section of 0.92 fb we expect about 3.2 events with an integrated 
luminosity of 100 fb^^. This gives an encouraging result of signal-to-background ratio of 6.4. 
By tuning the event selection criteria a better significance can be expected. 

6. CONCLUSION 

The Little Higgs model predicts the existence of a top-like heavy quark of mass within the range 
of few TeV We have studied the experimental signature of single production of this particle at 
LHC and subsequent decays in modes T ^ ht, h ^ bb and T Zt, Z i'^i' where 
W from top-quark decays leptonically. Preliminary studies show that for the first channel, a 
requirement of 3 b-tagged jets lead to too low an efficiency for the signal. But the second 
channel is encouraging and could be even more promising when we consider the hadronic mode 
of W since the event will be triggered by the leptons from Z anyway. The study of the channel 
where W 2-jets is under progress. Of course the background is expected to be larger. We 
plan to study the potential of CMS experiment for the heavy bosons soon. 
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Part XX 

TJ studies at the LHC: an update 

M. Dittmar, A. Djouadi, A. -S. Nicollerat 
1. INTRODUCTION 

The LHC discovery potential for a Z' in the reaction pp ^ Z' ^ with i = e, ^ is well 
known. As shown in previous studies, a Z' with a mass up to 5 TeV could be discovered at 
the LHC with 100 fb^^ We make here a summary of the detailed work described in Ref. [181] 
showing how, after a Z' signal has been detected at the LHC, one could identify it. In contrast 
to previous studies, where the models were either analyzed from a more theoretical point of 
view or a particular model was analyzed within a certain experimental frame, we combine here 
different experimental observables in order to investigate the realistic potential of the LHC 
experiments to distinguish between models and determine their parameters. In this study, two 
classes of Z' models are considered: Eq models, parametrized with cos /? and left-right (LR) 
models, parametrized with a^R (see Ref. [182], that we will follow, for a theoretical account). 



2. OBSERVABLES SENSITIVE TO Z' PROPERTIES 

Future measurements of Z' properties at the LHC can use the following observables: 

The total decay width of the Z' which is obtained from a fit to the invariant mass dis- 
tribution of the reconstructed dilepton system using a non-relativistic Breit-Wigner function: 

ao/[(M2 - M|,)^ + ai] with ai = r|,M|,. 

The Z' cross section times leptonic branching ratio which is calculated from the number 
of reconstructed dilepton events lying within ±3r around the observed peak. 

The leptonic forward-backward asymmetry Ap^, which is defined from the lepton angular 
distribution with respect to the quark direction in the center of mass frame, as: 

^hi + cos^ e*) + A^FB cos e* (1) 

d cos u* 8 

Ap^ can be determined with an unbinned maximum likelihood fit to the cos 9* distribution. 
Unfortunately, as the original quark direction in a proton-proton collider is not known, Ap^ 
cannot be used directly. However, it can be extracted from the kinematics of the dilepton system, 
as was shown in detail in [183]. The initial quark direction is assumed to be the boost direction 
of the i£ system with respect to the beam axis. The probability to assign the correct quark 
direction increases for larger rapidities of the dilepton system. A purer, though smaller signal 
sample, can thus be obtained by introducing a rapidity cut. For the following study we will 
require \ Yu\ > 0.8. 

The Z' rapidity distribution: To complete the Z' analysis, one can obtain some informa- 
tion about the fraction of Z"s produced from uu and ddhy analyzing the Z' rapidity distribution. 
Assuming that the and Z rapidity distribution has been measured in detail, following the 
ideas given in [184], relative parton distribution functions for u and d quarks as well as for 
the corresponding sea quarks and antiquarks are well known. Thus, the rapidity spectra can be 
calculated separately for uu and dd as well as for sea quark anti-quark annihilation and for the 
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mass region of interest to analyze the Z' rapidity distribution [185]. Using these distributions 
a fit can be performed to the Z' rapidity distribution which allows to obtain the corresponding 
fractions of Z"s produced from uu, dd as well as for sea quark anti-quark annihilation. This 
will thus reveal how the Z' couples to different quark flavors in a particular model. 

3. DISTINCTION BETWEEN MODELS AND PARAMETER DETERMINATION 

In the present analysis, PYTHIA events of the type pp 7, Z, Z' ee, fxfi were simulated at a 
center of mass energy of 14 TeV, and for different Z' models. These events were analyzed using 
simple acceptance cuts following the design criteria of ATLAS and CMS. The SM background 
relative to the signal cross section is found to be essentially negligible for the considered Z' 
models. We thus reconfirm the known Z' boson LHC discovery potential, to reach masses up 
to about 5 TeV for a luminosity of 100 fb'^ [182]. 

Let us discuss how well one can distinguish experimentally the different Z' models using 
the observables as defined before: aff x F, Ap^ as well as Run as obtained from the rapidity 
distribution. As a working hypothesis, a luminosity of 100 fb~^ and a Z' mass of 1.5 TeV will 
be assumed in the following. 

A precise knowledge of the cross section and the total width allows to make a first good 
distinction between the different models as we will be discussed later. It is not obvious how 
accurate cross sections can be measured and interpreted at the LHC. Following the procedure 
outlined in [184], an accuracy of ±1% could be assumed. It is however necessary to consider 
the other observables. 

Very distinct forward backward charge asymmetries are expected as a function of the 
dilepton mass and for the different Z' models, as shown in Figure[Tt. One finds that additional 
and complementary informations are also obtained from Ap^ measured in the interference re- 
gion. The Z' rapidity distribution is also analyzed. Figure [IJ5 shows the expected rapidity 
distribution for the Z'^ model. A particular Z' rapidity distribution is fitted using a linear com- 
bination of the three pure quark-antiquark rapidity distributions. The fit output gives the uu, dd 
and sea quarks fraction in the sample. In order to demonstrate the analysis power of this method 
we also show the Z'^ rapidity distribution which has equal couplings to uu and quarks. 

In a next step, assuming that a particular model has been selected, one would like to 
know how well the parameter(s), like cos /5 or aiR can be constrained. Figure |2l shows how 
the previously defined observables vary as the model parameters are varied. In the case of the 
Eq model for instance, one finds that cos /5 can not always be determined unambiguously. Very 
similar results can be expected for different observables but using very different values for cos p. 
Obviously, the combination of the various measurements, helps to reduce some ambiguities. 

If the Z' mass is increased, the number of events decreases drastically and the differences 
between the models start to become covered within the statistical fluctuations. For the assumed 
luminosity of 100 fb^^, one could still distinguish a Z'^ from a over a large parameter range 
and the Ap^ measurements provide some statistical significance up to Mz' = 2 — 2.5 TeV. On 
the contrary, a Z'^ could be differentiated from a Z'^, only up to a Z' mass of at most 2 TeV, as 
in that case, the dependence of Ap^ on the Z' mass is almost identical for the two models. 

4. CONCLUSIONS 

A realistic simulation of the study of the properties of Z' bosons in Eq and LR models has been 
performed for the LHC. We have shown that, in addition to the Z' production cross section times 
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Fi gure 1; Apg (a) as a function of Mu for four Z' models. The rapidity of the dilepton system is required to be 
larger than 0.8. The observable rapidity distribution for two Z' models is shown in (b), including the fit results 
which determine the types of qq fractions. A simulation of the statistical errors, including random fluctuations of 
the Z'^ model and with errors corresponding to a luminosity of 100 fb~^ has been included in both plots. 

total decay width, the measurement of the forward-backward lepton charge asymmetry, both on 
the Z' peak and in the interference region, provide complementary information. We have also 
shown that a fit of the rapidity distribution can provide a sensitivity to the Z' couplings to 
up-type and down-type quarks. The combination of all these observables would allow us to 
discriminate between Z' bosons of different models or classes of models for masses up to 2-2.5 
TeV, if a luminosity of 100 fb^^ is collected. 
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